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Abstract 

These lecture notes, based on a course given at the Zurich Clay Sum- 
mer School (June 23-July 18 2008), review our current mathematical 
understanding of the global behaviour of waves on black hole exterior 
backgrounds. Interest in this problem stems from its relationship to the 
non-linear stability of the black hole spacetimes themselves as solutions 
to the Einstein equations, one of the central open problems of general rel- 
ativity. After an introductory discussion of the Schwarzschild geometry 
and the black hole concept, the classical theorem of Kay and Wald on 
the boundedness of scalar waves on the exterior region of Schwarzschild 
is reviewed. The original proof is presented, followed by a new more ro- 
bust proof of a stronger boundedness statement. The problem of decay of 
scalar waves on Schwarzschild is then addressed, and a theorem proving 
quantitative decay is stated and its proof sketched. This decay statement 
is carefully contrasted with the type of statements derived heuristically in 
the physics literature for the asymptotic tails of individual spherical har- 
monics. Following this, our recent proof of the boundedness of solutions 
to the wave equation on axisymmetric stationary backgrounds (including 
slowly-rotating Kerr and Kerr-Newman) is reviewed and a new decay re- 
sult for slowly-rotating Kerr spacetimes is stated and proved. This last 
result was announced at the summer school and appears in print here 
for the first time. A discussion of the analogue of these problems for 
spacetimes with a positive cosmological constant A > follows. Finally, 
a general framework is given for capturing the red-shift effect for non- 
extremal black holes. This unifies and extends some of the analysis of the 
previous sections. The notes end with a collection of open problems. 
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1 Introduction: General relativity and evolu- 
tion 

Black holes are one of the fundamental predictions of general relativity. At the 
same time, they are one of its least understood (and most often misunderstood) 
aspects. These lectures intend to introduce the black hole concept and the 
analysis of waves on black hole backgrounds {M.,g) by means of the example 
of the scalar wave equation 

□sV' = 0. (1) 
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We do not assume the reader is familiar with general relativity, only ba- 
sic analysis and differential geometry. In this introductory section, we briefly 
describe general relativity in outline form, taking from the beginning the evolu- 
tionary point of view which puts the Cauchy problem for the Einstein equations- 
the system of nonlinear partial differential equations (see ([2]) below) governing 
the theory-at the centre. The problem ([T]) can be viewed as a poor man's lin- 
earisation for the Einstein equations. Study of ([T|) is then intimately related 
to the problem of the dynamic stability of the black hole spacetimes {^4,g) 
themselves. Thus, one should view the subject of these lectures as intimately 
connected to the very tenability of the black hole concept in the theory. 

1.1 General relativity and the Einstein equations 

General relativity postulates a 4-dimensional Lorentzian manifold {M, g)-space- 
time-which is to satisfy the Einstein equations 



Here, Rp,i,, R denote the Ricci and scalar curvature of g, respectively, and T^^ 
denotes a symmetric 2-tensor on M termed the stress-energy-momentum ten- 
sor of matter. (Necessary background on Lorentzian geometry to understand 
the above notation is given in Appendix [A]) The equations {[21) in of them- 
selves do not close, but must be coupled to "matter equations" satisfied by 
a collection {^'i} of matter fields defined on Ai, together with a constitutive 
relation determining T^^ from {g,'^i}. These equations and relations arc stipu- 
lated by the relevant continuum field theory (clcctromagnctism, fluid dynamics, 
etc.) describing the matter. The formulation of general relativity represents the 
culmination of the classical field-theoretic world-view where physics is governed 
by a closed system of partial differential equations. 

Einstein was led to the system in 1915, after a 7-year struggle to in- 
corporate gravity into his earlier principle of relativity. In the field-theoretic 
formulation of the "Newtonian" theory, gravity was described by the Newto- 
nian potential satisfying the Poisson equation 



where fi denotes the mass-density of matter. It is truly remarkable that the 
constraints of consistency were so rigid that encorporating gravitation required 
finally a complete reworking of the principle of relativity, leading to a theory 
where Newtonian gravity, special relativity and Euclidean geometry each emerge 
as limiting aspects of one dynamic geometrical structure-the Lorentzian mctric- 
naturally living on a 4-dimensional spacctimc continuum. A second remarkable 
aspect of general relativity is that, in contrast to its Newtonian predecessor, the 
theory is non-trivial even in the absence of matter. In that case, we set T^;^ = 
and the system takes the form 



(2) 



A = 47r/i, 



(3) 



R^v = 0. 



(4) 



5 



The equations are known as the Einstein vacuum equations. Whereas ([3]) is 
a hnear eUiptic equation, ([4]) can be seen to form a closed system of non-Unear 
(but quasilinear) wave equations. Essentially all of the characteristic features 
of the dynamics of the Einstein equations arc already present in the study of 
the vacuum equations ([^ . 

1.2 Special solutions: Minkowski, Schwarzschild, Kerr 

To understand a theory like general relativity where the fundamental equations 
Q are nonlinear, the first goal often is to identify and study important explicit 
solutions, i.e., solutions which can be written in closed formfj Much of the early 
history of general relativity centred around the discovery and interpretation of 
such solutions. The simplest explicit solution to the Einstein vacuum equations 
(|4]) is Minkowski space M'^^^. The next simplest solution of ([4|) is the so-called 
Schwarzschild solution, written down |139| already in 1916. This is in fact a 
one-parameter family of solutions (A^, gj\/), the parameter M identified with 
mass. See ([S]) below for the metric form. The Schwarzschild family lives as 
a subfamily in a larger two-parameter family of explicit solutions (A^,5A/,a) 
known as the Kerr solutions, discussed in Section 15.11 These were discovered 
only much later ^ (1963). 

When the Schwarzschild solution was first written down in local coordinates, 
the necessary concepts to understand its geometry had not yet been developed. 
It took nearly 50 years from the time when Schwarzschild was first discovered for 
its global geometry to be sufficiently well understood so as to be given a suitable 
name: Schwarzschild and Kerr were examples of what came to be known as black 
hole spacetimefH. The Schwarzschild solution also illustrates another feature of 
the Einstein equations, namely, the presence of singularities. 

We will spend Section [5] telling the story of the emergence of the black hole 
notion and sorting out what the distinct notions of "black hole" and "singular- 
ity" mean. For the purpose of the present introductory section, let us take the 
notion of "black hole" as a "black box" and make some general remarks on the 
role of explicit solutions, whatever might be their properties. These remarks 
are relevant for any physical theory governed by an evolution equation. 

1.3 Dynamics and the stability problem 

Explicit solutions arc indeed suggestive as to how general solutions behave, but 
only if they are appropriately "stable" . In general relativity, this notion can 
in turn only be understood after the problem of dynamics for 1^ has been 
formulated, that is to say, the Cauchy problem. 

In contrast to other non-linear field theories arising in physics, in the case 
of general relativity, even formulating the Cauchy problem requires addressing 
several conceptual issues (e.g. in what sense is ([4|) hyperbolic?), and these took a 
long time to be correctly sorted out. Important advances in this process include 

^The traditional terminology in general relativity for such solutions is exact solutions. 
■^Tliis name is due to John Wheeler. 
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the identification of the harmonic gauge by de Bonder [70], the existence and 
uniqueness theorems for general quasihnear wave equations in the 1930's based 
on work of Friedrichs, Schauder, Sobolev, Pctrovsky, Leray and others, and 
Leray's notion of global hypcrbolicity [112j . The well-poscdncss of the appro- 
priate Cauchy problem for the vacuum equations dH) was finally formulated and 
proven in celebrated work of Choquet-Bruhat [32 (1952) and Choquet-Bruhat- 
Geroch [35] (1969). See Appendix IB] for a concise survey of these developments 
and the precise statement of the existence and uniqueness theorems and some 
comments on their proof. 

In retrospect, much of the confusion in early discussions of the Schwarzschild 
solution can be traced to the lack of a dynamic framework to understand the 
theory. It is only in the context of the language provided by [35j that one can 
then formulate the dynamical stability problem and examine the relevance of 
various explicit solutions. 

The stability of Minkowski space was first proven in the monumental work 
of Christodoulou and Klaincrman [5T] . See Appendix IB. 51 for a formulation of 
this result. The dynamical stability of the Kerr family as a family of solutions 
to the Cauchy problem for the Einstein equations, even restricted to parameter 
values near Schwarzschild, i.e. \a\ « M|3 is yet to be understood and poses 
an important challenge for the mathematical study of general relativity in the 
coming years. See Section lS^ for a formulation of this problem. In fact, even the 
most basic linear properties of waves (e.g. solutions of ([TJ) on Kerr spacetime 
backgrounds (or more generally, backgrounds near Kerr) have only recently 
been understood. In view of the wave-like features of the Einstein equations (|4]) 
(see in particular Appendix IB. 4p . this latter problem should be thought of as a 
prerequisite for understanding the non-linear stability problem. 

1.4 Outline of the lectures 

The above linear problem will be the main topic of these lectures: We shall here 
develop from the beginning the study of the linear homogeneous wave equa- 
tion dl]) on fixed black hole spacetime backgrounds {Ai,g). We have already 
referred in passing to the content of some of the later sections. Let us give 
here a complete outline: Section |2] will introduce the black hole concept and 
the Schwarzschild geometry in the wider context of open problems in general 
relativity. Section [3] will concern the basic boundedncss properties for solutions 
tp of ([T]) on Schwarzschild exterior backgrounds. Section [4] will concern quan- 
titative decay properties for ip. Section [5] will move on to spacctimes {Ai,g) 
"near" Schwarzschild, including slowly rotating Kerr, discussing boundedness 
and decay properties for solutions to ([T]) on such (A^,g), and ending in Sec- 
tion l5.6l with a formulation of the non-linear stability problem for Kerr, the open 
problem which in some sense provides the central motivation for these notes. 
Section [6] will consider the analogues of these problems in spacetimes with a 

■^Notc that without symmetry assumptions one cannot study the stability problem for 
Schwarzschild per se. Only the larger Kerr family can be stable. 
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positive cosmological constant A, Section [7] will give a multiplier- type estimate 
valid for general non-degenerate Killing horizons which quantifies the classical 
red-shift effect. The importance of the red-shift effect as a stabilising mechanism 
for the analysis of waves on black hole backgrounds will be a common theme 
throughout these lectures. The notes end with a collection of open problems in 
Section [51 

The proof of Theorem 15.21 of Section [5] as well as all results of Section [7] 
appear in print in these notes for the first time. The discussion of Section 13.31 
as well as the proof of Theorem 14.11 have also been streamlined in comparision 
with previous presentations. We have given a guide to background literature in 
Sections 15:1 Ol [531 and [Ql 

We have tried to strike a balance in these notes between making the dis- 
cussion self-containted and providing the necessary background to appreciate 
the place of the problem in the context of the current state of the art of 
the Cauchy problem for the Einstein equations or and the main open 
problems and conjectures which will guide this subject in the future. Our so- 
lution has been to use the history of the Schwarzschild solution as a starting 
point in Section [2l for a number of digressions into the study of gravitational 
collapse, singularities, and the weak and strong cosmic censorship conjectures, 
deferring, however, formal development of various important notions relating 
to Lorentzian geometry and the well-posedness of the Einstein equations to a 
series of Appendices. We have already referred to these appendices in the text. 
The informal nature of Section [2l should make it clear that the discussion is 
not intended as a proper survey, but merely to expose the reader to important 
open problems in the field and point to some references for further study. The 
impatient reader is encouraged to move quickly through Section[2lat a first read- 
ing. The problem is itself rather self-contained, requiring only basic analysis 
and differential geometry, together with a good understanding of the black hole 
spacetimcs, in particular, their so-called causal geometry. The discussion of 
Section [21 should be more than enough for the latter, although the reader may 
want to supplement this with a more general discussion, for instance |55] . 

These notes accompanied a series of lectures at a summer school on "Evolu- 
tion Equations" organized by the Clay Mathematics Institute, June-July 2008. 
The centrality of the evolutionary point of view in general relativity is often 
absent from textbook discussions. (See however the recent |133| .) We hope that 
these notes contribute to the point of view that puts general relativity at the 
centre of modern developments in partial differential equations of evolution. 

2 The Schwarzschild metric and black holes 

Practically all concepts in the development of general relativity and much of its 
history can be told from the point of view of the Schwarzschild solution. We 
now readily associate this solution with the black hole concept. It is important 
to remember, however, that the Schwarzschild solution was first discovered in a 
thoroughly classical astrophysical setting: it was to represent the vacuum region 



8 



outside a star. The black hole interpretation-though in some sense inevitable- 
historically only emerged much later. 

The most efficient way to present the Schwarzschild solution is to begin at 
the onset with Kruskal's maximal extension as a point of departure. Instead, 
we shall take advantage of the informal nature of the present notes to attempt 
a more conversational and "historical" presentation of the Schwarzschild metric 
and its interpretation^ Although certainly not the quickest route, this approach 
has the advantage of highlighting the themes which have become so important 
in the subject-in particular, singularities, black holes and their event horizons- 
with the excitement of their step-by-step unravelling from their origin in a model 
for the simplest of general relativistic stars. The Schwarzschild solution will 
naturally lead to discussions of the Oppenheimer-Snyder collapse model, the 
cosmic censorship conjectures, trapped surfaces and Penrose's incompleteness 
theorems, and recent work of Christodoulou on trapped surface formation in 
vacuum collapse, and we elaborate on these topics in Sections 12.6112.81 (The 
discussion in these three last sections was not included in the lectures, however, 
and is not necessary for understanding the rest of the notes.) 

2.1 Schwarzschild's stars 

The most basic self-gravitating objects are stars. In the most primitive stellar 
models, dating from the 19th century, stars are modeled by a self-gravitating 
fluid surrounded by vacuum. Moreover, to a first approximation, classically 
stars are spherically symmetric and static. 

It should not be surprising then that early research on the Einstein equations 
(l2|) would address the question of the existence and structure of general rela- 
tivistic stars in the new theory. In view of our above discussion, the most basic 
problem is to understand spherically symmetric, static metrics, represented in 
coordinates (i,r, 0,0), such that the spacetimc has two regions: In the region 
r < i?o-the interior of the star-the metric should solve a suitable Einstein-matter 
system ([2]) with appropriate matter, and in the region r > i?o^the exterior of 

*This in no way should be considered as a true attempt at the history of the solution, 
simply a pedagogical approach to its study. See for example |76| . 
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the star-the spacetime should be vacuum, i.e. the metric should solve (2]). 
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This is the problem first addressed by Schwarzschild [1391 1140| , already in 
1916. Schwarzschild considered the vacuum region first [139] and arriveqfl at 
the one-parameter family of solutions: 

1 _ ^ j dt^ + (l - dr' + r\de^ + sin2 (5) 

Every student of this subject should explicitly check that this solves ([3]) (Ex- 
ercise). 

In |140j . Schwarzschild found interior metrics for the darker shaded region 
r < Rq above. In this region, matter is described by a perfect fluid. Wc shall not 
write down explicitly such metrics here, as this would require a long digression 
into fluids, their equations of state, etc. See [33]. Suffice it to say here that the 
existence of such solutions required that one take the constant M positive, and 
the value Rq marking the boundary of the star always satisfied Rq > 2M. The 
constant M could then be identified with the total mass of the star as measured 
by considering the orbits of far-away test particles In fact, for most reasonable 
matter models, static solutions of the type described above only exist under a 
stronger restriction on Rq (namely Rq > 9A'//4) now known as the Buchdahl 
inequality. See [Ml [2 [97]. 

The restriction on i?o necessary for the existence of Schwarzschild's stars 
appears quite fortuitous: It is manifest from the form {[5]) that the components 
of g are singular if the (t, r, 0, 0) coordinate system for the vacuum region is 
extended to r = 2M. But a natural (if perhaps seemingly of only academic 
interest) question arises, namely, what happens if one does away completely 

^As is often the case, the actual history is more comphcatcd. Schwarzschild based his 
work on an earlier version of Einstein's theory which, while obtaining the correct vacuum 
equations, imposed a condition on admissible coordinate systems which would in fact exclude 
the coordinates of ([5|). Thus he had to use a rescaled r as a coordinate. Once this condition 
was removed from the theory, there is no reason not to take r itself as the coordinate. It is in 
this sense that these coordinates can reasonably be called "Schwarzschild coordinates". 

^Test particles in general relativity follow timelike geodesies of the spacetime metric. Ex- 
ercise: Explain the statement claimed about far-away test particles. See also ApDendix lB.2.3l 
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with the star and tries simply to consider the expression ([SD for all values of 
r? This at first glance would appear to be the problem of understanding the 
gravitational field of a "point particle" with the particle removedQ 

For much of the history of general relativity, the degeneration of the met- 
ric functions at r = 2M, when written in these coordinates, was understood as 
meaning that the gravitational field should be considered singular there. This 
was the famous Schwarzschild "singularity" @ Since "singularities" were consid- 
ered "bad" by most pioneers of the theory, various arguments were concocted 
to show that the behaviour of g where r = 2M is to be thought of as "patholog- 
ical", "unstable", "unphysical" and thus, the solution should not be considered 
there. The constraint on i?o related to the Buclidahl inequality seemed to give 
support to this point of view. See also [75] . 

With the benefit of hindsight, we now know that the interpretation of the 
previous paragraph is incorrect, on essentially every level: neither is r = 2M 
a singularity, nor are singularities-which do in fact occurl-necessarily to be 
discarded! Nor is it true that non-existence of static stars renders the behaviour 
at r = 2M-whatever it is- "unstable" or "unphysical"; on the contrary, it was 
an early hint of gravitational collapse! Let us put aside this hindsight for now 
and try to discover for ourselves the geometry and "true" singularities hidden 
in dS]), as well as the correct framework for identifying "physical" solutions. In 
so doing, we are retracing in part the steps of early pioneers who studied these 
issues without the benefit of the global geometric framework we now have at 
our disposal. All the notions referred to above will reveal themselves in the next 
subsections. 

2.2 Extensions beyond the horizon 

The fact that the behaviour of the metric at r = 2M is not singular, but simply 
akin to the well-known breakdown of the coordinates ([5]) at 6* = 0,7r (this latter 
breakdown having never confused anyone. . . ), is actually quite easy to see, and 
there is no better way to appreciate this than by doing the actual calculations. 
Let us see how to proceed. 

First of all, before even attempting a change of coordinates, the following is 
already suggestive: Consider say a futurc-directecfl ingoing radial null geodesic. 

''Hence the title of [139] . 

*Let the reader keep in mind that there is a good reason for the quotation marks here and 
for those that follow. 

^We time-orient the metric by dt- See Appendix lAl 
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The image of such a null ray is in fact depicted below: 



II 

Y 

One can compute that this has finite afRne length to the future, i.e. these 
null geodesies are future-incomplete, while scalar curvature invariants remain 
bounded as s ^ oo. It is an amusing exercise to put oneself in this point of 
view and carry out the above computations in these coordinates. 

Of course, as such the above doesn't show anything^ But it turns out that 
indeed the metric can be extended to be defined on a "bigger" manifold. One 
defines a new coordinate 

r = t + 2Mlog(r-2A/). 
This metric then takes the form 

g = _ _ ^ j (dt^f + Mirfr dr + (l + ^) dr^ + r^das2 (6) 

on r > 2M. Note that ^ = each interpreted in its respective coordinate 
system. But now ([6]) can clearly be defined in the region r > 0, -oo < t* < oo, 
and, by explicit computation or better, by analytic continuation, the metric ([5]) 
must satisfy (0]) for all r > 0. 

Transformations similar to the above were already known to Eddington and 
Lemaitre jlllj in the early 1930's. Nonetheless, from the point of view of that 
time, it was difficult to interpret their significance. The formalisation of the 
manifold concept and associated language had not yet become common knowl- 
edge to physicists (or most mathematicians for that matter), and in any case, 
there was no selection principle as to what should the underlying manifold M 
be on which a solution g to dH) should live, or, to put it another way, the domain 
of g in (|4|) is not specified a priori by the theory. So, even if the solutions ([6]) 
exist, how do we know that they are "physical"? 

This problem can in fact only be clarifed in the context of the Cauchy problem 
for (l2|) coupled to appropriate matter. Once the Cauchy problem for (|4]) is 
formulated correctly, then one can assign a unique spacetime to an appropriate 
notion of initial data set. This is the maximal development of Appendix [b1 
It is only the initial data set, and the matter model, which can be judged for 

^"Consider for instance a cone with the vortex removed. . . 
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"physicality" . One cannot throw away the resulting maximal development just 
because one does not like its properties! 

From this point of view, the question of whether the extension ([S]) was 
"physical" was resolved in 1939 by Oppcnheimer and Snyder |125| . Specifically, 
they showed that the extension for t >0 arose as a subset of the solution to 
the Einstein equations coupled to a reasonable (to a first approximation at least) 
matter model, evolving from physically plausible initial data. With hindsight, 
the notion of black hole was born in that paper. 

Had history proceeded differently, we could base our futher discussion on |125| . 
Unfortunately, the model [125| was ahead of its time. As mentioned in the in- 
troduction, the proper language to formulate the Cauchy problem in general 
only came in 1969 |35| . The interpretation of explicit solutions remained the 
main route to understanding the theory. We will follow thus this route to the 
black hole concept-via the geometric study of so-called maximally extended 
Schwarzschild-even though this spacetime is not to be regarded as "physical" . 
It was through the study of this spacetime that the relevant notions were first 
understood and the important Penrose diagrammatic notation was developed. 
We shall return to [125] only in Section [2.5.3l 

2.3 The maximal extension of Synge and Kruskal 

Let us for now avoid the question of what the underlying manifold "should" be, 
a question whose answer requires physical input (sec paragraphs above), and 
simply ask the purely mathematical question of how big the underlying manifold 
"can" be. This leads to the notion of a "maximally extended" solution. In the 
case of Schwarzschild, this will be a spacetime which, although not to be taken 
as a model for anything per se, can serve as a reference for the formulation of 
all important concepts in the subject. 

To motivate this notion of "maximally extended" solution, let us examine 
our first extension a little more closely. The light cones can be drawn as follows: 

A A 




I IM 
I II 

Let us look say at null geodesies. One can see (Exercise) that future di- 
rected null geodesies either approach r = or are future-complete. In the for- 
mer case, scalar invariants of the curvature blow up in the limit as the affine 
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parameter approaches its supremum (Exercise). The spacetime is thus "sin- 
gular" in this sense. It thus follows from the above properties that the above 
spacetime is future null geodesically incomplete, but also future null geodesically 
inextendible as a Lorentzian metric, i.e. there does not exist a larger 4- 
dimensional Lorentzian manifold with metric such that the spacetime above 
embeds isometrically into the larger one such that a future null geodesic passes 
into the extension. 

On the other hand, one can see that past-directed null geodesies are not all 
complete, yet no curvature quantity blows up along them (Exercise). Again, 
this suggests that something may still be missing! 

Synge was the first to consider these issues systematically and construct 
"maximal extensions" of the original Schwarzschild metric in a paper |146| of 
1950. A more concise approach to such a construction was given in a celebrated 
1960 paper [107j of Kruskal. Indeed, let ^4 be the manifold with differentiable 
structure given by W x §^ where U is the open subset - i?^ < 1 of the (T, R)- 
plane. Consider the metric g 

r 

where r is defined implicitly by 




The region U is depicted below: 




This is a spherically symmetric 4-dimcnsional Lorentzian manifold satisfying 
dl]) such that the original Schwarzschild metric is isometric to the region R > \T\ 
(where t is given by tanh(^) = T/R), and our previous partial extension is 
isometric to the region T > -R (Exercise). It can be shown now (Exercise) 
that {A4,g) is inextendible as a (in fact C^) Lorentzian manifold, that is to 
say, if 

i:{M,g)^{M,g) 

is an isometric embedding, where {M,g) is a (in fact C") 4-dimensional 
Lorentzian manifold, then necessarily z(A^) = ^4. 
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The above property defines the sense in which our spacetime is "maximally" 
extended, and thus, {M,g) is called sometimes maximally- extended Schwarz- 
schild. In later sections, we will often just call it "the Schwarzschild solution" . 

Note that the form of the metric is such that the light cones are as depicted. 
Thus, one can read off much of the causal structure by sight. 

It may come as a surprise that in maximally-extended Schwarzschild, there 
are two regions which are isometric to the original r > 2M Schwarzschild region. 
Alternatively, a Cauchy surfac^Hl will have topology §^ x M with two asymptoti- 
cally flat ends. This suggests that this spacetime is not to be taken as a physical 
model. We will discuss this later on. For now, let us simply try to understand 
better the global geometry of the metric. 

2.4 The Penrose diagram of Schwarzschild 

There is an even more useful way to represent the above spacetime. First, let 
us define null coordinates U = T- R, V = T + R. These coordinates have infinite 
range. We may rcscale them by u = u(U), v = v{V) to have finite range. (Note 
the freedom in the choice of u and vl) The domain of (it,u) coordinates, when 
represented in the plane where the axes are at 45 and 135 degrees with the 
horizontal, is known as a Penrose diagram of Schwarzschild. Such a Penrose 
diagram is depicted belo-wF^: 




In more geometric language, one says that a Penrose diagram corresponds 
to the image of a bounded conformal map 

A^/S0(3) = Q^Ri+i, 

where one makes the identification v = t + x, u = t - x where (t, x) are now the 
standard coordinates M}^^ represented in the standard way on the plane. We 
further assume that the map preserves the time orientation, where Minkowski 
space is oriented by dt ■ (In our application, this is a fancy way of saying that 
u'{U),v'{V) > 0). It follows that the map preserves the causal structure of 

^^See Appendix IAI 

^^How can (n,f) be chosen so that the r = boundaries are horizontal lines? (Exercise) 
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Q. In particular, we can "read off" the radial null geodesies of Ai from the 
depiction. 

Now wc may turn to the boundary induced by the causal embedding. We 
define X"^ to be the boundary components as depicted0 These are characterized 
geometrically as follows: are limit points of future-directed null rays in Q 
along which r ^ oo. Similarly. X" are limit points of past-directed null rays for 
which r ->■ oo. We call future null infinity and X" past null infinity. The 
remaining boundary components i'^ and i"^ depicted are often given the names 
spacelike infinity and future (past) timelike infinity, respectively. 

In the physical application, it is important to remember that asymptoti- 
cally flacj spacetimes like our {M.,g) are not meant to represent the whole 
universq^, but rather, the gravitational field in the vicinity of an isolated self- 
gravitating system. X"*" is an idealization of far away observers who can receive 
radiation from the system. In this sense, "we"~as astrophysical observers of 
stellar collapse, say-are located at X^. The ambient causal structure of M}^^ 
allows us to talk about J" (p) n Q for p e xj^ and this will lead us to the black 
hole concept. Therein lies the use of the Penrose diagram representation. 

The systematic use of the conformal point of view to represent the global ge- 
ometry of spacetimes is one of the many great contributions of Penrose to general 
relativity. These representations can be traced back to the well-known "space- 
time diagrams" of special relativity, promoted especially by Synge |147j . The 
"formal" use of Penrose diagrams in the sense above goes back to Carter [55] , in 
whose hands these diagrams became a powerful tool for determining the global 
structure of all classical black hole spacetimes. It is hard to overemphasize how 
important it is for the student of this subject to become comfortable with these 
representations. 

2.5 The black hole concept 

With Penrose diagram notation, we may now explain the black hole concept. 
2.5.1 The definitions for Schwarzschild 

First an important remark: In Schwarzschild, the boundary component X^ en- 
joys a limiting affine completeness. More specifically, normalising a sequence 
of ingoing radial null vectors by parallel transport along an outgoing geodesic 
meeting X^ , the affine length of the null geodesies generated by these vectors, 

^^Our convention is that open endpoint circles are not contained in the intervals they bound, 
and dotted lines are not contained in the regions they bound, whereas solid lines arc. 
^*See Appendix IB . 2 . 3] for a definition. 

^^The study of that problem is what is known as "cosmology". See Section [6] 
l^Refer to Appendix |X] for J*. 



16 



parametrized by their parallel transport (restricted to J (I"*")), tends to infinity: 



/K 

d 

This has the interpretation that far-away observers in the radiation zone can 
observe for all time. (This is in some sense related to the presence of timelike 
geodesies near infinity of infinite length, but the completeness is best formulated 
with respect to X'*^.) A similar statement clearly holds for X". 

Given this completeness property, we define now the black hole region to 
be Q \ J"(X^), and the white hole region to be Q \ J^(X"). Thus, the black 
hole corresponds to those points of spacetime which cannot "send signals" to 
future null infinity, or, in the physical interpetation, to far-away observers who 
(in view of the completeness property!) nonetheless can observe radiation for 
infinite time. 

The future boundary of J" (X^ ) in Q (alternatively characterized as the past 
boundary of the black hole region) is a null hypersurfacc known as the future 
event horizon, and is denoted by 7^"*". Exchanging past and future, we obtain 
the past event horizon "HT . In maximal Schwarzschild, {r = 2M} = 'W' u "HT . 
The subset J^(X^) n J^(X") is known as the domain of outer communications. 

2.5.2 Minkowski space 

Note that in the case of Minkowski space, Q = IR'^^^/S0(3) is a manifold with 
boundary since the S0(3) action has a locus of fixed points, the centre of sym- 
metry. A Penrose diagram of Minkowski space is easily seen to be: 

o 

O 
II 

A 

a 

Here X^ and X^ are characterized as before, and enjoy the same completeness 
property as in Schwarzschild. One reads off immediately that J"(X^) n Q = Q, 
i.e. M?^^ does not contain a black hole under the above definitions. 
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2.5.3 Oppenheimer-Snyder 

Having now the notation of Penrose diagrams, we can concisely describe the ge- 
ometry of the Oppenheimer-Snyder solutions referred to earlier, without giving 
explicit forms of the metric. Like Schwarzschild's original picture of the grav- 
itational field of a spherically symmetric star, these solutions involve a region 
r < Ro solving ([2]) and r > Rq satisfying (|4]) . The matter is described now by a 
pressurelcss fluid which is initially assumed homogeneous in addition to being 
spherically symmetric. The assumption of staticity is however dropped, and 
for appropriate initial conditions, it follows that Ro{t*) with respect to a 
suitable time coordinate t* . (In fact, the Einstein equations can be reduced to 
an o.d.e. for Ro{t*).) We say that the star "collapses" A Penrose diagram 
of such a solution (to the future of a Cauchy hypcrsurface) can be seen to be of 
the form: 




The lighter shaded region is isometric to a subset of maximal Schwarzschild, 
in fact a subset of the original extension of Section 12.21 In particular, the 
completeness property of holds, and as before, we identify the black hole 
region to be Q\ J~{X^). 

In contrast to maximal Schwarzschild, where the initial configuration is un- 
physical (the Cauchy surface has two ends and topology RxS^), here the initial 
configuration is entirely plausible: the Cauchy surface is topologically , and its 
geometry is not far from Euclidean space. The Oppenheimer-Snyder model [125j 
should be viewed as the most basic black hole solution arising from physically 
plausible regular initial dataF^ 

It is traditional in general relativity to "think" Oppenheimer-Snyder but 
"write" maximally-extended Schwarzschild. In particular, one often imports 
terminology like "collapse" in discussing Schwarzschild, and one often reformu- 
lates our definitions replacing with one of its connected components, that is 
to say, we will often write J~{I*) n J^(I+) meaning J^{I\) n J"^(X^), etc. In 
any case, the precise relation between the two solutions should be clear from the 
above discussion. In view of Cauchy stability results [51], sufficiently general 
theorems about the Cauchy problem on maximal Schwarzschild lead immedi- 
ately to such results on Oppenheimer-Snyder. (See for instance the exercise in 
Section fS. 2. 61 ) One should always keep this relation in mind. 

^^Note that _Ro(t*) -> does not mean that the star collapses to "a point", merely that the 
spheres which foliate the interior of the star shrink to area. The limiting singular boundary 
is a spacelike hypcrsurface as depicted. 

Note however the end of Scction l2.6.2l 



18 



2.5.4 General definitions? 

The above definition of black hole for the Schwarzschild metric should be thought 
of as a blueprint for how to define the notion of black hole region in general. 
That is to say, to define the black hole region, one needs 

1 . some notion of future null infinity , 

2. a way of identifying J"(I^), and 

3. some characterization of the "completeness" of X^F^ 

If is indeed complete, we can define the black hole region as 

"the complement in M of J^{I^)" . 

For spherically symmetric spacetimes arising as solutions of the Cauchy prob- 
lem for ([2]), one can show that there always exists a Penrose diagram, and thus, a 
definition can be formalised along precisely these lines (see [HOI)- For spacetimes 
without symmetry, however, even defining the relevant asymptotic structure so 
that this structure is compatible with the theorems one is to prove is a main 
part of the problem. This has been accomplished definitively only in the case 
of perturbations of Minkowski space. In particular, Christodoulou and Klainer- 
man |51| have shown that spacetimes arising from perturbations of Minkowski 
initial data have a complete in a well defined sense, whose past can be iden- 
tified and is indeed the whole spacctime. See Appendix IB. 51 That is to say, 
small perturbations of Minkowski space cannot form black holes. 

2.6 Birkhoff's theorem 

Formal Penrose diagrams are a powerful tool for understanding the global causal 
structure of spherically symmetric spacetimes. Unfortunately, however, it turns 
out that the study of spherically symmetric vacuum spacetimes is not that rich. 
In fact, the Schwarzschild family parametrizes all spherically symmetric vacuum 
spacetimes in a sense to be explained in this section. 

2.6.1 Schvifarzschild for M <0 

Before stating the theorem, recall that in discussing Schwarzschild we have 
previously restricted to parameter value M > 0. For the uniqueness statement, 
we must enlarge the family to include all parameter values. 

^^The characterization of completeness can be formulated for general asymptotically flat 
vacuum space times using the results of |51j . This formulation is due to Christodoulou |47) . 
Previous attempts to formalise these notions rested on "asymptotic simplicity" and "weak 
asymptotic simplicity" . See 1911 . Although the qualitative picture suggested by these notions 
appears plausible, the detailed asymptotic behaviour of solutions to the Einstein equations 
turns out to be much more subtle, and Christodoulou has proven |48) that these notions 
cannot capture even the simplest generic physically interesting systems. 
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If we set M = in dSJ , we of course obtain Minkowski space in spherical polar 
coordinates. A suitable maximal extension is Minkowski space as we know it, 
represented by the Penrose diagram of Section 12.5.21 

On the other hand, we may also take M < in ([S]) . This is so-called negative 
mass Schwarzschild. The metric element dS]) for such M is now regular for all 
r > 0. The limiting singular behaviour of the metric at r = is in fact essential, 
i.e. one can show that along incxtcndible incomplete geodesies the curvature 
blows up. Thus, one immediately arrives at a maximally extended solution 
which can be seen to have Penrose diagram: 

o ' 

II I o 
t. 1 / 

A 

V 

Note that in contrast to the case of M'^^^, the boundary r = is here depicted 
by a dotted line denoting (according to our conventions) that it is not part of 
Q! 

2.6.2 Naked singularities and v^reak cosmic censorship 

The above spacetime is interpreted as having a "naked singularity" . The tra- 
ditional way of describing this in the physics literature is to remark that the 
"singularity" B = {r = 0} is "visible" to , i.e., J~{I'^) nB t 0. From the point 
of view of the Cauchy problem, however, this characterization is meaningless 
because the above maximal extension is not globally hyperbolic^i.e. it is not 
uniquely characterized by an appropriate notion of initial datacj From the 
point of view of the Cauchy problem, one must not consider maximal exten- 
sions but the maximal Cauchy development of initial data, which by definition 
is globally hyperbolic (see Theorem IB. 41 of Appendix IB|) . Considering an in- 
extendible spacelike hypersurface S as a Cauchy surface, the maximal Cauchy 

^"See Appendix |A] for the definition of global hypcrbolicity. 
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development of S would be the darker shaded region depicted below: 




The proper characterization of "having a naked singularity" , from the point of 
view of the darker shaded spacetime, is that its is incomplete. Of course, 
this example does not say anything about the dynamic formation of naked 
singularities, because the inital data hypersurface S is already in some sense 
"singular" , for instance, it is geodesically incomplete, and the curvature blows 
up along incomplete geodesies. The dynamic formation of a naked singularity 
from regular, complete initial data would be pictured by: 



q' 

o 




where we are to understand also in the above that is incomplete. The conjec- 
ture that for generic asympotically flal[^ initial data for "reasonable" Einstein- 
matter systems, the maximal Cauchy development "possesses a complete I^" 
is known as weak cosmic censorship}^ 

In light of the above conjecture, the story of the Oppenheimer-Snyder so- 
lution and its role in the emergence of the black hole concept docs have an 
interesting epilogue. Recall that in the Oppenheimer-Snyder solutions, the re- 
gion r < Rq, in addition to being spherically symmetric, is homogeneous. It 
turns out that by considering spherically symmetric initial data for which the 
"star" is no longer homogeneous, Christodoulou has proven that one can arrive 
at spacetimes for which "naked singularities" form [39j with Penrose diagram 
as above and with incomplete. Moreover, it is shown in [39j that this occurs 
for an open subset of initial data within spherical symmetry, with respect to a 
suitable topology on the set of spherically symmetric initial data. Thus, weak 
cosmic censorship is violated in this model, at least if the conjecture is restricted 
to spherically symmetric data. 

^^Sec Appendix IB . 2 . 3l for a formulation of this notion. Note that asymptotically flat data 
are in particular complete. 

■^■^This conjecture is originally due to Penrose |127l . The present formulation is taken from 
Christodoulou [47]. 
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The fact that in the Oppenheimer-Snyder solutions black holes formed ap- 
pears thus to be a rather fortuitous accident! Nonetheless, we should note that 
the failure of weak cosmic censorship in this context is believed to be due to the 
inappropriateness of the pressureless model, not as indicative of actual phenom- 
ena. Hence, the restriction on the matter model to be "reasonable" in the formu- 
lation of the conjecture. In a remarkable series of papers, Christodoulou [151 H7] 
has shown weak cosmic censorship to be true for the Einstein-scalar field system 
under spherical symmetry. On the other hand, he has also shown [53] that the 
assumption of genericity is still necessary by explicitly constructing solutions of 
this system with incomplete and Penrose diagram as depicted above 

2.6.3 Birkhoff's theorem 

Let us understand now by "Schwarzschild solution with parameter Af" (where 
M € M) the maximally extended Schwarzschild metrics described above. 
We have the so-called Birkhoff's theorem: 

Theorem 2.1. Let {M.,g) be a spherically symmetric solution to the vacuum 
equations ([7]) . Then it is locally isometric to a Schwarzschild solution with 
parameter M , for some M € R. 

In particular, spherically symmetric solutions to Q possess an additional 
Killing field not in the Lie algebra so(3). (Exercise: Prove Theorem 12.11 
Formulate and prove a global version of the result.) 

2.6.4 Higher dimensions 

In 3 + 1 dimensions, spherical symmetry is the only symmetry assumption com- 
patible with asymptotic flatness (see Appendix IB. 2. 3[) . such that moreover the 
symmetry group acts transitively on 2-dimensional orbits. Thus, Birkhoff's the- 
orem means that vacuum gravitational collapse cannot be studied in a 1 + 1 
dimensional setting by imposing symmetry. The simplest models for dynamic 
gravitational collapse thus necessarily involve matter, as in the Oppenheimer- 
Snyder model |125j or the Einstein-scalar field system studied by Christodoulou 

[nmg 

Moving, however, to 4 + 1 dimensions, asymptotically fiat manifolds can ad- 
mit a more general SU{2) symmetry acting transitively on 3-dimensional group 
orbits. The Einstein vacuum equations ((4|) under this symmetry admit 2 dy- 
namical degrees of freedom and can be written as a nonlinear system on a 
1 + 1-dimensional Lorentzian quotient Q = M./SU{2), where the dynamical de- 
grees of freedom of the metric are refiected by two nonlinear scalar fields on Q. 
This symmetry-known as "Triaxial Bianchi IX" -was first identified by Bizon, 
Chmaj and Schmidt [HI [T7| who derived the equations on Q and studied the 
resulting system numerically. The symmetry includes spherical symmetry as a 

^•^The discovery 1431 of these naked singularities led to the discovery of so-called critical 
collapse phenomena |37| which has since become a popular topic of investigation |87| . 
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special case, and thus, is admitted in particular by 4 + 1-dimensional Schwarz- 
schilcQ- The nonlinear stability of the Schwarzschild family as solutions of the 
vacuum equations Q can then be studied-within the class of Triaxial Bianchi 
IX initial data-as a 1 + 1 dimensional problem. Asymptotic stability for the 
Schwarzschild spacetime in this setting has been recently shown in the thesis of 
Holzegel j93[ [SH [M] , adapting vector field multiplier estimates similar to Sec- 
tion [4] to a situation where the metric is not known a priori. The construction of 
the relevant mutipliers is then quite subtle, as they must be normalised "from 
the future" in a bootstrap setting. The thesis [53] is a good reference for under- 
standing the relation of the linear theory to the non-linear black hole stability 
problem. See also Open problem [T3] in Section 

2.7 Geodesic incompleteness and "singularities" 

Is the picture of gravitational collapse as exhibited by Schwarzschild (or bet- 
ter, Oppcnhcimcr-Snydcr) stable? This question is behind the later chapters 
in the notes, where essentially the considerations hope to be part of a future 
understanding of the stability of the exterior region up to the event horizon, 
i.e. the closure of the past of null infinity to the future of a Cauchy surface. 
(See Section [STB] for a formulation of this open problem.) What is remarkable, 
however, is that there is a feature of Schwarzschild which can easily be shown 
to be "stable", without understanding the p.d.e. aspects of 1^: its geodesic 
incompleteness . 

2.7.1 Trapped surfaces 

First a definition: Let {M,g) be a time-oriented Lorentzian manifold, and S a 
closed spacelike 2-surface. For any point p & S, we may define two null mean 
curvatures trx and trx, corresponding to the two future-directed null vectors 
n(a;), n(x), where n, n are normal to S at x. Wc say that S is trapped if trx < 0, 
trx < 0. 

Exercise: Show that points p e Q \ clos(J"(I^)) correspond to trapped 
surfaces of Ai. Can there be other trapped surfaces? (Refer also for instance 

to [n].) 

2.7.2 Penrose's incompleteness theorem 

Theorem 2.2. (Penrose 1965 UMI) Let {M,g) he globally hyperboli^ with 
non- compact Cauchy surface E, where g is a metric, and let 

R^^V^'V > (7) 

for all null vectors V . Then if A4 contains a closed trapped two-surface S , it 
follows that {A4,g) is future causally geodesically incomplete. 

^''Exercise: Work out explicitly the higher dimensional analogue of the Schwarzschild 
solution for all dimensions. 
^^See Appendix [A] 
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This is the celebrated Penrose incompleteness theorem. 

Note that solutions of the Einstein vacuum equations (|4]) satisfy d?]) . (In- 
equality d?]), known as the null convergence condition, is also satisfied for so- 
lutions to the Einstein equations coupled to most plausible matter models, 
specifically, if the energy momentum tensor T^j, satisfies Tf^^V^V^ > for all 
null V^.) On the other hand, by definition, the unique solution to the Cauchy 
problem (the so-called maximal Cauchy development of initial data) is glob- 
ally hyperbolic (see Appendix IB.3[) . Thus, the theorem applies to the maximal 
development of (say) asymptotically flat (see Appendix IB.2.3[) vacuum initial 
data containing a trapped surface. Note finally that by Cauchy stability [OT] . 
the presence of a trapped surface in Ai is clearly "stable" to perturbation of 
initial data. 

From the point of view of gravitational collapse, it is more appropriate to 
define a slightly different notion of trapped. We restrict to iS c S a Cauchy 
surface such that S bounds a disc in S. We then can define a unique outward null 
vector field n along S, and we say that S is trapped if trx < and antitrappecQ 
if trx < 0, where trx denotes the mean curvature with respect to a conjugate 
"inward" null vector field. The analogue of Penrose's incompleteness theorem 
holds under this definition. One may also prove the interesting result that 
antitrapped surface cannot not form if they are not present initially. See [49] . 

Note finally that there are related incompleteness statements due to Penrose 
and Hawking [21] relevant in cosmological (see Section (5]) settings. 

2.7.3 "Singularities" and strong cosmic censorship 

Following [49] , we have called Theorem 12.21 an "incompleteness theorem" and 
not a "singularity theorem" . This is of course an issue of semantics, but let 
us further discuss this point briefly as it may serve to clarify various issues. 
The term "singularity" has had a tortuous history in the context of general 
relativity. As we have seen, its first appearance was to describe something that 
turned out not to be a singularity at all-the "Schwarzschild singularity" . It was 
later realised that behaviour which could indeed reasonably be described by the 
word "singularity" did in fact occur in solutions, as exemplified by the r = 
singular "boundary" of Schwarzschild towards which curvature scalars blow up. 
The prcscnse of this singular behaviour "coincides" in Schwarzschild with the 
fact that the spacetime is future causally geodesically incomplete-in fact, the 
curvature blows up along all incomplete causal geodesies. In view of the fact 
that it is the incompleteness property which can be inferred from Theorem 12. 2|, 
it was tempting to redefine "singularity" as geodesic incompleteness (see [HI]) 
and to call Theorem 12.21 a "singularity theorem" . 

This is of course a perfectly valid point of view. But is it correct then to 
associate the incompleteness of Theorem 12.21 to "singularity" in the sense of 
"breakdown" of the metric? Breakdown of the metric is most easily understood 
with curvature blowup as above, but more generally, it is captured by the no- 

^^Note that there exist other conventions in the hterature for this terminology. See | 12| . 
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tion of "inextendibility" of the Lorentzian manifold in some regularity class. We 
have already remarked that maximally-extended Schwarzschild is inextendible 
in the strongest of senses, i.e. as a C" Lorentzian metric. It turns out, how- 
ever, that the statement of Theorem I2.2[ even when applied to the maximal 
development of complete initial data for (UD, is compatible with the solution 
being extendible as a C°° Lorentzian metric such that every incomplete causal 
geodesic of the original spacetime enter the extension! This is in fact what hap- 
pens in the case of Kerr initial data. (See Sect ion [5TT] for a discussion of the Kerr 
metric.) The reason that the existence of such extensions does not contradict 
the "maximality" of the "maximal development" is that these extensions fail to 
be globally hyperbolic, while the "maximal development" is "maximal" in the 
class of globally hyperbolic spacetimes (see Theorem IB. 41 of Appendix [B|) . In 
the context of Kerr initial data. Theorem 12. 21 is thus not saying that breakdown 
of the metric occurs, merely that globally hyperbolicity breaks down, and thus 
further extensions cease to be predictable from initial dataFI 

A similar phenomenon is exhibited by the Reissner-Nordstrom solution of the 
Einstein-Maxwell equations |91| . which, unlike Kerr, is spherically symmetric 
and thus admits a Penrose diagram representation: 




What is drawn above is the maximal development of S. The spacetime is future 
causally geodcsically incomplete, but can be extended smoothly to a {M.,g) such 
that all inextendible geodesies leave the original spacetime. The boundary of 
(A^,g) in the extension corresponds to CT-C above. Such boundaries are known 
as Cauchy horizons. 

^^Further confusion can arise from the fact that "maximal extensions" of Kerr constructed 
with the help of analyticity are still geodesically incomplete and inextendible, in particular, 
with the curvature blowing up along all incomplete causal geodesies. Thus, one often talks of 
the "singularities" of Kerr, referring to the ideal singular boundaries one can attach to such 
extensions. One must remember, however, that these extensions arc of no relevance from the 
point of view of the Cauchy problem, and in any case, their singular behaviour in principle 
has nothing to do with Theorem 12.21 
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The strong cosmic censorship conjecture says that the maximal development 
of generic asymptotically flat initial data for the vacuum Einstein equations 
is inextendible as a suitably regular Lorentzian metric!^ One can view this 
conjecture as saying that whenever one has geodesic incompleteness, it is due 
to breakdown of the metric in the sense discussed above. (In view of the above 
comments, for this conjecture to be true, the behaviour of the Kerr metric 
described above would have to be unstable to pcrturbationQ) Thus, if by the 
term "singularity" one wants to suggest "breakdown of the metric", it is only 
a positive resolution of the strong cosmic censorship conjecture that would in 
particular (generically) make Theorem 12.21 into a true "singularity theorem" . 



2.8 Christodoulou's work on trapped surface formation in 
vacuum 

These notes would not be complete without a brief discussion of the recent 
breakthrough by Christodoulou [53j on the understanding of trapped surface 
formation for the vacuum. 

The story begins with Christodoulou's earlier }41j . where a condition is given 
ensuring that trapped surfaces form for spherically symmetric solutions of the 
Einstein-scalar field system. The condition is that the difference in so-called 
Hawking mass m of two concentric spheres on an outgoing null hupersurface be 
sufficiently large with respect to the difference in area radius r of the spheres. 
This is a surprising result as it shows that trapped surface formation can arise 
from initial conditions which are as close to dispersed as possible, in the sense 
that the supremum of the quantity 2m /r can be taken arbitrarily small initially. 

The results of [JT] lead immediately (see for instance HI] ) to the existence of 
smooth spherically symmetric solutions of the Einstein-scalar field system with 
Penrose diagram 

, r = 0„ ^ 




^*As with weak cosmic censorship, the original formulation of this conjecture is due to 
Penrose |128) . The formulation given here is from |47j . Related formulations are given in |54l 
I118| . One can also pose the conjecture for compact initial data, and for various Einstein- 
matter systems. It should be emphasized that "strong cosmic censorship" does not imply 
"weak cosmic censorship". For instance, one can imagine a spacetime with Penrose diagram 
as in the last diagram of Section 12.6.21 with incomplete X+ , but still inextendible across the 
null "boundary" emerging from the centre. 

^^Note that the instability concerns a region "far inside" the black hole interior. The black 
hole exterior is expected to be stable (as in the formulation of Section l5.6l l. hence these notes. 
See 1581 1591 for the resolution of a spherically symmetric version of this problem, where the 
role of the Kerr metric is played by Reissner-Nordstrom metrics. 
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where the point p depicted corresponds to a trapped surface, and the spacetime 
is past geodesically complete with a complete past null infinity, whose future is 
the entire spacetime, i.e., the spacetime contains no white holes F"l Thus, black 
hole formation can arise from spacetimes with a complete regular pastF^ 

In [S31 , Christodoulou constructs vacuum solutions by prescribing a charac- 
teristic initial value problem with data on (what will be) I". This I" is taken 
to be past complete, and in fact, the data is taken to be trivial to the past of 
a sphere on X^. Thus, the development will include a region where the met- 
ric is Minkowski, corresponding precisely to the lower lighter shaded triangle 
above. It is shown that-as long as the incoming energy per unit solid angle 
in all direction^ is sufficiently large in a strip of I" right after the trivial 
part, where sufficiently large is taken in comparison with the affine length of 
the generators of I~~a trapped surface arises in the domain of development of 
the data restricted to the past of this strip. Comparing with the spherically 
symmetric picture above, this trapped surface would arise precisely as before in 
the analogue of the darker shaded region depicted. 

In contrast to the spherically symmetric case, where given the lower triangle, 
existence of the solution in the darker shaded region (at least as far as trapped 
surface formation) follows immediately, for vacuum collapse, showing the exis- 
tence of a sufficiently "big" spacetime is a major difficulty. For this, the results 
of [53] exploit a hierarchy in the Einstein equations in the context of what 
is there called the "short pulse method". This method may have many other 
applications for nonlinear problems. 

One could in principle hope to extend |53j to show the formation of black 
hole spacetimes in the sense described previously. For this, one must first extend 
the initial data suitably, for instance so that I" is complete. If the resulting 
spacetime can be shown to possess a complete future null infinity I^, then, 
since the trapped surface shown to form can be proven (using the methods of 
the proof of Theorem l2.2[) not to be in the past of null infinity, the spacetime will 
indeed contain a black hole region!^ Of course, resolution of this problem would 
appear comparable in difficulty to the stability problem for the Kerr family (see 
the formulation of Section [5TB|) . 

^"The triangle "under" the darker shaded region can in fact be taken to be Minkowski. 

^^The singular boundary in general consists of a possibly empty null component emanating 
from the regular centre, and a spacelike component where r = in the limit and across which 
the spacetime is inextendible as a C" Lorentzian metric. (This boundary could "bite off" the 
top corner of the darker shaded rectangle.) The null component arising from the centre can be 
shown to be empty gcncrically after passing to a slightly loss regular class of solutions, for which 
well-poscdness still holds. See Christodoulou's proof of the cosmic censorship conjectures 1451 
for the Einstcin-scalar field system. 

■^■^This is defined in terms of the shear of X". 

■^''in spherical symmetry, the completeness of null infinity follows immediately once a single 
trapped surface has formed, for the Einstein equations coupled to a wide class of matter 
models. See for instance 1601 . For vacuum collapse, Christodoulou has formulated a statement 
on trapped surface formation that would imply weak cosmic censorship. See |47| . 



27 



3 The wave equation on Schwarzschild I: uni- 
form boundedness 



In the remainder of these lectures, we will concern ourselves solely with linear 
wave equations on black hole backgrounds, specifically, the scalar linear homo- 
geneous wave equation ([1]). As explained in the introduction, the study of the 
solutions to such equations is motivated by the stability problem for the black 
hole spacetimes themselves as solutions to (|4]) . The equation can be viewed 
as a poor man's linearisation of neglecting tensorial structure. Other linear 
problems with a much closer relationship to the study of the Einstein equations 
will be discussed in Section [H] 

3.1 Preliminaries 

Let {A4,g) denote (maximally-extended) Schwarzschild with parameter M > 0. 
Let S be an arbitrary Cauchy surface^ that is to say, a hypersurfacc with the 
property that every inextendible causal geodesic in A4 intersects S precisely 
once. (See Appendix [Al) 

Proposition 3.1.1. Ift]) e ij;' e then there is a unique ^ with 

tpls e H^^^{S), ngipls e H^^^{S), for all spacelike S c M, satisfying 

where denotes the future unit normal of S. For m > 1, i/ il) e H^^^ , e 
H{^^, then € ns^\s e HIS^{S). Moreover, if^\>iA'i, and ij^a,^^ 

are as above and xj^i = 1^)2, '^'i = '^'2 in an open set lA c E, then -01 = ■02 in 
M\J^iE\c\os{U)). 

We will be interested in understanding the behaviour of in the exterior of 
the black hole and white hole regions, up to and including the horizons. It is 
enough of course to understand the behaviour in the region 

V = c\os{J-{lX) n J^{Ia)) n Q 

where Ij^ denote a pair of connected components of I'"' , respectively, with a 
common limit point Ifj 

Moreover, it suffices (Exercise: Why?) to assume that S n = 0, and 
that we are interested in the behaviour in J"(I^) n J+(I]). Note that in this 
case, by the domain of dependence property of the above proposition, we have 
that the solution in this region is determined by ij'lpns, il'lcns- In the case 

^*We will sometimes be sloppy with distinguishing between n^^(p) and p, where tt : A4 Q 
denotes the natural projection, distinguishing J^(p) and J (p) n Q, etc. The context should 
make clear what is meant. 
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where E itself is spherically symmetric, then its projection to Q will look like: 

■r = 




r = 



If E is not itself spherically symmetric, then its projection to Q will in general 
have open interior. Nonetheless, wc shall always depict S as above. 

3.2 The Kay— Wald boundedness theorem 

The most basic problem is to obtain uniform boundedness for ip. This is resolved 
in the celebrated: 

Theorem 3.1. Let tp, ijj, t|)' be as in Proposition \3JJl with ^\> 6 iJj'^+i(E), 
ij)' € i7j™p(E) for a sufficiently high m, and such that t|), \\>' decay suitably at i^ . 
Then there is a constant D depending on xJj, t|)' such that 

in v. 

The proof of this theorem is due to Wald HM] and Kay- Wald [SHI- The 
"easy part" of the proof (Section I3.2.3P is a classic application of vector field 
commutators and multipliers, together with elliptic estimates and the Sobolcv 
inequality. The main difficulties arise at the horizon, and these arc overcome by 
what is essentially a clever trick. In this section, wc will go through the original 
argument, as it is a nice introduction to vector field multiplier and commutator 
techniques, as well as to the geometry of Schwarzschild. We will then point out 
(Section I3.2.7P various disadvantages of the method of proof. Afterwards, we 
give a new proof that in fact achieves a stronger result (Theorem 13. 2p . As we 
shall see, the techniques of this proof will be essential for future applications. 

3.2.1 The KilHng fields of Schwarzschild 

Recall the symmetries of {Ai,g): {M,g) is spherically symmetric, i.e. there 
is a basis of Killing vectors {ili}'^^^ spanning the Lie algebra so(3). These are 
sometimes known as angular momentum operators. In addition, there is another 
Killing field T (equal to dt in the coordinates dH) ) which is hypersurface orthog- 
onal and future directed timelike near iq. This Killing field is in fact timelike 
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everywhere in J~(I^) n J^(I~)^ becoming null and tangent to the horizon, van- 
ishing at T-L^ nT-l' . We say that the Schwarzschild metric in J^{I^) n J^(I^) is 
static. T is spacelike in the black hole and white hole regions. 

Note that whereas in Minkowski space R'^'*'^, the Killing fields at any point 
span the tangent space, this is no longer the case for Schwarzschild. Wc shall 
return to this point later. 

3.2.2 The current and its energy estimate 

Let (ft denote the 1-paramctcr group of diffcomorphisms generated by the Killing 
field T. Define S,- = (/3t(S n V). Wc have that {St-}t->o defines a spacelike 
foliation of 

Define 

and 

7?.(0,t) = Uo<t<tSi=. 

Let denote the future directed unit normal of E, and let n!^ define a null 
generator of l-L^ , and give 'W' the associated volume formF^ 

Let denote the energy current defined by applying the vector field T 

as a multiplier, i.e. 

with its associated current K'^{ip), 

where T^i^ denotes the standard energy momentum tensor of i/j (sec Appendix [D|i. 
Since T is Killing, and V^T^y = 0, it follows that K'^{ip) = 0, and the divergence 
theorem (See Appendix [C| applied to jj in the region TZ{0,t) yields 

L '"^^^""-^ ' L(o.) -^^'^^^"^ - L 

Recall that for null surfaces, the definition of a volume form relies on the choice of a 
normal. All integrals in what follow will always be with respect to the natural volume form, 
and in the case of a null hypersurface, with respect to the volume form related to the given 
choice of normal. See Appendix [C] 
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See 

r = 




r = 



Since T is future-directed causal in T), we have 

JJ(^X > 0, > 0. (9) 

Let us fix an tq > 2M. It follows from ([SD, (HD that 

-/STnlrSro} -'So 

As long as -g{T,n^^) < B for some constant we have 

i?(ro,S)((aiV)' + (a.V)' + ll^V'l')> jJ(V')ri^>K^o,S)((5tV)' + (5.V0' + ll^V'l')- 

Here, l^f'V'P denotes the induced norm on the group orbits of the 50(3) action, 
with 'p the gradient of the induced metric on the group orbits. We thus have 

[ {dti,f + {dri'f + \n?<B{r,,T.)f jl{^)n^^^. 

JT.^n{r>ro} -/So 

3.2.3 T as a commutator and pointwise estimates away from the 
horizon 

We may now commute the equation with T (See Appendix [E]) . i.e., since [□, T] = 
0, if Ugi) = then ng{Tip) = 0. We thus obtain an estimate 

[ (df^P)' + {drdt^f + \fdM^ < B{ro,^) [ JJ(TV'X. (10) 

JS-rn{r>ro} JSo " 

Exercise: By elliptic estimates and a Sobolev estimate show that if ^\>{x) 
as X ^ io, then (jlOp implies that for r > rg, 

< B{ro,^) (X^ Jj(^)<„ + X„ ^J(^^)<o) ' (11) 

for solutions oi nigip = 0. 

The right hand side of (|lip is finite under the assumptions of Theorem 13. 1[ 
for m = 1. Thus, proving the estimate of Theorem 13.11 awav from the horizon 

■^^For definiteness, one could choose S to bo a surface of constant t* defined in Section l2.2l 
or alternatively, require that it be of constant t for large r. 
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poses no difficulty. The difficulty of Theorem 13.11 is obtaining estimates which 
hold up to the horizon. 

Remark: The above argument via elliptic estimates clearly also holds for 
Minkowski space. But in that case, there is an alternative "easier" argument, 
namely, to commute with all translations We see thus already that the lack 
of Killing fields in Schwarzschild makes things more difficult. We shall again 
return to this point later. 

3.2.4 Degeneration at the horizon 

As one takes ro 2M, the constant i3(ro,S) provided by the estimate (|lip 
blows up. This is precisely because T becomes null on "H^ and thus its control 
over derivatives of "0 degenerates. Thus, one cannot prove uniform boundedness 
holding up to the horizon by the above. 

Let us examine more carefully this degeneration on various hypcrsurfaces. 

On E^, we have only 

^J(^X. ^ + (1 - 2M/r){dr^f + m^). (12) 

Wc see the degeneration in the presence of the factor (1 - 2AI/r). Note that 
(Exercise) 1 - 2M/r vanishes to first order on Ti^ \ TiT . Alternatively, one 
can examine the flux on the horizon T-i^ itself. For definiteness, let us choose 
nu* =T inUnn^. We have 

J^{i')T^ = im'. (13) 

Comparing with the analogous computation on a null cone in Minkowski space, 
one sees that a term \'fip\'^ is "missing". 

Are estimates of the terms (|12p . (|13p enough to control ipl It is a good 
idea to play with these estimates on your own, allowing yourself to commute 
the equation with T and Hi to obtain higher order estimates. Exercise: Why 
docs this not lead to an estimate as in (fTT|) ? 

It turns out that there is a way around this problem and the degeneration 
on the horizon is suggestive. For suppose there existed a ijj such that 

□ gV; = 0, T4,= (14) 

Let us see immediately how one can obtain estimates on the horizon itself. For 
this, we note that 

Commuting now with the whole Lie algebra of isometrics . we obtain 

i 



Easier, but not necessarily better. . . 
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Clearly, by a Sobolev estimate applied on the horizon, together with the estimate 

for T = T,fli (here (a) denotes a multi-index of arbitrary order), we would 
obtain 

\^\'^B EEL ^J(r("^^)<o (15) 

r=T,n, |(q)|<2 -^^0 

on n n. 

It turns out that the estimate (|15p can be extended to points not on the 
horizon by considering t = c surfaces. Note that these hypersurfaces all meet 
at 'H'^ n "W". It is an informative calculation to examine the nature of the de- 
generation of estimates on such hypersurfaces because it is of a double nature, 
since, in addition to T becoming null, the limit of (subsets of) these space- 
like hypersurfaces approaches the null horizon Ti^. We leave the details as an 
exercise. 

3.2.5 Inverting an elliptic operator 

So can a satisfying (|14p actually be constructed? We have 

Proposition 3.2.1. Suppose m is sufficiently high, iji, decay suitably at , 
and '^p\■^+n^H- = 0? ■^■V'lw+n'H- = for some spherically symmetric timelike vector 
field S defined along 'W' n TiT . Then there exists a ip satisfying dgip = with 
Tip = ip in T>, and moreover, the right hand side of {\15^ is finite. 

Formally, one sees that on t = c say, if we let g denote the induced Riemannian 
metric, and if we impose initial data 

Tip\t=c = ip, 

V'|t=c = A[l-2M/r)-ig^V', 

and let ip solve the wave equation with this data, then 

T-tp = -0 

as desired. 

So to use the above, it suffices to ask whether the initial data for -ip above 
can be constructed and have sufficient regularity so as for the right hand side of 
(fT5|) to be defined. To impose the first condition, since T = along H'^n'H^, one 
must have that ip vanish there to some order. For the second condition, note 
first that the metric (1 - 2M/r)^^g has an asymptotically hyperbolic end and 
an asymptotically flat end. Thus, to construct ^jl_2M/r)-^g'^^ suitably well- 
behavec(ff|, one must have that Tip decays appropriately towards the ends. We 
leave to the reader the task of verifying that the assumptions of the Proposition 
are sufficient. 

•^^so that we may apply to this quantity the arguments of Section 13.2.41 
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3.2.6 The discrete isometry 



Proposition 13.2.11 together with estimates (|15|) and (jlip , yield the proof of 
Theorem 13. II in the special case that the conditions of Proposition [21211] happen 
to be satisfied. In the original paper of Wald [151] . one took Sq to coincide with 
t = and restricted to data ij), which were supported in a compact region not 
containing nTL^. Clearly, however, this is a deficiency, as general solutions 
will be supported in n H^. (See also the last exercise below.) 

It turns out, however, that one can overcome the restriction on the support 
by the following trick: Note that the previous proposition produces a ijj such 
that Tip = ip on all of V. We only require however that Tip = ip onTZ. The idea 
is to define a new tJ), on E, such that i|) = tJ;, i];' = on Eg and, denoting 
by ■0 the solution to the Cauchy problem with the new data, ip\f{+nn- = 0, 
S'i/'lw+n-H- = 0. By the previous proposition and the domain of dependence 
property of Proposition 13.1. 11 we will have indeed constructed a ip with Tip = ip 
in TZ for which the right hand side of (jlSp is finite. 

Remark that Schwarzschild admits a discrete symmetry generated by the 
map X -X in Kruskal coordinates. Define ip, ij;' so that t|)(X, ■) = •), 
^'{X,.) = -^'{-X,.). 

Proposition 3.2.2. Under the above assumptions, it follows that 

iP{X,-) = -iP{-X,-). 

The proof of the above is left as an exercise in preservation of symmetry for 
solutions of the wave equation. It follows immediately that 

and that 

duip = -dyip, 

and thus {djj + dv)ip = 0. In view of the above remarks and Proposition 13.2.11 
with 'E. = djj + dv, we have shown the full statement of Theorem 13. II 

Exercise: Work out explicit regularity assumptions and quantitative depen- 
dence on initial data in Theorem l3.1l describing in particular decay assumptions 
necessary at ig. 

Exercise: Prove the analogue of Theorem 13.11 on the Oppcnheimcr-Snydcr 
spacetime discussed previously Hint: One need not know the explicit form of 
the metric, the statement given about the Penrose diagram suffices. Convince 
yourself that the original restricted version of Theorem 13.11 due to Wald |151j , 
where the support of ip is restricted near TL'^ n TL^ , is not sufficient to yield this 
result. 

3.2.7 Remarks 

The clever proof described above successfully obtains pointwise boundedness 
for Ip up to the horizon "H"^. Does this really close the book, however, on the 
boundedness question? From various points of view, it may be desirable to go 
further. 
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1. Even though one obtains the "correct" pointwise result, one does not 
obtain boundedness at the horizon for the energy measured by a local 
observer, that is to say, bounds for 

This indicates that it would be difficult to use this result even for the 
simplest non-linear problems. 

2. One does not obtain boundedness for transverse derivatives to the horizon, 
i.e. in {t*,r) coordinates, drip, d^^J, etc. (Exercise: Why not?) 

3. The dependence on initial data is somewhat unnatural. (Exercise: Work 
out explicitly what it is.) 

As far as the method of proof is concerned, there are additional short- 
comings when the proof is viewed from the standpoint of possible future 
generalisations: 

4. To obtain control at the horizon, one must commute (sec (|15p ) with all an- 
gular momentum operators Qi. Thus the spherical symmetry of Schwarz- 
schild is used in a fundamental way. 

5. The exact staticity is fundamental for the construction of ip. It is not clear 
how to generalise this argument in the case say where T is not hypersurface 
orthogonal and Killing but one assumes merely that its deformation tensor 
"^77^^ decays. This would be the situation in a bootstrap setting of a non- 
linear stability problem. 

6. The construction of ip requires the discrete isometry of Schwarzschild, 
which again, cannot be expected to be stable. 

3.3 The red-shift and a new proof of boundedness 

We give in this section a new proof of boundedness which overcomes the short- 
comings outlined above. In essence, the previous proof limited itself by relying 
solely on Killing fields as multipliers and commutators. It turns out that there is 
an important physical aspect of Schwarzschild which can be captured by other 
vector-field multipliers and commutators which are not however Killing. This 
is related to the celebrated red-shift effect. 

3.3.1 The classical red-shift 

The red-shift effect is one of the most celebrated aspects of black holes. It is 
classically described as follows: Suppose two observers, A and B are such that 
A crosses the event horizon and B does not. If A emits a signal at constant 
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frequency as he measures it, then the frequency at which it is received by B is 
"shifted to the red" . 



The consequences of this for the appearance of a coUapsing star to far-away 
observers were first explored in the seminal paper of Oppenheimer- Snyder [125j 
referred to at length in Section [2l For a nice discussion, see also the classic 
textbook [TT7] . 

The red-shift effect as described above is a global one, and essentially de- 
pends only on the fact that the proper time of B is infinite whereas the proper 
time of A before crossing is finite. In the case of the Schwarzschild black 
hole, there is a "local" version of this red-shift: If B also crosses the event 
horizon but at advanced time later than A: 



then the frequency at which B receives at his horizon crossing time is shifted to 
the red by a factor depending exponentially on the advanced time difference of 
the crossing points of A and B. 

The exponential factor is determined by the so-called surface gravity^ a quan- 
tity that can in fact be defined for all so-called Killing horizons. This localised 
red-shift effect depends only on the positivity of this quantity. We shall under- 
stand this more general situation in Section [T] Let us for now simply explore 
how we can "capture" the red-shift effect in the Schwarzschild geometry. 

3.3.2 The vector fields N, Y, and Y 

It turns out that a "vector field multiplier" version of this localised red-shift 
effect is captured by the following 

Proposition 3.3.1. There exists a ipt-invariant smooth future- directed timelike 
vector field N on TZ and a positive constant b> such that 





K^{^)>bj;^{m'' 



on T-L^ for all solutions "0 of Hgip = 0. 
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(See Appendix ID] for the , notation.) 

Proof. Note first that since T is tangent to T-L^ , it follows that given any cr < oo, 
there clearly exists a vector field Y ox\TZ such that 

1. Y IS Lpt invariant and spherically symmetric. 

2. Y is future-directed null on T-L^ and transverse to T-L^ , say g{T,Y) = -2. 

3. OnU^, 

VyY =-(j{Y + T). (16) 
Since T is tangent to , along which Y is null, we have 

g(VTr,r) = o. (17) 

From properties[T]and[21 and the form of the Schwarzschild metric, one computes 
(Exercise) 

g{VTY,T)=2K>Q (18) 
on T-L^ . Defining a local frame E2 for the S0(3) orbits, we note 

g{VEXY) = ]^E,g{Y,Y) = Q, 

g{VE^Y,E2) = -g{Y,VE,E2) = -g{Y,VE,Ei) = g{VE,Y,Ei). 
Writing thus 

VtY = -kY + a^Ei+ E2 (19) 
VyY = -aT-aY (20) 

Ve,Y =h\Ei + hlE2-]^a^Y (21) 

Ve^Y =h\Ei+hlE2-hi^Y (22) 
with {h\ = h^), we now compute 

= -i(T(y,y)(-«) + iT(T,y)(-a) + iT(T,r)(-a)) 

-^{T{Ei,Y)a' + T{E2,Y)a^) 

+ T{EuEi)hl+T{E2,E2)hl + T{EuE2){hl + hl) 

where we denote the energy momentum tensor by T, to prevent confusion with 
T. (Note that, in view of the fact that Q imbeds as a totally geodesic sub- 
manifold oi M, we have in fact = = 0. This is of no importance in our 
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computations, however.) It follows immediately in view again of the algebraic 
properties of T, that 

- cT{T, Y + T)- cv'T(r,r + r)T(r,y) 

where c is independent of the choice of cr. It follows that choosing a large 
enough, we have 

So just set TV = T + r, noting that = + = . □ 

The computation (|18p represents a well known property of stationary black 
holes holes and the constant k is the so-called surface gravity. (See |148j .) Note 
that since Y is (/9f-invariant and T is Killing, we have 

giVrY, T) = g( VyT, T) = -giVrT, Y) 

on H^. On the other hand 

g{yTT,E,) = -g{VE,T,T) = 0, 

since T is null on "H^. Thus, k is alternatively characterized by 

VtT=kT 

on "H^ . We will elaborate on this in Section [71 where a generalisation of Propo- 
sition 13.3.11 will be presented. 

Exercise: Relate the strength of the red-shift with the constant k, for the 
case where observers A and B both cross the horizon, but B at advanced time 
V later than A. 

If one desires an explicit form of the vector field, then one can argue as 
follows: Define first the vector field Y by 

y = ^—r9u- (23) 

(See Appendix [FT) Note that this vector field satisfies g(Vyi^,T) = 0. Define 
Y = {l + 5i{r-2M))Y + 52{r-2M)T. 

It suffices to choose (5i, 82 appropriately. 

The behaviour of N away from the horizon is of course irrelevant in the 
above proposition. It will be useful for us to have the following: 

Corollary 3.1. Let S he as before. There exists a ipt-invariant smooth future- 
directed timelike vector field N on TZ, constants b > 0, B > 0, and two values 
2M < ro < ri < 00 such that 

1. yhJ^ni^ forr<ro, 

2. N = T for r>ri, 

3. \K^\ < BJ^n^, and jf^ni^ ~ j'^n'^ for r^<r<ri. 
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3.3.3 iV as a multiplier 

Recall the definition of TZ{0,t). Applying the energy identity with the current 
in this region, we obtain 

Js^ ' ^ -VW+(0,r) ^ " J{r<ro}mZ{0,T) 

= f {-K'')+ f Jun^. (24) 

The reason for writing the above identity in this form will become apparent in 
what follows. Note that since N is timelike at 'H*, we see all the "usual terms" 
in the flux integrals, i.e. 

J>^~(9*.^)2 + |1(^VP, 

and 

j>^^~(9,.^)2 + (a.^)2 + |^^|2. 

The constants in the ~ depend as usual on the choice of the original Eq and the 
precise choice of N. 

Now the identity (|24p also holds where Sq is replaced by S^/, 7{^(0,r) is 
replaced by H'^{t\t), and TZ{0,t) is replaced by TZ{t',t), for an arbitrary 
< r' < T. 

We may add to both sides of (|24p an arbitrary multiple of the spacctime 
integral f{r>r„}nn{T'.T) "^J"-!;- ^^'^'^ ^hc fact that 

f «~ [If 

J {r>r'}nn{T' ,t) ' Jt' \J{r>r'}n^^ ^ 

for any r' > 2M (where ~ depends on Eqj -^)j from the inequalities shown and 
property [3] of Corollary 13. II we obtain 

On the other hand, in view of our previous dS]), ©, we have 

X:(X..'J"£)<'^^(--')X/Jn£- (25) 



Setting 



we have that 



fir) + b £ f{f)df < BD{t -t') + fir') (26) 

for aU r > r' > 0, from which it follows (Exercise) that / < B{D + /(O)). (We 
use the inequality with D = f.^^ '^J'^So') view of the trivial inequahty 



X,/J»-^^X„ 



Jfj. "-So' 
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we obtain 



X J>k<-Bj^ J>k- (27) 



Wc have obtained a "local observer's" energy estimate. This addresses 
point [U of Section 

3.3.4 y as a commutator 

It turns out (Exercise) that from (p7|) . one could obtain pointwise bounds as 
before on ij) by commuting with angular momentum operators No construc- 
tion of '0, '0, etc., would be necessary, and this would thus address points [31 [SJ 
iniof Section [SX71 

Commuting with f2i clearly would not address however point [J] Moreover, 
it would not address point [H Exercise: Why not? 

It turns out that one can resolve this problem by applying N not only as a 
multiplier, but also as a commutator. The calculations are slightly easier if we 
more simply commute with Y defined in (|23p . 

Proposition 3.3.2. Let ip satisfy Hgip = 0. Then we may write 

HgiYij) = -Y{Y{^)) - -{Y{T,p)) + Pi0 (28) 
r r 

where Pi is the first order operator Piip = -^{Tip - Ytp). 

This is proven easily with the help of Appendix [E] As we shall see, the sign 
of the first term on the right hand side of (|28p is important. We will interpret 
this computation geometrically in terms of the sign of the surface gravity in 
Theorem Oof Section [71 

Let us first note that our boundedness result gives us in particular 

f K^{iI;)<BDt (29) 

J {r<ro}nH{0,r) 

where D comes from initial data. (Exercise: Why?) Commute now the wave 
equation with T and apply the multiplier N. See Appendix [El One obtains in 
particular an estimate for 

f {YTi:y<Bf K^{Ti')<BDT, (30) 

J{r<ro}n7?.(0,r) J{r<ro}nK(0.r) 

where again D refers to a quantity coming from initial data. Commuting now 
the wave equation with Y and applying the multiplier N, one obtains an energy 
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identity of the form 

= [ (-i^^(y(V0) 

_/{ro<r<ri}n7?.(0,r) 

+ [ [ £''{Y^^) 

J{r<ro}n7?.(0,r) J{r>7-o}n7l(0,T) 

where J^(FV'), ^ST^C^V') are defined by dm]) . (fT^ . respectively, with 
replacing -0, and 



{Y^) = -{NY^)(-Y{Y{^))--{Y{Ti^)) + Pi'il;] 
r 

- \{N - Y)Yi>){YYi}) + -(7Vyi/;)(r(rV')) 
r r 

- (iVy7/')Pii/'. 

The first term on the right hand side has a good sign! Applying Cauchy-Schwarz 
and the fact that - F = T on 'H'^, it follows that choosing 7'o accordingly, one 
obtains that the second two terms can be bounded in r < ro by 

tK^{Y^) + t-^(YTi,f 
whereas the last term can be bounded by 

tK^(Yi)) + t-^K^{i,). 
In view of (l29l) and (1301). one obtains 



r e^{Yi,)<e [ K^{Yil:) + Be-^DT. 

J{r<ro}nK(0,r) J{r<ro}nTJ.(0,T) 

Exercise: Show how from this one can arrive again at an inequality (|26p . 

Commuting repeatedly with T, F, the above scheme plus elliptic estimates 
yield natural 7J™ estimates for all m. Pointwise estimates for all derivatives 
then follow by a standard Sobolev estimate. 



3.3.5 The statement of the boundedness theorem 

We obtain finally 

Theorem 3.2. Let T, be a Cauchy hypersurface for Schwarzs child such that 
E n T-i^ = 0, let TiQ = T) n T, , let T,r denote the translation of T,q, let n-^^ denote 
the future normal of T,r, and let TZ = Ut->oS7-. Assume -g{nso,T) is uniformly 
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bounded. Then there exists a constant C depending only on Sq such that the 
following holds. Let ip, '[\>, t|)' be as in Proposition ISJjl with ^\) e 
^\>'eH|:^^(E), and 

for < m < k. Then 

If k> 1 , then we have 

E E |(V^)™^r."^ VI < C( lim 1^1 + |V^(°)il;U.(So) + WIhh^o)] 

0<Tn<fc-l rni+Tn2=m,?rii>0 / 
in TZ. 

Note that {y^)"'^n"^^'ip denotes an mi-tcnsor on the Ricmannian manifold 
Et-, and I • I on the left hand side of the last inequality above just denotes the 
induced norm on such tensors. 

3.4 Comments and further reading 

The first discussion of the wave equation on Schwarzschild is perhaps the work 
of Regge and Wheeler [131| . but the true mathematical study of this problem 
was initiated by Wald |151| . who proved Theorem l3.1l under the assumption that 
ip vanished in a neighborhood of n TL^ . The full statement of Theorem 13. II 
and the proof presented in Section [3T2l is due to Kay and Wald [98]. The present 
notes owe a lot to the geometric view point emphasized in the works |151[ [98] . 

Use of the vector field F as a multiplier was first introduced in our [55] . 
and its use is central in [66j and [67j . In particular, the property formalised 
by Proposition 13 . 3 . Il was discovered there. It appears that this may be key to 
a stable understanding of black hole event horizons. See Section below, as 
well as Section [T] for a generalisation of Proposition 13.3.11 

It is interesting to note that in [66U67j . Y had always been used in conjunc- 
tion with vector fields X of the type to be discussed in the next section (which 
require a more delicate global construction) as well as T. This meant that one 
always had to obtain more than boundedness (i.e. decay!) in order to obtain 
the proper boundedness result at the horizon. Consequently, one had to use 
many aspects of the structure of Schwarzschild, particularly, the trapping to 
be discussed in later lectures. The argument given above, where boundedness 
is obtained using only N and T as multipliers is presented for the first time 
in a self-contained fashion in these lectures. The argument can be read off, 
however, from the more general argument of |68] concerning perturbations of 
Schwarzschild including Kerr. The use of Y as a commutator to estimate higher 
order quantities also originates in [68j . The geometry behind this computation 
is further discussed in Section [7] 
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Note that the use of Y together with T is of course equivalent to the use of 

and T. We have chosen to give a name to the vector field N = T + Y merely 
for convenience. Timelike vector fields are more convenient when perturbing. . . 

Another remark on the use of K as a commutator: Enlarging the choice of 
commutators has proven very important in previous work on the global analysis 
of the wave equation. In a seminal paper, Klainerman [100] showed improved 
decay for the wave equation on Minkowski space in the interior region by com- 
mutation with scaling and Lorentz boosts. This was a key step for further 
developments for long time existence for quasilinear wave equations [lOlj . 

The distinct role of multipliers and commutators and the geometric consid- 
erations which enter into their construction is beautifully elaborated by Christo- 
doulou [52]. 

3.5 Perturbing? 

Can the proof of Thcorcm l3.2l bc adapted to hold for spacetimes "near" Schwarz- 
schild? To answer this, one must first decide what one means by the notion of 
"near" . Perhaps the simplest class of perturbed metrics would be those that 
retain the same differentiable structure of TZ, retain 71^ as a null hypersurface, 
and retain the Killing field T. One infers (without computation!) that the state- 
ment of Proposition [3?3TT] and thus Corollary [SH] is stable to such perturbations 
of the metric. Therein lies the power of that Proposition and of the multiplier 
TV. (In fact, see Section [T)) Unfortunately, one easily sees that our argument 
proving Theorem 13.21 is still unstable, even in the class of perturbations just 
described. The reason is the following: Our argument relies essentially on an 
a priori estimate for Jjin^ (sec (|25p '). which requires T to be non-spacelike 
in IZ. When one perturbs, T will in general become spacelike in a region of TZ. 
(As we shall see in Section 15.11 this happens in particular in the case of Kerr. 
The region where T is spacelikc is known as the ergoregion.) 

There is a sense in which the above is the only obstruction to perturbing the 
above argument, i.e. one can solve the following 

Exercise: Fix the differentiable structure of TZ and the vector field T. Let g 
be a metric sufficiently close to Schwarzschild such that Ti* is null, and suppose 
T is Killing and non-spacclike in TZ, and T is null on H^. Then Theorem 13.21 
applies. (In fact, one need not assume that T is non-spacelike in TZ, only that 
T is null on the horizon.) See also Section [Tj 

Exercise: Now do the above where T is not assumed to be Killing, but ^Tr^y 
is assumed to decay suitably. What precise assumptions must one impose? 

This discussion may suggest that there is in fact no stable boundedness 
argument, that is to say, a "stable argument" would of necessity need to prove 
more than boundedness, i.e. decay. We shall see later that there is a sense in 
which this is true and a sense in which it is not! But before exploring this, let us 
understand how one can go beyond boundedness and prove quantitative decay 
for waves on Schwarzschild itself. It is quantitative decay after all that we must 
understand if we are to understand nonlinear problems. 
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4 The wave equation on Schwarzschild II: quan- 
titative decay rates 

Quantitative decay rates are central for our understanding of non-linear prob- 
lems. To discuss energy decay for solutions ijj of □gV' = on Schwarzschild, one 
must consider a different foliation. Let Eg be a spacelike hypersurface termi- 
nating on null infinity and define St- by translation. 




The main result of this section is the following 

Theorem 4.1. There exists a constant C depending only on Eg such that the 
following holds. Let i\> e H^^^, € H^^^, and suppose lim^-^io tj) = and 

|(q)|<3 r={n,} 

where uq denotes the unit normal of the hypersurface {t = 0}. Then 

j;^Wn^l<CE,r-\ (31) 

where N is the vector field of Section \3.3.2\ Now let -\\> € -ffj^^., ^\)' € H^^^, 
lim^^io t]) = 0, and suppose 

E2= Z E f rV;«(r(")V')<<oo. 

|(Q)|<6r={f2.} -''-o 

Then 

sup < CV^T"\ supr\iP\<C^E^T-'^'^. (32) 

The fact that (|3ip "loses derivatives" is a fundamental aspect of this problem 
related to the trapping phenomenon, to be discussed in what follows, although 
the precise number of derivatives lost above is wasteful. Indeed, there are several 
aspects in which the above results can be improved. See Proposition 14.2.11 and 
the exercise of Section 14.31 
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We can also express the pointwise decay in terms of advanced and retarded 
null coordinates u and v. Defining! v = 2{t + r*) =2{t + T + 2Mlog(r - 2M)), 
u = 2{t-r*) = 2{t-r-2Mlog{r-2M)), it follows in particular from 1^ that 

\^P\<CE2{\v\ + l)-\ H\<C{ro)E2U-K (33) 

where the first inequality applies in I?nclos({i > 0}), whereas the second applies 
only in Vn{t > 0}n{r > rg}, with C(ro) ^- oo as tq ^- 2A1. See also Appendix iFl 
Note that, as in Minkowski space, the first inequality of (|33p is sharp as a 
uniform decay rate in v. 

4.1 A spacetime integral estimate 

The zero'th step in the proof of Theorem 14.11 is an estimate for a spacetime 
integral whose integrand should control the quantity 

X<(V)"|^ (34) 

where x is a (/s^-invariant weight function such that x degenerates only at infinity. 
Estimates of the spacetime integral (|34p have their origin in the classical virial 
theorem, which in Minkowski space essentially arises from applying the energy 
identity to the current with V = 

Naively, one might expect to be able to obtain an estimate of the form say 

f xJ^{i^)nl <B f J^nl , (35) 
Jh(o,t) ^ Jso ^0 

for such a x- It turns out that there is a well known high-frequency obstruc- 
tion for the existence of an estimate of the form (|35p arising from trapped null 
geodesies. This problem has been long studied in the context of the wave equa- 
tion in Minkowski space outside of an obstacle, where the analogue of trapped 
null geodesies are straight lines which reflect off the obstacle's boundary in such 
a way so as to remain in a compact subset of space. In Schwarzschild, one can 
easily infer from a continuity argument the existence of a family of null geodesies 
with i"'' as a limit point|3 But in view of the integrability of geodesic flow, one 
can in fact understand all such geodesies explicitly. 

Exercise: Show that the hypersurface r = 3M is spanned by null geodesies. 
Show that from every point in TZ, there is a codimcnsion-one subset of future 
directed null directions whose corresponding geodesies approach r = 3A/, and 
all other null geodesies either cross Ti^ or meet I^. 

The timelike hypersurface r = 3M is traditionally called the photon sphere. 
Let us first see how one can capture this high frequency obstruction. 

■^^The strange convention on the factor of 2 is chosen simply to agree with 1651 . 
^''This can be thought of as a very weak notion of what it would mean for a null geodesic 
to be trapped from the point of view of decay results with respect to the foliation Sr . 
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4.1.1 A multiplier X for high angular frequencies 

We look for a multiplier with the property that the spacetime integral it gener- 
ates is positive definite. Since in Minkowski space, this is provided by the vector 
field we will look for simple generalizations. Calculations are slightly easier 
when one considers dr* associated to Regge- Wheeler coordinates {r*,t). See 
Appendix IF. 21 for the definition of this coordinate svstemF^ For X = f{r*)dr*, 
where / is a general function, we obtain the formula 

Here /' denotes We can now define a "modified" current 



associated to the vector field X and the function w. Let 



Choosing 



we have 



Recall that in view of the spherical symmetry of A4, we may decompose 

£>0,\m\<i 

where Ym,i are the so-called spherical harmonics, each summand satisfies again 
the wave equation, and the convergence is in of the S0(3) orbits. 

Let us assume that '0^ = for spherical harmonic number £ < L for some L to 
be determined. We look for K^ "" such that /g^ K^''" > 0, but also | J^'"'n^| < 



Remember, when considering coordinate vector fields, one has to specify the entire coordi- 
nate system. When considering dr, it is here to be understood that we are using Schwarzschild 
coordinates, and when considering d^* , it is to be understood that we are using Rcgge-Wheeler 
coordinates. The precise choice of the angular coordinates is of course irrelevant. 
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BJ^n^. Here /g2 denotes integration over group orbits of the S0(3) action. 
For such tp, in view of the resulting inequality 

it follows that taking L sufficiently large and < 6 <1 sufficiently small so that 
1 ~ ^'^2r3^ - 5' it clearly suffices to construct an / with the following properties: 

1- \f\<B, 

2. /' > 5(1 - 2M/r)r-^, 

3. f{r = 3M) = 0, 

4. -in, (2/' + 41^/ + 2'-i^ (1 - 3M) /) (, = ui) > 0, 

5. (2/' + 41^/ + 2^il0^ (1 - ^) /) < Br-^ 

for some constant B. Exercise. Show that one can construct such a function. 
Note the significance of the photon sphere! 

4.1.2 A multiplier X for all frequencies 

Constructing a multiplier for all spherical harmonics, so as to capture in addition 
"low frequency" effects, is more tricky. It turns out, however, that one can ac- 
tually define a current which does not require spherical harmonic decomposition 
at all. The current is of the form: 

i 

1 rify (r-2M (r*-a-a''') 



2 pir - 2M) 



lr-2M {r*-a-a^l'^) \ 
\ + {r* - a - a^/'^)^ J 



Here, N is as in Section[S321 X'' = pdr^^, = fdr*, the warped current J^-"" 
is defined as in Section 14.1.11 

pa _ £a 

f=U ll^^Z^ _ tan-i(-l - a'l^-)] , 
a \ a I 

and e, Ca, Ca, a are positive parameters which must be chosen accordingly. 
With these choices, one can show (after some computation) that the divergence 
= S/^J^ controls in particular 

l^K'^Wybxl^Jl^Wn^, (36) 
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where x is non- vanishing but decays (polynomially) as r ^ oo . Note that in view 
of the normahsation (|123p of the r* coordinate, = precisely at r = 3M. The 
left hand side of the incquahty (|36p controls also second order derivatives which 
degenerate however at r = 3M. We have dropped these terms. It is actually 
useful for applications that the part of the current is not "modified" by 

a function w;", and thus ■0 itself does not occur in the boundary terms. That is 
to say 




On the event horizon 'H^, we have a better one-sided bound 




For details of the construction, see [67] . 

In view of (j36|) , (l37|) and (l38|) , together with the previous boundcdncss result 
Theorem 13.21 one obtains in particular the estimate 

^J^i^>t ^ B f^^^^^ ( + E J^cav-) j , (39) 

for some nonvanishing (ySfinvariant function x which decays polynomially as 
r -* oo. 

On the other hand, considering the current jJ^{P<Lip) + J^'^{{I - P<l)^), 
where J^'™ is the current of Section [4. 1.11 and P<l''P denotes the projection to 
the space spanned by spherical harmonics with £ < L, we obtain the estimate 

f xhJ^Wn'^^ <B f ^ J^{^)nl , (40) 

where h is any smooth nonnegativc function < h < 1 vanishing at r = 3M, and 
B depends also on the choice of function h. 

4.2 The Morawetz conformal Z multiplier and energy de- 
cay 

How does the estimate (|39p assist us to prove decay? 

Recall that energy decay can be proven in Minkowski space with the help of 
the so-called Morawetz current. Let 

Z = u^du + v^d^ (41) 

and define 

t7»,/,x .7/ ,x tr*(l - 2M/r) , ^ , r*(l-2M/r) 
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(Here {u,v), {r*,t) are the coordinate systems of Appendix |F]) Setting M = 
0, this corresponds precisely to the current introduced by Morawetz [119] on 
Minkowski space. 

It is a good exercise to show that (for AI > 0!) the coefficients of this 
current are C° but not C"^ across u H^. 

To understand how one hopes to use this current, let us recall the situation 
in Minkowski space. There, the signifance of (|4ip arises since it is a conformal 
Killing field. Setting M = 0, r* = r in the above one obtain^ 

f Jfrn^>^, (42) 

Jt=T 

K^"" = 0. (43) 

The inequality (|42p remains true in the Schwarzschild case and one can 
obtain exactly as before 

^ ^ J^'^n'^ >b£^ u^duiJ^ + v^id^i^y + (l - ^) (u^ + v^)\nf- (44) 

(In fact, we have dropped positive O'th order terms from the right hand side of 
(I44p . which will be useful for us later on in Section [4731 ) Note that away from 
the horizon, we have that 

f J^'-n'^>biro,Ry f ^ ^ ,J^n^. (45) 

Jt=T ^ J {t=T}n{ro<r<R} ^ 

Thus, if the left hand side of ([45[) could be shown to be bounded, this would 
prove the first statement of Theorem 14.11 where Sr is replaced however with 
{t = T}n{ro<r< R}. 

In the case of Minkowski space, the boundedness of the left hand side of ([44[) 
follows immediately by (I43p and the energy identity 

f Jf^"'+ [ K^''"= [ Jf'"' (46) 
Jt=T ^ Ja<t<T Jt=o 

as long as the data arc suitably regular and decay so as for the right hand side 
to be bounded. For Schwarzschild, one cannot expect ([43p to hold, and this is 
why we have introduced the X-related currents. 

First the good news: There exist constants vq < R such that 

for r < rg , and in fact 



K^'"" > h^iP^ (47) 



*^The reason for introducing the O'th order terms is because the wave equation is not confor- 
mally invariant. It is remarkable that one can nonetheless obtain positive definite boundary 
terms, although a slightly unsettling feature is that this positivity property (I42II requires 
looking specifically at constant t = t surfaces and integrating. 
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for r > i? and some constant b. These terms have the "right sign" in the energy 
identity (|46p . In {vq <r<R}, however, the best we can do is 



-K^^'" <Bt{\f,p\' + \i,\'). 

This is the bad news, although, in view of the presence of trapping, it is to be 
expected. Using also (|47p . we may estimate 



Jo<t<T J {0<t<T}n{ro<r<R} 

< Bt [ Jl^n''. (48) 

J {0<t<T}n{ro<r<R} 

In view of the fact that the first integral on the right hand side of (|48|) is 
bounded by (|39|) . and the weight in (|45|) . applying the energy identity of 
the current J'^'™ in the region < t < r, we obtain immediately a preliminary 
version of the first statement of the Theorem 14. 1[ but with replaced by r, and 
the hypcrsurfaces Sr replaced by {t = t} n {r' < r < R'} for some constants r', 
R', but where B depends on these constants. (Note the geometry of this region. 
All {t = constant} hypersurfaces have common boundary H"^ n . Exercise: 
Justify the integration by parts (|46p , in view of the fact that Z and w are only 
C° at WuH^.) 

Using the current and an easy geometric argument, it is not difficult to 
replace the hypersurfaces {t = t} n {r' < r < R'} above with Et- n {r > r'}@ 
obtaining 

Js^n{r>r'} \Jt=0 JSq ^ i^l ) 

(49) 

To obtain decay for the nondcgenerate energy near the horizon, note that 
by the pigeonhole principle in view of the boundedness of the left hand side of 
(p9|) and what has just been proven, there exists (exercise) a dyadic sequence 
TiTi for which the first statement of Theorem 14.11 holds, with r"^ replaced by 
T^^ . Finally, by Theorem 13. 2i we immediately (exercise; why?) remove the 
restriction to the dyadic sequence. 

We have thus obtained 

f J^Wn'^<Br-'( f J,f"'(^)n^+ f Wn'^ + {n,^/j)nA . 

JYIt \Jt=0 -'So j_2 / 

(50) 

The statement (I50p loses one power of r in comparison with the first state- 
ment of Theorem 13.21 How do we obtain the full result? First of all, note that, 
commuting once again with flj, it follows that (I50p holds for ip replaced with 
fljip. Now we may partition TZ{0,t) dyadically into subrcgions TZ{Ti,Ti+i) and 



*^Hint: Use l|44|) to estimate the energy on {t = to} n J"'^(St) with weights in r. Send 
to -> cx) and estimate backwards to T,t using conservation of the flux. 
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revisit the X-estimate (|39p on each such region. In view of (|50p apphed to both 
ip and ^jip, the estimate (|39p gives 

f xJ^ii;)nl<BDr.\ (51) 

where Z? is a quantity coming from data. Summing over i, this gives us that 

f txJ!:'Wn'^<BD{l + \og\T + l\) 
Jiz(o,t) ^ 

This estimates in particular the first term on the right hand side of the first 
inequahty of (|48p . Applying this inequality, we obtain as before (|49p . but 
with T~^(l + log|T+ 1|) replacing r. Using (|5ip and a pigeonhole principle, one 
improves this to (|50p . with t^^(1 + logjr + 1|) now replacing r. iterating this 
argument again one removes the log (exercise). 

Note that this loss of derivatives in (I3ip simply arises from the loss in (|39p . 
If f2i could be replaced by fif in (I39p . then the loss would be 3e. The latter 
refinement can in fact be deduced from the original (I3ip using in addition work 
of Blue-Soffer |21| . Running the argument of this section with the e-loss version 
of (|3ip . we obtain now 

Proposition 4.2.1. For any e > 0, statement {\31^ holds with 3 replaced by e 
in the definition of Ei and C replaced by . 

4.3 Point wise decay 

To derive pointwise decay for ij: itself, we should remember that we have in fact 
dropped a good O'th order term from the estimate (I44p . In particular, we have 
also 

/ Jf-(^)n^>& / (rV^ + l)^'^. 

Jt=T J {t=T}n{r>ro} 

From this and the previously derived bounds, pointwise decay can be shown 
easily by applying VLt as commutators and Sobolev estimates. See [65] for details. 

Exercise: Derive pointwise decay for all derivatives of '0, including trans- 
verse derivatives to the horizon of any order, by commuting in addition with Y 
as in the proof of Theorem 13.21 

4.4 Comments and further reading 
4.4.1 The X-estimate 

The origin of the use of vector field multipliers of the type X (as in Section [47T|) 
for proving decay for solutions of the wave equation goes back to Morawetz. 
(These identities are generalisations of the classical virial identity, which has it- 
self a long and complicated history.) In the context of Schwarzschild black holes, 
the first results in the direction of such estimates were in Laba and Soffer [llOj 
for a certain "Schrodinger" equation (related to the Schwarzschild t- function). 
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and, for the wave equation, in Blue and Soffer [T^]. These results were incom- 
plete (see [20|), however, and the first estimate of this type was actually obtained 
in our [SS] , motivated by the original calculations of [TH IllOj . This estimate re- 
quired decomposition of ip into individual spherical harmonics and choosing 
the current J^ "^ separately for each ■0^. A slightly different approach to this 
estimate is provided by [501. ^ somewhat simpler choice of current J^ "^ which 
provides an estimate for all sufficiently high spherical harmonics was first pre- 
sented by Alinhac [T]. Our Section [4.1.11 is similar in spirit. The first estimate 
not requiring a spherical harmonic decomposition was obtained in [S7]. This 
is the current of Section [4.1.21 The problem of reducing the loss of derivatives 
in (|39p has been addressed in Blue-Soffer [lT]|f3 The results of [H] in fact also 
apply to the Reissner-Nordstrom metric. 

A slightly different construction of a current as in Section [4.1.2l has been given 
by Marzuola and collaborators [116j . This current does not require commuting 
with Vli. In their subsequent |115j . the considerations of [116j are combined 
with ideas from [651 167| to obtain an estimate which docs not degenerate on the 
horizon: One includes a piece of the current of Section 13.3.21 and exploits 
Proposition [331 

4.4.2 The Z-estimate 

The use of vector-field multipliers of the type Z also goes back to celebrated work 
of Morawetz, in the context of the wave equation outside convex obstacles [119] . 
The geometric interpretation of this estimate arose later, and the use of Z 
adapted to the causal geometry of a non-trivial metric first appears perhaps in 
the proof of stability of Minkowski space [51] . The decay result Theorem 14.11 
was obtained in our [65j . A result yielding similar decay away from the horizon 
(but weaker decay along the horizon) was proven independently in a nice paper 
of Blue and Sterbenz [52]. Both [55] and [BS] make use of a current based on 
the vector field Z. In [22], the error term analogous to K^'"' of Section [4?2l 
was controlled with the help of an auxiliary collection of multipliers with linear 
weights in t, chosen at the level of each spherical harmonic, whereas in [65| . these 
error terms are controlled directly from (|39p by a dyadic iteration scheme similar 
to the one we have given here in Section 14.21 The paper [55] does not obtain 
estimates for the non-degenerate energy flux (|3ip : moreover, a slower pointwise 
decay rate near the horizon is achieved in comparison to Theorem l4.1l Motivated 
by [65] , the authors of [22] have since given a different argument [23] to obtain 
just the pointwise estimate (I32p on the horizon, exploiting the "good" term in 
j^z,w jj^gg^j. |;jjg horizon. The proof of Theorem 14.11 presented in Section 14.21 is 
a slightly modified version of the scheme in |65] . avoiding spherical harmonic 
decompositions (for obtaining (|39p ) by using in particular the result of [67] . 

''■'A related refinement, where h of (1401) is replaced by a function vanishing logarithmically 
at 3Af, follows from |115j referred to below. 
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4.4.3 Other results 

Statement (|32p of Theorem 14.11 has been generaUsed to the Maxwell case by 
Blue [IH]. In fact, the Maxwell case is much "cleaner", as the current need 
not be modified by a function w, and its flux is pointwise positive through any 
spacelike hypersurface. The considerations near the horizon follow [23| and thus 
the analogue of (|3ip is not in fact obtained, only decay for the degenerate flux of 

. Nevertheless, the non-dcgcncratc (|3ip for Maxwell can be proven following 
the methods of this section, using in particular currents associated to the vector 
field Y (Exercise). 

To our knowledge, the above discussion exhausts the quantitative pointwise 
and energy decay-type statements which are known for general solutions of the 
wave equation on SchwarzschildF^ The best previously known results on gen- 
eral solutions of the wave equation were non-quantitative decay type statements 
which we briefly mention. A pointwise decay without a rate was first proven in 
the thesis of Twainy |149| . Scattering and asymptotic completeness statements 
for the wave, Klein-Gordon, Maxwell and Dirac equations have been obtained 
by [111 [73l m m 1122] . These type of statements are typically insensitive to the 
amount of trapping. See the related discussion of Section 14. 6[ where the state- 
ment of Theorem 14.11 is compared to non-quantitative statements heuristically 
derived in the physics literature. 

4.5 Perturbing? 

Use of the current "stabilises" the proof of Theorem 14.11 with respect to 
considerations near the horizon. There is, however, a sense in which the above 
argument is still fundamentally attached to Schwarzschild. The approach taken 
to derive the multiplier estimate (I36p depends on the structure of the trapping 
set, in particular, the fact that trapped null geodesies approach a codimension-1 
subset of spacetime, the photon sphere. Overcoming the restrictiveness of this 
approach is the fundamental remaining difficulty in extending these techniques 
to Kerr, as will be accomplished in Section [5?3l Precise implications of this fact 
for multiplier estimates are discussed further in [T] . 

4.6 Aside: Quantitative vs. non- quantitative results and 
the heuristic tradition 

The study of wave equations on Schwarzschild has a long history in the physics 
literature, beginning with the pioneering Reggc and Wheeler |131| . These stud- 
ies have all been associated with showing "stability" . 

A seminal paper is that of Price |130| . There, insightful heuristic arguments 
were put forth deriving the asymptotic tail of each spherical harmonic tpi evolv- 
ing from compactly supported initial data, suggesting that for r > 2M , 

Mr,t) - Cet-^^^^'\ (52) 

*^For fixed spherical harmonic £ = 0, there is also the quantitative result of | 63| . to be 
mentioned in Section 14.61 
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These arguments were later extended by Gundlach et al [5S] to suggest 

i^t\'H--Ctv-^^^^'\ r^,\x. - Ctu-^^^'\ (53) 

Another approach to these heuristics via the analytic continuation of the Green's 
function was followed by |3T]. The latter approach in principle could perhaps 
be turned into a rigorous proof, at least for solutions not supported on T-i^ nT-L" . 
See [m [106] for just 1^ for the £ = case. 

Statements of the form (I52p are interesting because, if proven, they would 
give the fine structure of the tail of the solution. However, it is important to 
realise that statements like (|52p in of themselves would not give quantitative 
bounds for the size of the solution at all later times in terms of initial data. In 
fact, the above heuristics do not even suggest what the best such quantitative 
result would be, they only give a heuristic lower bound on the best possible 
quantitative decay rate in a theorem like Theorem 14. II 

Let us elaborate on this further. For fixed spherical harmonic, by compact- 
ness a statement of the form (|52p would immediately yield some bound 

\M{r,t) < D{r,i^,)t-\ (54) 

for some constant D depending on r and on the solution itself. It is not clear, 
however, what the sharp such quantitative inequality of the form (|54p is sup- 
posed to be when the constant is to depend on a natural quantity associated 
to data. It is the latter, however, which is important for the nonlinear stability 
problem. 

There is a setting in which a quantitative version of (I54p has indeed been 
obtained: The results of |63j (which apply to the nonlinear problem where the 
scalar field is coupled to the Einstein equation, but which can be specialised to 
the decoupled case of the ^ = harmonic on Schwarzschild) prove in particular 
that 

|"s>o| + |^o|<C,W,^')t"'^', \r^|^o\<CD{^\>,^\>')T-^ (55) 

where depends only on e, and _D(i|),t|)') is a quantity depending only on the 
initial J"^ energy and a pointwise weighted norm. In view of the relation 
between r, u, and u, (|55p includes also decay on the horizon and null infinity 
as in the heuristically derived (|53p . The fact that the power 3 indeed appears 
in both in the quantitative (I55p and in (I54p may be in part accidental. See 
also [15] . 

For general solutions, i.e. for the sum over spherical harmonics, the situation 
is even worse. In fact, a statement like (|52p a priori gives no information 
whatsoever of any sort, even of the non-quantitative kind. It is in principle 
compatible with limsup^^^ ^'(''j ^) = oo|3 It is well known, moreover, that to 
understand quantitative decay rates for general solutions, one must quantify 
trapping. This is not, however, captured by the heuristics leading to 
essentially because for fixed £, the effects of trapping concern an intermediate 

*^Of course, given the quantitative result of Thcorem l3.2l and the statement (I52II . one could 
then infer that for each r > 2M, then limt->oo <f){r,t) = 0, without however a rate (exercise). 
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time interval not reflected in the tail. It should thus not be surprising that these 
heuristics do not address the fundamental problem at hand. 

Another direction for heuristic work has been the study of so-called quasi- 
normal modes. These are solutions with time dependence e~"^* for uj with 
negative imaginary part, and appropriate boundary conditions. These occur 
as poles of the analytic continuation of the resolvent of an associated elliptic 
problem, and in the scattering theory literature are typically known as reso- 
nances. Quasinormal modes are discussed in the nice survey article of Kokkotas 
and Schmidt |104] . Rigorous results on the distribution of resonances have been 
achieved in Bachelot-Motet-Bachelot [7] and Sa Barreto-Zworski |135| . The 
asymptotic distribution of the quasi-normal modes as £ ^ oo can be thought 
to reflect trapping. On the other hand, these modes fo not reflect the "low- 
frequency" effects giving rise to tails. Thus, they too tell only part of the story. 
See, however, the case of Schwarzschild-de Sitter in Section [6l 

Finally, we should mention Stewart |144j . This is to our knowledge the first 
clear discussion in the physics literature of the relevance of trapping on the 
Schwarzschild metric in this context and the difference between quantitative 
and non-quantitative decay rates. It is interesting to compare Section 3 of [144j 
with what has now been proven: Although the predictions of |144| do not quite 
match the situation in Schwarzschild (it is in particular incompatible with (|52p ) , 
they apply well to the Schwarzschild-de Sitter case developed in Section [HI 

The upshot of the present discussion is the following: Statements of the 
form (|52p . while interesting, may have little to do with the problem of non- 
linear stability of black holes, and are perhaps more interesting for the lower 
bounds that they sugg est0 In fact, in view of their non-quantitative nature, 
these results are less relevant for the stability problem than the quantitative 
boundedness theorem of Kay and Wald. Even the statement of Section 13.2.31 
cannot be derived as a corollary of the statement (|52p , nor would knowing ()52|) 
simplify in any way the proof of Section 13.2.31 

5 Perturbing Schwarzschild: Kerr and beyond 

We now turn to the problem of perturbing the Schwarzschild metric and prov- 
ing boundedness and decay for the wave equation on the backgrounds of such 
perturbed metrics. Let us recall our dilemma: The boundedness argument of 
Section [3] required that T remains causal everwhere in the exterior. In view 
of the comments of Section 13.51 this is clearly unstable. On the other hand, 
the decay argument of Section 3] requires understanding the trapped set and in 
particular, uses the fact that in Schwarzschild, a certain codimcnsion-1 subset 
of spacetime-the photon sphere-plays a special role. Again, as discussed in 
Section [4751 this special structure is unstable. 

It turns out that nonetheless, these issues can be addressed and both bound- 
edness (see Theorem 15. ip and decay (see Theorem 15. 2p can be proven for the 
wave equation on suitable perturbations of Schwarzschild. As we shall see, the 

*^See for instance the relevance of this in |59j . 
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boundedness proof (See Section 15. 2p turns out to be more robust and can be 
applied to a larger class of metrics-but it too requires some insight from the 
Schwarzschild decay argument! The decay proof (See Section l5.3p will require 
us to restrict to exactly Kerr spacetimes. 

Without further delay, perhaps it is time to introduce the Kerr family. . . 

5.1 The Kerr metric 

The Kerr metric is a 2-parameter family of metrics first discovered |99] in 1963. 
The parameters are called mass M and specific angular momentum a, i.e. an- 
gular momentum per unit mass. In so-called Boycr-Lindquist local coordinates, 
the metric element takes the form: 




The vector fields dt and 9^ arc Killing. We say that the Kerr family is stationary 
and axisvmmetricW\ Traditionally, one denotes 

A = r^-2A/r + a^. 

If a = 0, the Kerr metric clearly reduces to Schwarzschild ([5]). 

Maximal extensions of the Kerr metric were first constructed by Carter [29] . 
For parameter range < |a| < M, these maximal extensions have black hole 
regions and white hole regions bounded by future and past event horizons 'H'"' 
meeting at a bifurcate sphere. The above coordinate system is defined in a 
domain of outer communications, and the horizon will correspond to the limit 
r ^ r+, where r+ is the larger positive root of A = 0, i.e. 

r+ = M + \/M2-a2. 

Since the motivation of our study is the Cauchy problem for the Einstein 
equations, it is more natural to consider not maximal extensions, but maximal 
developments of complete initial data. (See Appendix |B]) In the Schwarzschild 
case, the maximal development of initial data on a Cauchy surface S as described 
previously coincides with maximally-extended Schwarzschild. In Kerr, if we are 
to take an asymptotically flat (with two ends) hypersurface in a maximally 
extended Kerr for parameter range < \a\ < M , then its maximal development 
will have a smooth boundary in maximally-extended Kerr. This boundary is 

**There are various conventions on the meaning of the words "stationary" and "axisym- 
metric" depending on the context. Let us not worry about this here. . . 
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what is known as a Cauchy horizon. We have aheady discussed this phenomenon 
in Section 12.7.31 in the context of strong cosmic censorship. The maximally 
extended Kerr solutions are quite bizarre, in particular, they contain closed 
timelikc curves. This is of no concern to us here, however. By definition, for us 
the term "Kerr metric (Al, (7A/,a)" will always denote the maximal development 
of a complete asymptotically flat hypersurface E, as above, with two ends. One 
can depict the Penrose-diagramatic representation of a suitable two-dimensional 
timelike slice of this solution as below: 




This depiction coincides with the standard Penrose diagram of the spherically 
symmetric Reissner- Nordstrom metric [911 1148] . 

With this convention in mind, we note that the dependence of gM,a on a 
is smooth in the range < |a| < M. In particular, Kerr solutions with small 
\a\ « M can be viewed as close to Schwarzschild. 

One can see this explicitly in the subregion of interest to us by passing to a 
new system of coordinates. Define 

t* =t + t{r) 

(/)* = (/) + (j){r) 



where 



^"(r) = (r2 + a^)/A^ ^(r) = a/A. 
ar dr 



(These coordinates are often known as Kerr- star coordinates.) These coor- 
dinates are regular across \ We may finally define a coordinate 
'"Schw = ''Schw(7'7a) such that which takes [r+,oo) [2A/, oo) with smooth de- 
pendence in a and such that rschwir, 0) is the identity map. In particular, if we 
define Sq by r> = {t* = 0}, and define 11 = V n {t* > 0}, and fix rschw, t* , ((>* 



^^Of course, one again needs two coordinate systems in view of the breakdown of spherical 
coordinates. We shall suppress this issue in the discussion that follows. 
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Schwarzschild coordinates, then the metric functions of gM,a written in terms 
of these coordinates as defined previously depend smoothly on a for < \a\ < M 
in 7?., and, for a = 0, reduce to the Schwarzschild metric form in {r,t*,(j>,6) 
coordinates where t* is defined from Schwarzschild t as above. 

Wc note that dt = dt* in the intersection of the coordinate systems. We 
immediately note that dt is spacelike on the horizon, except where = 0,7r, 
i.e. on the axis of symmetry. Note that we shall often abuse notation (as we 
just have done) and speak of dt on the horizon or at = 0, where of course 
the (r, t, 9, (f) coordinate system breaks down, and formally, this notation is 
meaningless. 

In general, the part of the domain of outer communications plus horizon 
where dt is spacelike in known as the ergoregion. It is bounded by a hypersurface 
known as the ergosphere. The ergosphere meets the horizon on the axis of 
symmetry = 0, tt. 

The ergosphere allows for a particle "process" , originally discovered by Pen- 
rose [127] ■ for extracting energy out of a black hole. This came to be known as 
the Penrose process. In his thesis, Christodoulou |38| discovered the existence 
of a quantity-the so-called irreducible mass of the black hole-which he should 
to be always nondecreasing in a Penrose process. The analogy between this 
quantity and entropy led later to a subject known as "black hole thermodynam- 
ics" [HUH]- This is currently the subject of intense investigation from the point 
of view of high energy physics. 

In the context of the study of Dgtp = 0, we have already discussed in Section[5] 
the effect of the ergoregion: It is precisely the presence of the ergoregion that 
makes our previous proof of boundedness for Schwarzschild not immediately 
generalise for Kerr. Moreover, in contrast to the Schwarzschild case, there 
is no "easy result" that one can obtain away from the horizon analogous to 
Scction [3.2.3l In fact, the problem of proving any sort of boundedness statement 
for general solutions to Ugip = on Kerr had been open until very recently. We 
will describe in the next section our recent resolution |68| of this problem. 

5.2 Boundedness for axisymmetric stationary black holes 

Wc will derive a rather general boundedness theorem for a class of axisymmetric 
stationary black hole exteriors near Schwarzschild. The result (Theorem 15. ip 
will include slowly rotating Kerr solutions with parameters \a\ « M . 

Wc have already explained in what sense the Kerr metric is "close" to 
Schwarzschild in the region TZ. Let us note that with respect to the coordinates 
fSchw, t* , 0*, in TZ, then dt* and d^* are Killing for both the Schwarzschild 
and the Kerr metric. The class of metrics which will concern us here are met- 
rics defined on TZ such that the metric functions are close to Schwarzschild in a 
suitable sensj^. and dt*^ dip* are Killing, where these are defined with respect 
to the ambient Schwarzschild coordinates. 

There is however an additional geometric assumption we shall need, and this 

^''This requires moving to an auxiliary coordinate system. See |68j . 
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is motivated by a geometric property of the Kerr spacetime, to be described in 
the section that follows immediately. 

5.2.1 Killing fields on the horizon 

Let us here remark a geometric property of the Kerr spacetime itself which turns 
out to be of utmost importance in what follows: Let V denote a null generator 
of 7^+. Then 

F€Span{at»,a^»}. (56) 

There is a deep reason why this is true. For stationary black holes with 
non-degenerate horizons, a celebrated argument of Hawking retrieves a second 
Killing field in the direction of the null generator V. Thus, if df and 9^* span 
the complete set of Killing fields, then V must evidently be in their span. 
In fact, choosing V accordingly we have 

V = dt* + ia/2Mr^)d^. (57) 

(For the Kerr solution, we have that there exists a timelike direction in the 
span of dt* and d^* for all points outside the horizon. We shall not explicitly 
make reference to this property, although in view of Section [3 one can infer 
this property (exercise) for small perturbations of Schwarzschild of the type 
considered here, i.e., given any point p outside the horizon, there exists a Killing 
field V (depending on p) such that V{p) is timelike.) 

5.2.2 The axisymmetric case 

From (|57p . it follows that there is a constant wq > 0, depending only on the 
parameters a and M, such that if 

|9t.Vl'>t^o|90*^|', (58) 

on TL'*^ , then the flux satisfies 

Jj(^)n^, > 0. (59) 
Note also that, for fixed M, we can take 

Wo ^ 0, as a ^ 0. (60) 

There is an immediate application of (|58p . Let us restrict for the moment to 
axisymmetric solutions, i.e. to V' such that d^ip = 0. It follows that (|58p trivially 
holds. As a result, our argument proving boundedness is stable, i.e. Theorem l3.2l 
holds for axisymmetric solutions of the wave equation on Kerr spacetimes with 
\a\ « M. (See the exercise of Section [5751 ) In fact, the restriction on a can be 
be removed (Exercise, or go directly to Section [7]). 

Let us note that the above considerations make sense not only for Kerr but 
for the more general class of metrics on TZ close to Schwarzschild such that dt» , 
d(f,* are Killing, T-L^ is null and (|56p holds. In particular, (|58|) implies (I59p . 
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where in (|60p , the condition a is replaced by the condition that the metric 
is taken suitably close to Schwarzschild. The discussion which follows will 
refer to metrics satisfying these assumptions For simplicity, the reader 
can specialise the discussion below to the case of a Kerr metric with \a\ « M. 

5.2.3 Superradiant and non-superradiant frequencies 

There is a more general setting where we can make use of (|58|) . Let us suppose 
for the time being that we could take the Fourier transform V'(^) of our solution 
ijj in t* and then expand in azimuthal modes tpm, i-C. modes associated to the 
Killing vector field dcf,* . 

If we were to restrict tJj to the frequency range 

\uj\'^ > Uorr? , (61) 

then (|58|) and thus (|59p holds after integrating along Ti^ . In view of this, 
frequencies in the range (|6ip are known as nonsuperradiant frequencies. The 
frequency range 

\ujf < uorr? (62) 

determines the so-called superradiant frequencies. In the physics literature, the 
main difficulty of this problem has traditionally been perceived to "lie" with 
these frequencies. 

Let us pretend for the time being that using the Fourier transform, we could 
indeed decompose 

= "0(1 + V'b (63) 

where tpj is supported in (|61|) . whereas tjj^ is supported in (I62p . 

In view of the discussion immediately above and the comments of Sec- 
tion 15.2.21 it is plausible to expect that one could indeed prove boundedness 
for ^/jjj in the manner of the proof of Theorem 13.21 In particular, if one could 
localise the integrated version of (|59p to arbitrary sufficiently large subsegments 
?^(r',r"), one could obtain 

f jr^(^p)<^ < B f jm^px,- (64) 

This would leave Since this frequency range does not suggest a direct 
boundedness argument, it is natural to revisit the decay mechanism of Schwarz- 
schild. We have already discussed (see Section |43)) the instability of the decay 
argument; this instability arose from the structure of the set of trapped null 
geodesies. At the heuristic level, however, it is easy to see that, if one can take 
ujQ sufficiently small, then solutions supported in (|62p cannot be trapped. In 
particular, for \a\ « M. superradiant frequencies for Ug'il' = on Kerr 
are not trapped. This will be the fundamental observation allowing for the 
boundedness theorem. Let us sec how this statement can be understood from 
the point of view of energy currents. 

^^They are summarised again in the formulation of Theorem 15.11 
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5.2.4 A stable energy estimate for superradiant frequencies 

We continue here our heuristic point of view, where we assume a decomposition 
(|63p where ip\, is supported in (|62|) . In particular, one has an inequahty 

/OO /^27r OO 27r 

L0l{d^4,^f d(l>* dt* > J {dt^^f dcj)* dt* (65) 

for aU (r, 9). We shaU see below that (I65p allows us easily to construct a suitable 
stable current for Schwarzschild. 

It may actually be a worthwhile exercise for the reader to come up with 
a suitable current for themselves. The choice is actually quite flexible in com- 
parison with the considerations of Section [iTT] Our choice (see [55]) is defined 

by 

where here, N is the vector field of Section r3. 3. 2[ Xa = fad,-*, with 

fa = -r'^{ro)^, forr<ro 
fa = -1, forro<r<i?i, 

fa = -1+ f 7^ for Ri<r<R2, 
fa = 0forr>i?2, 

Xh = fbdr* with 

/h = X(r-*)7r"^ f " 2 
Jo a;^ + a'' 

and xC''*) is a smooth cutoff with x = for r* < and x = 1 for r* > 1. Here r 
and r* arc Schwarzschild coordinates^ The function Wh is given by 

Wb = fl, + -{l-2M/r){l-M/r)fi,. 
r 

The parameters e, a, tq, i?2 must be chosen accordingly! 

Restricting to the range (|62p , using (|65p , with some computation we would 
obtain 

/OO 27r OO 27r 

K'^{i,^)d<P*de>bj_^j^ Xj;^{4'^)n'^d<p*dt\ (67) 

for all {r,e). 

The above inequality can immediately be seen to be stable to smallfl ax- 
isymmetric, stationary perturbations of the Schwarzschild metric. That is to 

Since we are dealing now with general perturbations of Schwarzschild, wc shall now use 
r for what we previously denoted by rgchw Note that in the special case that our metric is 
Kerr, this r is different from the Boyer-Lindquist r. 

^•^Of course, in view of the degeneration towards , it is important that smallncss is under- 
stood in a weighted sense. 
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say, for such metrics, if tjj^ is supported in (|62|) (where frequencies here are de- 
fined by Fourier transform in coordinates t* , (j)*), then the inequaUty (|67p holds 
as before. In particular, (|67p holds for Kerr for small \a\ « M. 

How would (|67p give boundcdncss for i/j^l We need in fact to suppose some- 
thing slightly stronger, namely that (|67p holds localised to 72.(0, r). Consider 
the currents 

where 62 is a positive parameters, and is the current of Section [B. 3. 21 Then, 
for metrics g close enough to Schwarzschild, and for 62 sufficiently small, we 
would have from a localised (|57|) that 

1(0,.) ^W^O, 

and thus 

Moreover, for g sufficiently close to Schwarzschild and ei, 62 suitably defined, 
we also have (exercise) 

f Jll^- {i,^)n^^ < B f JM\>X^- 

We thus would obtain 

[ jT^{i^^)n^<B f JrnV^b)"^ (68) 

Adding (fSS]) and (f64|) . we would obtain 

provided that we could also estimate say 

[ jr°(V'H)n^<S [ (69) 

5.2.5 Cutoff and decomposition 

Unfortunately, things arc not so simple! 

For one thing, to take the Fourier transform necessary to decompose in 
frequency, one would need to know a priori that ip{t*,-) is in L^{t*). What we 
want to prove at this stage is much less. A priori, can grow exponentially in 
t* . In order to apply the above, one must cut off the solution appropriately in 
time. 

This is achieved as follows. For definiteness, define Sq to be t* = 0, and E,- 
as before. We will also need two auxiliary families of hypcrsurfaccs defined as 
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follows. (The motivation for considering these will be discussed in Section [5.2.6n 
Let X be a cutoff such that xi^) = for x > and x = 1 for x < -1, and define 
by 

= t* - x{-r + R){1 + r - R)^'^ 

and 

r = t* + x{-r + R){1 + r - R)^'^ 

where i? is a large constant, which must be chosen appropriately. Let us define 
then 

Finally, we define 

n{T,,T2)= U s(t), 

n <r<T2 

n <T<T2 

7^"(n,T2)= u s-(t). 

n <r<r2 

Let ^ now be a cutoff function such that ^ = 1 in J^(I]j) n J"(I]^_j), and 
^ = in J+(S+) n J"(Eo). We may finally define 

= ^V- 

The function ?/'5.< is a solution of the inhomogencous equation 

Note that F is supported in W uTV {t - 1,t) . 

Another problem is that sharp cutoffs in frequency behave poorly under 
localisation. We thus do the following: Let C be a smooth cutoff supported in 
[-2, 2] with the property that C = f in [-1, f ], and let > be a parameter to 
be determined later. For an arbitrary ^' of compact support in i*, define 

*a(r,-) = *o+ E / (l-C((^om)-ia;)) #™(a;,.)e'-* do;. 



Note of course that ^fji + = ^. Wc shall use the notation ijjt^ for {■4}^)\^ and 
■i/'jj for (V's<)jj- Note that Vn satisfy 

n<,V^l, = ^^l„ n3Vp = i^p. (70) 
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5.2.6 The bootstrap 

With ^jj^, ip^ well defined, we now try to fill in the argument heuristically outlined 
before. 

We wish to show the boundedness of 

qi sup f J^n^'. (71) 

0<-f<r J^f 

We will argue by continuity in r. We have already seen heuristically how to 
obtain a bound for q in Sections 15.2.31 and 15.2.41 When interpreted for the , 
defined above, these arguments produce error terms from: 

• the inhomogeneous terms Fi,, _Fjj from (ITOj) 

• the fact that we wish to localise estimates (I59p and (I65p to subregions 
H+(T',r") and TZir' ,t") resepectively 

• the fact that (|69p is not exactly true. 

These error terms can be controlled by q itself. For this, one studies carefully 
the time-decay of , Fa away from the cutoff region 7?,"(0, 1)u7?,^(t-1, t) using 
classical properties of the Fourier transform. An important subtlety arises from 
the presence of O'th order terms in and it is here that the divergence of the 
region TV" from 7?.(0, r) is exploited to exchange decay in r and r. 

To close the continuity argument, it is essential not only that the error terms 
be controlled by q itself, but that a small constant is retrieved, i.e. that the error 
terms are controlled by eq, so that they can be absorbed. For this, use is made 
of the fact that for metrics in the allowed class sufficiently close to Schwarzschild 
(in the Kerr case, for \a\ « M), one can control a priori the exponential growth 
rate of (I7ip to be small. See |68) . 

5.2.7 Pointwise bounds 

Having proven the uniform boundedness of (I7ip , one argues as in the proof of 
Theorem l3.2l to obtain higher order energy and pointwise bounds. In particular, 
the positivity property in the computation of Proposition 13.3.21 is stable. (It 
turns out that this positivity property persists in fact for much more general 
black hole spacetimes and there is in fact a geometric reason for this! See 
Chapter El) 

5.2.8 The boundedness theorem 

We have finally 

Theorem 5.1. Let g be a metric defined on the differ entiahle manifold TZ with 
stratified boundary T-C" u Sq, and let T and $ = 17i be Schwarzschild Killing 
fields. Assume 

1. g is sufficiently close to Schwarzschild in an appropriate sense 
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2. T and $ are Killing with respect to g 

3. l-L^ is null with respect to g and T and $ span the null generator of T-L^ . 
Then the statement of Theorem \3.2\ holds. 

See [SS] for the precise formulation of the closeness assumption [TJ 

Corollary 5.1. The result applies to Kerr, and to the more general Kerr- 
Newman family (solving Einstein-Maxwell), for parameters \a\ « M (and also 
\Q\ « M in the Kerr-Newman case). 

Thus, we have quantitative pointwisc and energy bounds for ip and arbitrary 
derivatives on slowly rotating Kerr and Kerr-Newman exteriors. 

5.3 Decay for Kerr 

To obtain decay results analogous to Theorem 14.11 one needs to understand 
trapping. For general perturbations of Schwarzschild of the class considered 
in Theorem 15.11 it is not a priori clear what stability properties one can infer 
about the nature of the trapped set, and how these can be exploited. But 
for the Kerr family itself, the trapping structure can easily be understood, in 
view of the complete integrability of geodesic flow discovered by Carter [55] . The 
codimensionality of the trapped set persists, but in contrast to the Schwarzschild 
case where trapped null geodesies all approach the codimensional-1 subset r = 
SM of spacetime, in Kerr, this codimensionality must be viewed in phase space. 

5.3.1 Separation 

There is a convenient way of doing phase space analysis in Kerr spacetimes, 
namely, as discovered by Carter [3D], the wave equation can be separated. 
Walker and Penrose [153] later showed that both the complete integrability 
of geodesic flow and the separability of the wave equation have their fundamen- 
tal origin in the presence of a Killing tensor^^ In fact, as we shall see, in view of 
its intimate relation with the integrability of geodesic flow. Carter's separation 
of Ug immediately captures the codimensionality of the trapped set. 

The separation of the wave equation requires taking the Fourier transform, 
and then expanding into oblate spheroidal harmonics. As before, taking the 
Fourier transform requires cutting off in time. We shall here do the cutoff, 
however, in a somewhat different fashion. 

Let Et- be defined specifically as t* = r. Given r' < r, define TZ{t',t) as 
before, and let ^ be a cutoff fimction as in Section [5. 2. 51 but with Er'+i replacing 
Ej, T,ri replacing Eq, and Et- replacing Y,^, Et-_i replacing E^_]^. Define as 
before 

= ^V- 

^*See [321 1108| for recent higher-dimensional generalisations of these properties. 
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The function ipi.< is a solution of the inhomogeneous equation 

Note that F is supported in TI{t',t' + 1) u 7^(r - l,r). 

Since t/'s^ is compactly supported in t* we may consider its Fourier transform 
= ip^(uj,-). We may now decompose 



rntp* 



ml 

where Smi are the oblate spheroidal harmonics. For each m e Z, and fixed lo, 
these are a basis of eigenfunctions Smi satisfying 

~~: — (sin6'— S'm^) + , ^ ^ Smi - a^w^ cos^ 03^1 = \niSmt, 
sm 9 dO \ dO I sm 9 

and, in addition, satisfying the orthogonality conditions with respect to the 9 
variable, 

J dip J d{cos9)e""''^ Sme{auJ,cos9) e'""- Srn'i'{aiu,cos9) = Srnm'Sw- 

Here, the Xme{'-^) are the eigenvalues associated with the harmonics Smi- Each 
of the functions R^iir) is a solution of the following problem 

Note that if a = 0, we typically label Smi by € > |m| such that 

A„,(^)=^(^+l)/2. 

With this choice, Smt coincides with the standard spherical harmonics Y^i- 

Given any wi > 0, Ai > then we can choose a such that for \uj\ < uji, 
Am£ < Ai, then 

\Xmi-i{i+l)/2\<e. 
Rewriting the equation for the oblate spheroidal function 



ld/.„rf„ \ "1 22 2/10 

smt*— + . y = Kii^rni +a cos 6'5,„f, 
sm 61 d6' V d6' / sm 6* 



the smallest eigenvalue of the operator on the left hand side of the above equation 
is m{m + 1). This implies that 

Ato£ > m(TO+ 1) - a^w^. (72) 

This will be all that we require about Xmi- For a more detailed analysis of Xmi, 
see EH. 
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5.3.2 Frequency decomposition 

Let C be a sharp cutoff function such that C = 1 for kl ^ 1 a-nd C = for |a;| > f . 
Note that 

= 1. (73) 

Let uji, Ai be (potentiaUy large) constants to be determined, and A2 be a (po- 
tentially small) constant to be determined. 
Let us define 



V'b= f^a^/^i) E R^,{r)S,ni{'^iu,cose)e'"'^*e''^''dLu, 

Pc(Wt^i) E i?:^,,(r)5,„,(ac^,cos0)e*'"^*e'"**dt^, 

/oo ^ ^ 

(l-ac^M)) E i?;;,(r)5™,(ao;,cos^)e™'^ e'^-* do;, 



^(l-aW^^i)) E KAr) SmAato, cos e)e"'''^ e'^' dio. 

Wc have clearly 

For quick reference, we note: 

• "01, is supported in < wi, Xmt < Ai, 

• Tpi is supported in < ui, Xmt > Ai, 

• tp^ is supported in > uji, X„ie > X2U!'^ and 

• -011 is supported in \uj\ > uii, Xmt < X2UJ^. 

5.3.3 The trapped frequencies 

Trapping takes place in -0^ . We show here how to construct a multiplier for this 
frequency range. 

Defining a coordinate r* by 

dr* + c? 



dr A 

and setting 

u{r) = (r^ + a^'^R-Ur), H{r) - ^^f^f^,, . 
then u satisfies 
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where 



^mAn - (r2 + a2)2 + 

Consider the foUowing quantity 



df ^ ( du \ 1 d^f , ,2 



2dr* 



Then, with the notation ' = 



Q' = 2f'\uf - fV'\u\^ + Rc{2fHu' + f'Hu) - 



For i/j^ , we have 



We set 



Xjni + i^'^c? > (A2 + a^)w^ > (A2 + a^)tJi. 
Vo = (Am£ + w a ) 



(74) 



(75) 



(r^ + a2)2 



so that 



Vi = V-Vo = 1-^ ^75 + - 

(r^ + a^)^ 

Using ([72]) . (TTS]) . we easily see that 



(r2+a2)3 (r2 + a2)^ 



r'\v;\^ 



\ Ar2 V 



< CAr ^ (|aTOti;| + a^{\,ni + o^lo"^) + l) 

< eAr-2(A„,, + aV), (76) 



where e can be made arbitrarily small, if uji is chosen sufficiently large, and a 
is chosen a < e. On the other hand 



y' - 2 

° (r2 + a2)4 

Ar'" 



(Am^ + i^^CL^ ) ((r - M)(r2 + a^) - 2r{r'^ - 2Mr)) 



^^(A„, + ^V)(r- 



(r2 + a2) 



3M + 



r-M' 



(77) 



This computation implies that Vq has a simple zero in the neighborhood of 
r = 3M. Furthermore, 



< -Ar ^(A,„<? + tj^a^) 



From the above and (|76p , it follows that for ivi sufficiently large and a sufficiently 
small, we have 



\ Ar2 / 



+ V'y<-iAr-2(A,„,+c.V). 
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This alone implies that V' has at most a simple zero. 

To show that V' indeed has a zero we examine the boundary values at r+ 
and oo. From {\77\i we see that 

— Vo ~ 6 {Xrni + u; a ) 

for some positive constant C on the horizon and 

near r = oo. On the other hand, from the inequality as applied to the first term 
on the right hand side of (f75|) . it follows that 



< tr{\rnl +UJ a ), 



where e can be chosen arbitrarily small if uji is chosen sufficiently large and a 
sufficiently small. Thus, for suitable choice of lui, it follows that 



Ar' 



2 



Ar^ loo Ar^ i 



and thus y has a unique zero. Let us denote the r- value of this zero by rj^^. 
We now choose / so that 

1. /'>o, 

2. / < for r < r'^^^ and / > for r > rj^^ , 
3- -fV - i/'" > c. 

Property [3] can be verified by ensuring that f"'{r':^i) < as well as requiring 
that /"' < at the horizon. We may moreover normalise / to -1 on the horizon. 
Finally, we may assume that there exists an R such that for all r > R, f is of 
the form: 

/ = tan-i r*-a-^ _ _ 
a. 

In particular, for r > R, the function / will not depend on uj, £, m. 

Note the similarity of this construction with that of Section I4.1.1[ modulo 
the need for complete separation to centre the function / appropriately. 

Integrating the identity ([71)1 and using that m ^ as r ^ oo we obtain that 
for any compact set Ki in r* and a certain compact set K2 (which in particular 
does not contain r = 3Af), there exists a positive constant > so that 

bf {\u'\'' + \u\'')dr + b{X^e + u;^) f 

<(|uf + (w2-T/)|w|2)(r+) + J Re{2 f Hu' + fiiu) dr. 
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On the horizon r = r+, we have v! = {iuj + {iam/2Mr+))u and 

AMramuj - a^m^ 



{r1 + a?Y 
Therefore, we obtain 

bf {\uf + \u\'')dr* + b{X^i + uj^)f \u\''dr* 

<{lo'^ + em^)\u\^{r^) + j Rc{2 f Hu' + f'Hu)dr* . (78) 

We now wish to reinstate the dropped indices m,£,a;, and sum over to, £ 
and integrate over uj. Note that by the orthogonahty of the S^^g, it foUows that 
for any functions a and /? with coefficients defined by 

a{LL>,-) = a^j^(r)S',„£(aa;,cos6')e™''^ , /3(w,-) = ^ P^^i{r)Smi{auj, cos 6) e^""^ 
we have 

f {t*,r, 6,^) sin edifi d9 dt* = f°°Y.\<ntir)\'^duj, 
/ a- Psin 9 dip d9 dt* = / ^ aj^^ • /J^^dw. 

m,i 

Clearly, the summed and integrated left hand side of (|78p bounds 

Similarly, we read off immediately that the first term on the right hand side of 
(|78p upon summation and integration yields precisely 



Note that we can bound 

< B f Jl^Wn^^ + e f (d,.^f 

(Exercise: Why?) 

The "error term" of the right hand side of (|78p is more tricky. To estimate 
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the second summand of the integrand, note that 

E S~'\{frniKef{r) + S\^^fdco 

^ 5-'B\F-^,\\r) + 5\r^,fdio 

m,^:Am£(cj)>A2CJ^ 

= r^B J (F^f sine d<j)dedt* + s J H'^f sine d(t)dedt* 

< J F^ sine d(j)dedt* +S J sin e d(j) de dt* , 

where 6 can be chosen arbitrarily. In particular, this estimate holds for r < R. 
For r > i?, in view of the fact that / is independent of w, m, £, we have in fact 

^ oo _ 

m , £ : A jyi, ^ ( ) > A 2 o; ^ 

/oo _ 
E F-Arm,{r)du; 

= f'{r) J F^iP^sined(t)dedt* 
= f'{r) J F^i]:^sined4)dedt* , 

where for the last line we have used (I73p . The first summand of the error 
integrand of (I78p can be estimated similarly. 
We thus obtain 

<B f J^{^)n^^ + ef + F^ 

+ 5 f i!^ + (drll}f 

Jnn{r<R} 

+ [°°dt* [ (2/(r2 + a2)i/2Ra..((r2+a2)i/2^^) 

J-oo Jr>R 4 

+/Yr^ + a^)RV's-<) sin e d(j)* dedr*, (79) 

where x is a cutoff which degenerates at infinity and ft, is a function < h < 1 
which vanishes in a suitable neighborhood of r = 3Af. 

5.3.4 The untrapped frequencies 

Given A2 sufRciently small and any choice of lji, Ai, then, for a sufficiently small 
(where sufficiently small depends on these latter two constants), it follows that 
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for = b, <i, tt, wc may produce currents of type ^ as in Section 15.2.41 such 
that 

for X 9- suitable cutoff function degenerating at infinity, and x ^ suitable cutoff 
function degenerating at infinity and vanishing in a neighborhood of 'H'^ . These 
currents can in fact be chosen independently of a for such small a, and moreover, 
they can be chosen so that, defining 

we have on the one hand 

Jnn{r>R} J-oo Jr>R. 

+f' {r^ + a^)F>,ij^)—^— sine d(t)* dOdr* 
for the / of Section fS. 3. 31 and on the other hand, for the region r < R, we have 
f £'^><B5-^ f F^ + B5 f it + idri'^f+xJUiiJ^h" 

Jnn{r<R} JHn{r<R] Jnn{r<R} ^ 

where x is supported near the horizon and away from a neighborhood of r = 3Af . 
Moreover, one can show as in Section [5 . 2 . 61 that 



< B 



(Exercise: Prove the last inequality.) 
From the identity 

and the above remarks, one obtains finally an estimate 

<B f jf^{i')nt' + BS-^ f 

JS^, ^ Jnn{r<R} 

+ BS f + (drtl^f + xJuMN>' 

Jnn{r<R} ^ 

+ f°°dt* f {2f{r' + ay{F^ + F^+F.)dA(r' + aY'^i,^) 

J-oo J r>R ^ 



+f'{r^ + a^){F^ +F^+ FM'^) smO dcf)* dO dr* . (80) 
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5.3.5 The integrated decay estimates 

Now, we will add (|79p . (|80p and the energy identity of eJ^{ip) 



J^-r Jn(T',T)n{r<ro} 



7Z(r',T)n{ro<ri<ro} 



(81) 



for a small e with e « e, and where ro < ri < 3Af are as in Corollary 13.11 and ri 
is in the support of 1^2- 

In the resulting inequality, the left hand side bounds in particular 



where x is a cutoff decaying at infinity, x is a cutoff decaying at infinity and 
vanishing at and h is a function with < h < 1 such that h vanishes precisely 
in a neighborhood of r = 3M. (As a ^ 0, this neighborhood can be chosen 
smaller and smaller in the sense of the coordinate r.) 

Let us examine the right hand side of the resulting inequality. 

The second term of the first line of the right hand side of (I79p is absorbed 
by the first term on the right hand side of (ISip provided that e « e. 

The third term of the first line of the right hand side of (I79p and the second 
term of (|50|) are bounded by 



in view of Theorem 15.11 

The second line of the right side of (|79p and the third term of (|80p can be 
absorbed by (|82p , provided that S is chosen suitably small, whereas the second 
term of the right hand side of (jSip can be absorbed by (f82)) , provided that e is 
sufficiently small. 

The fourth terms of the right hand sides of (|79p and (|80p combine to yield 



Note where F is supported and how it decays. Using our boundedness Theo- 
rem 15.11 a Hardy inequality and integration by parts we may now bound this 
term by 



But the remaining terms on the right hand side of (|79p , (|80p and (|8ip are also 
of this form! We thus obtain 




(82) 




J -oo Jr>R 

+f'{r'^+a^)Fi,^) sin6>d0* dO dr* . 



+ a 
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Proposition 5.3.1. There exists a ipt-invariant weight degenerating only at 
io, a second ipt -invariant weight x, degenerating at iq and vanishing at 'H^ , a 
third tpt-invariant weight h, which vanishes on a neighborhood of r = 3M, and 
a constant B > such that the following estimates hold for all r' < r, 

f Xhj^{m' + X^^<B i J.fWK, 

JTCyr'.r) Ji^^i ^ 

JTZ{t',t) jt.^, ^ 

for all solutions dgip = on Kerr. 

Similar estimates could be shown on regions 1Z{t' ,t), E^, after having de- 
rived a priori suitable decay of ip in rF^ 

5.3.6 The Z-estimate 

To turn integrated decay as in Proposition l5.3.1] into decay of energy and point- 
wise decay, we must adapt the argument of Section 

Let y be a c/it-invariant vector field such that V = dt* for r > r+ + C2 and 
V = dt* + {a/2Mr+)dif,* for / < r+ + ci for some ci < C2, and such that V is 
timelike in 7?. \ 'H'^ . Note that V is Killing except in r+ + ci < r < r+ + C2. As 
a -> 0, we can construct such a V with C2 arbitrarily small. 

Now let us define u and v to be the SchwarzschilcQ coordinates 

U-L - fSchwJ 

V = t + Tg^.!-,^. 

With respect to the coordinates defining L = du, then L vanishes 

smoothly along the horizon. Define L = V - L. Finally, define the vector field 

Z = u^L + v^L. 

Note that under these choices Z is null on "H^. With w as before, the 
currents J^'™ together with can be used to control the energy fiuxes on St 
with weights. Use of the energy identities of J^'™ and leads to estimates of 
the form 

f Xi^^+ f , J^{^)nl <BDT-^+Br-^ f £, (83) 

where x is a cutoff function supported suitably, and where £ is an error term 
arising from the part of K^''^ which has the "wrong" sign; D arises from data. 
We may partition 

£ = £i+£2+£3 

where 

^®In the section that follows, we shall in fact localise the above estimate in a different way 
applying a cutoff function. The resulting O'th order terms which arise can be controlled using 
the "good" O'th order term in the boundary integrals of J^'™. 

Recall that we are considering both the Kerr and Schwarzschild metric on the fixed dif- 
ferentiable structure Ti as described in Section l5.ll 
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• £\ is supported in some region < r < i?o, 

• El is supported in r < rp and 

• £3 is supported in r > i?o- 

Recall that L+L is Killing for r > 2M+C2- It follows (Exercise) that choosing 
C2 <tq^ there are no terms growing quadratically in t for £1, £3. Moreover, by our 
construction, Z depends smoothly on a away from the horizon. The behaviour 
near the horizon is more subtle as Z itself is not smooth! We shall return to 
this when discussing £2- 

In view of our above remarks, we have that 

fi<i?t(J^(^)7V^ + V'), 

just like in the case of Schwarzschild. In view of Proposition 1 5 . 3 . l1 this leads to 
the following estimate: li ij^ = ij} in TI{t',t") n {r < Rq}, where -0 solves again 
□gV' = 0, then 

f £,iyj)= f SM)<Bt' f iJ^W + j;!im)nl . (84) 

Jn{T',T") Jtz{t',t") je^, ^ ^ 

The introduction of tl) is related to our localisation procedure we shall carry out 
in what follows. 

Recall that in the Schwarzschild case, for Rq suitably chosen, there is no 
£3 term, as the term X^'™ has a good sign in that region. (Sec Section 221) 
Examining the r-decay of error terms in the smooth dependence of Z in a, we 
obtain 

where e can be made arbitrarily small if a is small. If r" - t' ^ t' ~ t, this leads 
to an estimate 

f £M) < e{r"-T')iT" + r') f -J^^Wn^, 

Jn{T',T") JT.^,n{r<T"-T'} ^ 

+ e\og\T"-T'\ f J,r(V)<^,. (85) 

In the region r+ + ci < r < r+ + C2, then, choosing vq such that £2 is absent in 
Schwarzschild, we can argue without computation from the smooth dependence 
on a that 

£2<et\j;^{^P)N^ + ^P') 

where e can be made arbitrarily small by choosing a small. The necessity of a 
quadratically growing error term arises from the fact that i + i is not Killing 
in this region IfJ 

As we have already mentioned, an important subtely occurs near the horizon 
where Z fails to be . This means that £2 is not necessarily small in local 

coordinates, and one must understand how to bound the singular terms. It 

turns out that these singular terms have a structure: 

^^Alternatively, one can keep L + L Killing at the expense of Z failing to be causal on the 
horizon. This would lead to errors of a similar nature. 
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Proposition 5.3.2. Let V, Y, Ei, E2 extend V to a null frame in r < r+ + ci. 
We have 

£2 < ev\\og{r-r^)\P{T{Y,V) + T{V,V))+evJl^{'^)N^'. 

Proof. The warping function w can be chosen as in Schwarzschild near , and 
thus, the extra terms it generates are harmless. For the worst behaviour, it 
suffices to examine now K'^ itself. We must show that terms of the form: 

|log(r-r,)r(T(f,y)) 

do not appear in the computation for . 

The relevant property follows from examining the covariant derivative of Z 
with respect to the null frame: 

VyZ = 2u{Vu)L + 2vV{v)L + v^VyV - Ar^vVyL + 4:{r*fVyL, 

VyZ = 2u{Yu)L + 2v{Yv)L + v^Vy^ - 4r*vVyL + 4(r*)2VyL, 
VeiZ = 2u{Eiu)L + 2v{Eiv)L + v^V EiV - 4r*vy EiL + 4:{r*fv EiL, 
VE2Z = 2u{E2u)L + 2v{E2v)L + v^V e^V - 'ir*vV E2L + A{r*fv E2L. 

□ 

To estimate now £2, we first remark that with Proposition I5.3.T1 we can 
obtain the following refinement of the red-shift multiplier construction of Corol- 
lary O 

Proposition 5.3.3. If we weaken the requirement that N be smooth in Corol- 
larv \3.1\ with the statement that N is at TV' and smooth away from T-L^ , then 
given p > 0, we may construct an N as in Corollary \S.1\ where property [7] is 
replaced by the stronger inequality: 

K^'i^) > 6p| log(r - r,)r(T(y, V) + T{V, V)) 

for r< ro- 
ll now follows immediately from Proposition 15.3.11 that with and as 
before, we have 

To obtain energy decay from (fSS]) . ([M]) and dHS]), we argue now by 

continuity. Introduce the bootstrap assumptions 



L 



j';{i^)N>' + xr <C Dt-^^'\ (87) 



{T^)N^' < C' Dt-^^^" (88) 
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for a, S > 0. 

Now dyadically decompose the interval [0,r] by < r^+i. Using ([84|) and 
the above, we obtain 

f £i(tlj)<y f £iU) 

i 

<6-\CDt-^^^^ + C'Dt^^). (89) 

Here, ip is constructed separately on each dyadic region TZ{Ti,Ti+i) by throwing 
a cutoff on i/^ls^. eqaul to 1 in r < ti+i - Ti and vanishing in Ti+i - Ti < r, 
solving again the initial value problem in TZ(Ti,Ti+i), and exploiting the domain 
of dependence property. See the original |65| for this localisation scheme. The 
parameters of the "dyadic" decomposition must be chosen accordingly for the 
constants to work out. Similarly, using (|86p we obtain 



f f2(^)<E f ^2(^) 

j J^rin{r<Ti + i-Ti} 
i 

< eS-\^^CD (90) 

and using (|85p 

<eY.{r^Tf--^'CD^D\ogT') 

i 

< e5-^T^^CD. (91) 



For Tip we obtain 
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Y -/S,^n{r<r,+i-r,} 

i 

< e5-^T^^^^C'D + eS-'^T^^CD, (93) 
Jn{o.T) i Jn{r,,Ti^i) 

<eY{r!rr^^^'C'D + D\ogn) 

i 

< e6-\^^^^C"D. (94) 

Wc use here the algebra of constants where Be = e. The constant D is a quantity 
coming from data. Exercise: What is D and why is ([92]) true? 

For e « 6 and C sufficiently large, wc sec that from (I83p applied to Tip in 
place of ip, using (l92|) . (l93|) . wc improve (l88|) . 

On the other hand choosing C" « C and then r sufficiently large, we have 

t-^5-\CDt-^^^' + C'Dt^') < ^CDt-^^^' 

and thus, again for e « 5, using (|89p . (|90p wc can improve (|87p from (I83p . 

Once one obtains (|87p . then decay can be extended to decay in St- by the 
argument of Section [42l by applying conservation of the flux backwards!^ 

5.3.7 Pointwise bounds 

In any region r < _R, wc may now obtain pointwise decay bounds simply by 
further commutation with T, N as in Section [5.3.4l To obtain the correct point- 
wise decay statement towards null infinity, one must also commute the equation 
with a basis fli for the Lie algebra of the Schwarzschild metric, exploiting the 
r-weights of these vector fields. Defining Hi = C{r)ni, where C is a cutoff which 
vanishes for r < Rq, where 3M « Rq, and, setting ip = rtip, we have 

Ogip = Fid'^i; + F2dip 

**Note that in view of the fact that we argued by continuity to obtain (j87|l . we could not 
obtain this extended decay through earlier. This is why we have localised as in 1651 . not 
as in Section 14.21 
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where Fi = 0{r~^) and F2 = 0{r~^). Having estimates already for ip, Tip, one 
can may apply the X and Z estimates as before for tp, only, in view of the F2 
term, now one must exploit also the X-estimate in n {?■ < n+i - Ti}) n 

J"(S^;^j). Wc leave this as an exercise. 

5.3.8 The decay theorem 

Wc have obtained thus 

Theorem 5.2. Let {J\4,g) be Kerr for \a\ « M , T> he the closure of its domain 
of dependence, let Eg ^6 surface T) n {t* = 0}, let ^\), ^\)' be initial data 
on Eq swc/i that '\\> 6 H^^^{T.), \|)' e 7Jj-^"^^(E) for s > 1, and hm^-^io ijj = 0, 
and let ip be the corresponding unique solution of Ogip = 0. Let denote the 
1-parameter family of diffeomorphisms generated by T , let Eq be an arbitrary 
spacelike hypersurface in J^(Eo n where lA is an open neighborhood of the 
asymptotically fiat ewtF^. and define Tit = <Pt(So). Let s > 3 and assume 

-/So 

Then there exists a S > depending on a ( with S^Oasa^O) and a B 
depending only on Eq such that 

JSt S^ 
Now let s > 5 and assume 

E2^T. E f r2(j;<'(r» + j;o(r"TV') + j;°(r"rTV))<<oo 

\a\<2r={T.N.ni} -^^0 

where fli are the Schwarzschild angular momentum operators. Then 
sup\/f|i/'| < B\fE'2T'^^\ suprlV'l < Bv^-^^"^^''^^^- 

Onc can obtain decay for arbitrary derivatives, including transversal deriva- 
tives to H^, using additional commutation by N . See [5^ . 



5.4 Black hole uniqueness 

In the context of the vacuum equations ([3]) , the Kerr solution plays an important 
role not only because it is believed to be stable, but because it is believed to 
be the only stationary black hole solutionis This is the celebrated no-hair 
"theorem". In the case of the Einstein-Maxwell equations, there is an analogous 
no-hair "theorem" stating uniqueness for Kerr-Newman. A general reference 
is [92]. 

^^This is just the assumption that Sq "terminates" on null infinity 

^"A further extrapolation leads to the "belief" that all vacuum solutions eventually decom- 
pose into n Kerr solutions moving away from each other. 
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Neither of these results is close to being a theorem in the generality which 
they are often stated. Reasonably definitive statements have only been proven 
in the much easier static case, and in the case where axisymmetry is assumed 
a priori and the horizon is assumed connected, i.e. that there is one black hole. 
Axisymmetry can be inferred from stationarity under various special assump- 
tions, including the especially restrictive assumption of analyticity. See [57] for 
the latest on the analytic case, and [96] for new interesting results in the di- 
rection of removing the analyticity assumption in inferring axisymmetry from 
stationarity. 

Nonetheless, the expectation that black hole uniqueness is true reasonably 
raises the question: why the interest in more general black holes, allowed in 
Theorem [53? 

For a classical "astrophysical" motivation, note that black hole solutions can 
in principle exist in the presence of persistent atmospheres. Perhaps the simplest 
such constructions would involve solutions of the Einstein- Vlasov system, where 
matter is described by a distribution function on phase space invariant under 
geodesic flow. These black hole spacetimes would in general not be Kerr even 
in their vacuum regions. Recent speculations in high energy physics yield other 
possible motivations: There are now a variety of "hairy black holes" solving 
Einstein-matter systems for non-classical matter, like Yang-Mills fields |141| . 
and a large variety of vacuum black holes in higher dimensions [77] , many of 
which are currently the topic of intense study. 

There is, however, a second type of reason, which is relevant even when we 
restrict our attention to the vacuum equations (|4]) in dimension 4. The less 
information one must use about the spacetime to obtain quantitative control 
on fields, the better chance one has at obtaining a stability theorem. The 
essentially non-quantitativj^ aspect of our current limited understanding of 
black hole uniqueness should make it clear that these arguments probably will 
not have a place in a stability proof. Indeed, it would be an interesting problem 
to explore the possibility of obtaining a more quantitative version of uniqueness 
theorems (in a neighborhood of Kerr) following ideas in this section. 

5.5 Comments and further reading 

Theorem lS.ll was proven in [68j . In particular, this provided the first global result 
of any kind for general solutions of the Cauchy problem on a (non-Schwarzschild) 
Kerr background. Theorem 15 . 21 was first announced at the Clay Summer School 
where these notes were lectured. Results in the direction of Propositio n [ 5.3.11 
are independently being studied in work in progress by Tataru-TohaneaniJ^i and 



The best previous results concerning Kerr had been obtained by Finster and 
collaborators in an important series of papers culminating in [79| . See also |80| . 

^^As should be apparent by the role of analyticity or Carloman estimates, 
^^communication from Mihai Tohaneanu, a summer school participant who attended these 
lectures 

^^lecture of P. Blue, Mittag-Leffler, September 2008 
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The methods of [75] are spectral theoretic, with many pretty applications of 
contour integration and o.d.e. techniques. The results of [79] do not apply to 
general solutions of the Cauchy problem, however, only to individual azimuthal 
modes, i.e. solutions -0™ of fixed to. In addition, [73] imposes the restrictive 
assumption that T-L^ n not be in the support of the modes. (Recall the 
discussion of Section [3.2.61 ) Under these assumptions, the main result stated 
in [79] is that 

limVm(r,0 = (95) 

for any r > r+. Note that the reason that (|95p did not yield any statement 
concerning general solutions, i.e. the sum over m-not even a non-quantitative 
one-is that one did not have a quantitative boundedness statement as in Theo- 
rem l5.ll Moreover, one should mention that even for fixed m, the results of [TS] 
are in principle compatible with the statement 

supV'm = oo, 

i.e. that the azimuthal modes blow up along the horizon. See the comments in 
Section 14.61 It is important to note, however, that the statement of [79] need 
not restrict to |a| « M , but concerns the entire subextremal range \a\ < M . 
Thus, the statement (|95|) of [75] is currently the only known global statement 
about azimuthal modes on Kerr spacetimes for large but subextremal a. 

There has also been interesting work on the Dirac equation [751 1^0], for 
which superradiance does not occur, and the Klein-Gordon equation |89] . For 
the latter, see also Section 1531 

5.6 The nonlinear stability problem for Kerr 

We have motivated these notes with the nonlinear stability problem of Kerr. 
Let us give finally a rough formulation. 

Conjecture 5.1. Let {Yi,g,K) he a vacuum initial data set (see Avvendix \B.S\) 
sujficiently close (in a weighted sense) to the initial data on Cauchy hypersurface 
in the Kerr solution (A^,(7A/.a 

) for some parameters < |a| < M. Then the 
maximal vacuum development {M.,g) possesses a complete null infinity T'*' such 
that the metric restricted to J^(I^) approaches a Kerr solution (A1, gA/^.a^ ) 
in a uniform way ( with respect to a foliation of the type Sr of Section [^ with 
quantitative decay rates, where Mj, aj are near M , a respectively. 

Let us make some remarks concerning the above statement. Under the 
assumptions of the above conjecture, {M,g) certainly contains a trapped surface 
S by Cauchy stability. By Penrose's incompleteness theorem (Theorem l2.2|) . this 
implies that {Ai,g) is future causally gcodesically incomplete. By the methods 
of the proof of Theorem l2.21 it is easy to see that 5n J"(I^) = 0. Thus, as soon 
as is shown to be complete, it would follow that the spacetimc has a black 
hole region in the sense of Section I2.5.4(^ 

®^Let us also remark the obvious fact that the above conjecture impUes in particular that 
weak cosmic censorship holds in a neighborhood of Kerr data. 
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In view of this, one can also formulate the problem where the initial data 
are assumed close to Kerr initial data on an incomplete subset of a Cauchy 
hypcrsurfaco with one asymptotically flat end and bounded by a trapped surface. 
This is in fact the physical problenF^. but in view of Cauchy stability, it is 
equivalent to the formulation we have given above. Note also the open problem 
described in the last paragraph of Section 12.81 

In the spherically symmetric analogue of this problem where the Einstein 
equations are coupled with matter, or the Bianchi-triaxial IX vacuum problem 
discussed in Section [2.6.41 the completeness of null infinity can be inferred easily 
without detailed understanding of the geometry [SOIISI] • One can view this as an 
"orbital stability" statement. In this spherically symmetric case, the asymptotic 
stability can then be studied a posteriori, as in [63l[94]. This latter problem is 
much more difficult. 

In the case of Conjecture 15.11 in contrast to the symmetric cases mentioned 
above, one does not expect to be able to show any weaker stability statement 
than the asymptotic stability with decay rates as stated. Note that it is only the 
Kerr family as a whole-noi the Schwarzschild subfamily-which is expected to 
be asymptotically stable: Choosing a = certainly does not imply that a/ = 0. 
On the other hand, if \a\ « M, then by the formulation of the above conjecture, 
it would follow that \af\ « Mf. It is with this in mind that we have considered 
the \a\ « M case in this paper. 



6 The cosmological constant A and Schwarzschild- 
de Sitter 

Another interesting setting for the study of the stability problem arc black 
holes within cosmological spacetimes. Cosmological spacetimes-as opposed to 
asymptotically flat spacetimes (See Appendix IB.2.3P , which model spacetime 
in the vicinity of an isolated self-gravitating system-are supposed to model 
the whole universe. The working hypothesis of classical cosmology is that the 
universe is approximately homogeneous and isotropic (sometimes known as the 
Copernican principle [91|). In the Newtonian theory, it was not possible to 
formulate a cosmological model satisfying this hypothesis Iff] One of the major 
successes of general relativity was that the theory allowed for such solutions, 
thus making cosmology into a mathematical science. 

In the early years of mathematical cosmology, it was assumed that the uni- 
verse should be staticFI. To allow for such static cosmological solutions, Einstein 

^^Cf. the comments on the relation between maximally-extended Schwarzschild and 
Oppenheimer-Snyder. 

^^It is possible, however, if one geometrically reinterprets the Newtonian theory and al- 
lows space to be— say-the torus. See 11321 . These reinterpretations, of course, postdate the 
formulation of general relativity. 

^^much like in the early studies of asymptotically fiat spacetimes discussed in Section 12.11 
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modified his equations 1^) by adding a O'th order term: 

Rtj.iy - ^df^i^R + ^9tii^ = 8'!TTf_,^. (96) 

Here A is a constant known as the cosmological constant. When coupled with a 
perfect fluid, this system admits a static, homogeneous, isotropic solution with 
A > and topology S'^ x R. This spacetime is sometimes called the Einstein 
static universe. 

Cosmological solutions with various values of the parameter A were studied 
by Friedmann and Lcmaitre, under the hypothesis of exact homogeneity and 
isotropy. Static solutions arc in fact always unstable under perturbation of ini- 
tial data. Typical homogeneous isotropic solutions expand or contract, or both, 
beginning and or ending in singular configurations. As with the early studies 
(referred to in Sections 12. 2p illuminating the extensions of the Schwarzschild 
metric across the horizon, these were ahead of their timcF^ (See the forth- 
coming book [123j for a history of this fascinating early period in the history 
of mathematical cosmology.) These predictions were taken more seriously with 
Hubble's observational discovery of the expansion of the universe, and the sub- 
sequent evolutionary theories of matter, but the relevance of the solutions near 
where they are actually singular was taken seriously only after the incomplete- 
ness theorems of Penrose and Hawking-Penrosc were proven (see Section [2 .71) . 

We shall not go into a general discussion of cosmology here, nor tell the 
fascinating story of the ups and downs of A--from its adoption by Einstein to 
his subsequent well-known rejection of it, to its later "triumphant" return in 
current cosmological models, taking a very small positive value, the "explana- 
tion" of which is widely regarded as one of the outstanding puzzles of theoretical 
physics. Rather, let us pass directly to the object of our study here, one of the 
simplest examples of an inhomogeneous "cosmological" spacetime, where non- 
trivial small scale structure occurs in an ambient expanding cosmology. This is 
the Schwarzschild-de Sitter solution. 

6.1 The Schwarzschild-de Sitter geometry 

Again, this metric was discovered in local coordinates early in the history of 
general relativity, independently by Kottler |105| and Weyl |155j . Fixing A > oF^ 
Schwarzschild-de Sitter is a one-parameter family of solutions of the from 

- (1 - 2M/r - Ar^)dt^ + (1 - 2M/r - Ar^y^dr^ + r^dasi. (97) 

®*In fact, the two are very closely related! The interior region of the Oppenheimer-Snyder 
collapsing star is precisely isometric to a region of a Friedmann universe. See |117) . 

®^The expression (|97|l with A < defines Schwarzschild-anti-de Sitter. See Section 18.41 
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The black hole case is the case where < M < A maximally-extended 

solution (see PSIISS] ) then has as Penrose diagram the infinitely repeating chain: 



r = 




To construct "cosmological solutions" one often takes spatially compact quo- 
tients. (One can also glue such regions into other cosmological spacctimcs. 
See [S3]. For more on the geometry of this solution, see [TU].) 

6.2 Boundedness and decay 

The region "analogous" to the region studied previously for Schwarzschild and 
Kerr is the darker shaded region 2? above. The horizon T-l separates I? from an 
"expanding" region where the spacetimc is similar to the celebrated de-Sitter 
space. If S is a Cauchy surface such that E n = E n H =0, then let 
us define Eg = 2? n E, and let us define E,- to be the translates of Eq by the 
flow ipt generated by the Killing field T (= ^). Note that, in contrast to the 
Schwarzschild or Kerr case, Eg is compact. 



Theorem 6.1. The statement of Theorem \3. Si holds for these spacetimes, where 
E, Eq, Et- are as above, and lim^^jo is replaced by sup^^j^^^ \\\)\. 

Proof. The proof of the above theorem is as in the Schwarzschild case, except 
that in addition to the analogue of iV, one must use a vector field N which plays 
the role of N near the "cosmological horizon" . It is a good exercise for the 
reader to think about the properties required to construct such a iV. A general 
construction of such a vector field applicable to all non-extremal stationary black 
holes is done in Section [T) □ 

As for decay, wc have 

Theorem 6.2. For every k > 0, there exist constants Ck such that the following 
holds. Let ^\) e H^^^ , \|;' e Hy^^, and define 



We have 



(c,)\<kv={n,} -'^o 




Then 




(98) 
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For k > 1 we have 

sup IV' - 1 < Ck^Ei^T^ , (99) 

where tpQ denotes the 'th spherical harmonic, for which we have for instance 
the estimate 

sup iV'ol < sup ^\)o + Co^Eai-^o,-^'^). (100) 

The proof of this theorem uses the vector fields T, Y and Y (alternatively N, 
N), together with a version of X as multipliers, and requires commutation of the 
equation with Qi to quantify the loss caused by trapping. (Like Schwarzschild, 
the Schwarzschild-de Sitter metric has a photon sphere which is at r = 3A/ for 
all values of A in the allowed range. See [55] for a discussion of the optical 
geometry of this metric and its importance for gravitational lensing.) An esti- 
mate analogous to ()39|) is obtained, but without the x weight, in view of the 
compactness of Sq. The result of the Theorem follows essentially immediately, 
in view of Theorem 16. II and a pigeonhole argument. No use need be made of a 
vector field of the type Z as in Section Note that for tp = constant, Ek = 0, 
so removing the O'th spherical harmonic in (|99|) is necessary. See for details. 

Note that if fli can be replaced by fij in (|39p . then it follows that the loss 
in derivatives for energy decay at any polynomial rate k in (|98p can be made 
arbitrarily small. If fli could be replaced by logQi, then what would one obtain? 
(Exercise) 

It would be a nice exercise to commute with Y as in the proof of Theo- 
rem l6.1[ to obtain pointwise decay for arbitrary derivatives of k. See the related 
exercise in Section 14.31 concerning improving the statement of Theorem 14.11 

6.3 Comments and further reading 

Theorem 16.21 was proven in [66j . Independently, the problem of the wave equa- 
tion on Schwarzschild-de Sitter has been considered in a nice paper of Bony- 
Hafner [21] using methods of scattering theory. In that setting, the presence of 
trapping is manifest by the appearance of resonances, that is to say, the poles 
of the analytic continuation of the resolvent 1^ The relevant estimates on the 
distribution of these necessary for the analysis of [24] had been obtained earlier 
by Sa Barreto and Zworski [135j . 

In contrast to Theorem 16.21 the theorem of Bony-Hafner [24] makes the 
familiar restrictive assumption on the support of initial data discussed in Sec- 
tion [5231 For these data, however, the results of [23] obtain better decay than 
Theorem 16.21 awav from the horizon, namely exponential, at the cost of only an 
e derivative. The decay results of [23] degenerate at the horizon, in particular, 
they do not retrieve even boundedness for ip itself. However, using the result 
of [23] together with the analogue of the red-shift Y estimate as used in the 
proof of Theorem 16.21 one can prove exponential decay up to and including 

'^'^In the physics hterature, these are known as quasi-normal modes. See 11041 for a nice 
survey, as well as the discussion in Section 14.61 
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the horizon, i.e. exponential decay in the parameter r (Exercise). This still 
requires, however, the restrictive hypothesis of [24j concerning the support of 
the data. It would be interesting to sort out whether the restrictive hypothesis 
can be removed from [23] , and whether this fast decay is stable to perturbation. 
There also appears to be interesting work in progress by Sa Barreto, Melrose 
and Vasy |150| on a related problem. 

One should expect that the statement of Theorem 16.11 holds for the wave 
equation on axisymmetric stationary perturbations of Schwarzschild-de Sitter, 
in particular, slowly rotating Kerr-de Sitter, in analogy to Theorem 15. II 

Finally, we note that in many context, more natural than the wave equation 
is the conformally covariant wave equation Ogip - ^Rip = 0. For Schwarzschild- 
de Sitter, this is then a special case of Klein-Gordon (|106|) with /i > 0. The 
analogue of Theorem 16.11 holds by virtue of Section 17.21 Exercise: Prove the 
analogue of Theorem 16.21 for this equation. 

7 Epilogue: The red-shift effect for non-extremal 
black holes 

We give in this section general assumptions for the existence of vector fields 
Y and N as in Section 13.3.21 As an application, we can obtain the bounded- 
ness result of Theorem 13.21 or Theorem 16.11 for all classical non-extremal black 
holes for general nonnegative cosmological constant A > 0. See [HIl 11481 for 
discussions of these solutions. 

7.1 A general construction of vector fields Y and 

Recall that a Killing horizon is a null hypcrsurfacc whose normal is Killing [92[ 
I148j . Let "H be a sufficiently regular Killing horizon with (future-directed) gen- 
erator the Killing field V, which bounds a spacctimc V. Let (pY denote the 
one-parameter family of transformations generated by V, assumed to be glob- 
ally defined for all t > 0. Assume there exists a spatial hypersurface S c D 
transverse to V, such that nH = S is a compact 2-surface. Consider the 
region 

n' = ut>o(^r(i]) 

and assume that 7?.' n 2? is smoothly foliated by iptC^)- 
Note that 

VvV=kV 

for some function k ■■TL ^M.. 

Theorem 7.1. Let %, V, TV , S, V , ip( he. as above. Suppose k > 0. Then 
there exists a (pY -invariant future-directed timelike vector field N on TZ' and a 
constant b> such that 

in an open ipt-invariant (for t>0) subset U c TZ' containing "H n 72.'. 
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Proof. Define y on S' so that Y is future directed null, say 

g{Y,V) = -2, (101) 

and orthogonal to S. Moreover, extend F off 5 so that 

VyY = -a{Y +V) (102) 

on S. Now push Y forward by (pY to a vector field on U. Note that all the 
above relations still hold on 

It is easy to see that the relations (|19p - (|22p hold as before, where Ei, E2 are 
a local frame for TpLp({S). Now , are not necessarily 0, hence our having 
included them in the original computation! We define as before 

N =V + Y. 

Note that it is the compactness of S which gives the uniformity of the choice of 
6 in the statement of the theorem. □ 

We also have the following commutation theorem 

Theorem 7.2. Under the assumptions of the above theorem, ifip satisfies Hgip = 
0, then for all k > 1. 

0<|m|<fc, 0<m4<fe 

on T-L* , where Kk > 0- 

Proof. From (|19p -(|22 p . we deduce that relative to the null frame (on the hori- 
zon) V,Y, El, E2 the deformation tensor ^tt takes the form 

^TTyy = 2cr, ^TTvV = 2k, ^TTvY = C ^T^YEi = 0, ^T^VEi = a'l ^'^'EiEj = 

As a result the principal part of the commutator expression- the tcrm2'^7r"''i:)„i:);3V' 
can be written as follows 

2^7r"''VaV/3V = KVyy^J + Cr(Vyy + Vyy)V' " aVy^^^A + 2 /l} V s^. V' ■ 

The result now follows by induction on fc. □ 



7.2 Applications 

The proposition applies in particular to sub-extremal Kerr and Kerr-Newman, as 
well as to both horizons of sub-cxtrcmal Kerr-de Sitter, Kerr-Newman-de Sitter, 
etc. Let us give the following general, albeit somewhat awkward statement: 

Theorem 7.3. Let {TZ,g) be a manifold with stratified boundary Ti^ u S, such 
that TZ is globally hyperbolic with past boundary the Cauchy hypersurface 
where E andTi are themselves manifolds with (common) boundary S. Assume 
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where the unions are disjoint and each HI , Si is connected. Assume each 
satisfies the assumptions of Theorem \ 7. 1\ with future- directed Killing field Vi, 
some subset Si c S, and cross section a connected component Si of S . Let us 
assume there exists a Killing field T with future complete orbits, and ^pt is the 
one-parameter family of transformations generated by T . Let lAi be given by 
Theorem \7.1\ and assume that there exists a V as above such that 

7e = (^t(Ex V) u u,tiW«. 

and 

where Yi^- = (/^^-(S), ip^^ represents the one-parameter family of transformations 
generated by , and the last inequality is assumed for all values oft, r where 
the left hand side can be defined. Finally, let ip be a solution to the wave equation 
and assume that for any open neighborhood V of S in S, there exists a positive 
constant &v > such that 

Jl{T^i^)n^^>bvJ-HT^^)n'^ (103) 

m E \ V and 

= c,V^i}} (104) 

on "Hi . It follows that the first statement of Theorem \3.Si holds for ip. 

Assume in addition that S is compact or asymptotically flat, in the weak 
sense of the validity of a Sobolev estimate near infinity. Then the second 
statement of Theorem \3.2\ holds for -0. 

In the case where T is assumed timehke mTZ\ T-L^ , then (|104p is automatic, 
whereas (|103p holds if 

-g{T,T)>-bvg{n^,T) 

in E \ V. Thus we have 

Corollary 7.1. The above theorem applies to Reissner- Nordstrom, Reissner- 
Nordstrom-de Sitter, etc, for all subextremal range of parameters. Thus Theo- 
rem \3.S\ holds for all such metricsf^ 

On the other hand, (|104p . (|103p can be easily seen to hold for axisymmetric 
solutions tpo of DgV = on backgrounds in the Kerr family (see Section EH]) . 
We thus have 

Corollary 7.2. The statement of Theorem \3.S\ holds for axisymmetric solutions 
V'o of for Kerr-Newman and Kerr-Newman-de Sitter for the full subextremal 
range of parametersV^ 

Let us also mention that the the theorems of this section apply to the Klein- 
Gordon equation Ugi}} = IJ-'^ip, as well as to the Maxwell equations (Exercise). 

^^In the A = case this range is M > 0, < \Q\ < M. Exercise: What is it for A > 0? 

In the A - case this range is Af > 0, < IQI < n/M2 - a2 . Exercise: What is it for A > 0? 
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8 Open problems 



We end these notes with a discussion of open problems. Some of these are 
related to Conjecture 15. 1[ but all have independent interest. 

8.1 The wave equation 

The decay rates of Theorem 14.11 are sharp as uniform decay rates in v for any 
nontrivial class of initial data. On the other hand, it would be nice to obtain 
more decay in the interior, possibly under a stronger assumption on initial data. 

Open problem 1. Show that there exists a S > such that {\31\l holds with r 
replaced with t^^^^^^\ for a suitable redefinition of Ei. Show the same thing 
for Kerr spacetimes with \a\ « M . 

At the very least, it would be nice to obtain this result for the energy re- 
stricted to El- n {r < i?}. 

Recall how the algebraic structure of the Kerr solution is used in a funda- 
mental way in the proof of Theorem 15.21 On the other hand, one would think 
that the validity of the results should depend only on the robustness of the 
trapping structure. This suggests the following 

Open problem 2. Show the analogue of Theorem \5.2\ for the wave equation 
on metrics close to Schwarzschild with as few as possible geometric assumptions 
on the metric. 

For instance, can Theorem 15.21 be proven under the assumptions of Theo- 
rem [nnP Under even weaker assumptions? 

Our results for Kerr require \a\ « M . Of course, this is a "valid" assumption 
in the context of the nonlinear stability problem, in the sense that if this con- 
dition is assumed on the parameters of the initial reference Kerr solution, one 
expects it holds for the final Kerr solution. Nonetheless, one certainly would 
like a result for all cases. See the discussion in Section [5751 

Open problem 3. Show the analogue of Theorem 1 5. 2\ for Kerr solutions in the 
entire subextremal range < |a| < M . 

The extremal case \a\ = M may be quite different in view of the fact that 
Section [7] cannot apply: 

Open problem 4. Understand the behaviour of solutions to the wave equa- 
tion on extremal Reissner- Nordstrom, extremal Schwarzschild-de Sitter, and ex- 
tremal Kerr. 

Turning to the case of A > 0, we have already remarked that the analogue 
of Theorems 16.11 and 16.21 should certainly hold in the case of Kcrr-de Sitter. In 
the case of both Schwarzschild-de Sitter and Kcrr-de Sitter, another interesting 
problem is to understand the behaviour in the region C = J^iT-L^) ^ J^{'Hg), 
where "H^, "Wg are two cosmological horizons meeting at a sphere: 
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Open problem 5. Understand the behaviour of solutions to the wave equation 
in region C of Schwarzschild-de Sitter and Kerr-de Sitter, in particular, their 
behaviour along r = oo as is approached. 

Let us add that in the case of cosmological constant, in some contexts it is 
appropriate to replace Dg with the conformaUy covariant wave operator Dg-^R. 
In view of the fact that R is constant, this is a special case of the Klein-Gordon 
equation discussed in Section [8T3l below. 

8.2 Higher spin 

The wave equation is a "poor man's" linearisation of the Einstein equations 
(|4]). The role of linearisation in the mathematical theory of nonlinear partial 
differential equations is of a different nature than that which one might imagine 
from the formal "perturbation" theory which one still encounters in the physics 
literature. Rather than linearising the equations, one considers the solution of 
the nonlinear equation from the point of view of a related linear equation that 
it itself satisfies. 

In the case of the simplest nonlinear equations (say ()107p discussed in Sec- 
tion [8i6] below), typically this means freezing the right hand side, i.e. treating it 
as a given inhomogeneous term. In the case of the Einstein equations, the proper 
analogue of this procedure is much more geometric. Specifically, it amounts to 
looking at the so called Bianchi equations 

V[^i?,A]pa = 0, (105) 

which are already linear as equations for the curvature tensor when g is regarded 
as fixed. For more on this point of view, see [51]. The above equations for a 
field Sxfijyp with the symmetries of the Riemann curvature tensor are in general 
known as the spin- 2 equations. This motivates: 

Open problem 6. State and prove the spin-2 version of Theorems 15. Jl or \5.S\ 

( or Open problem [7p on Kerr metric backgrounds or more generally, metrics 
settling down to Kerr. 

In addition to [51], a good reference for these problems is [50], where this 
problem is resolved just for Minkowski space. In contrast to the case of Minkowski 
space, an additional difficulty in the above problem for the black hole setting 
arises from the presence of nontrivial stationary solutions provided by the cur- 
vature tensor of the solutions themselves. This will have to be accounted for 
in the statement of any decay theorem. From the "linearisation" point of view, 
the existence of stationary solutions is of course related to the fact that it is 
the 2-parameter Kerr family which is expected to be stable, not an individual 
solution. 

8.3 The Klein-Gordon equation 

Another important problem is the Klein-Gordon equation 

□ g'0 = /^V'- (106) 
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A large body of heuristic studies suggest the existence of a sequence of quasi- 
normal modes (see Section 14. 6p approaching the real axis from below in the 
Schwarzschild case. When the metric is perturbed to Kerr, it is thought that 
essentially this sequence "moves up" and produces exponentially growing solu- 
tions. See [1581 [7T] . This suggests 

Open problem 7. Construct an exponentially growing solution of (| 1 00^ on 
Kerr, for arbitrarily small ^ > and arbitrary small a. 

Interestingly, if one fixed m, then adapting the proof of Section [5. 2[ one can 
show that for /i > sufficiently small and a sufficiently small, depending on m, 
the statement of Theorem 15. II holds for (|106p for such Kerr's. This is consistent 
with the quasinormal mode picture, as one must take m ^ oo for the modes 
to approach the real axis in Schwarzschild. This shows how misleading fixed-m 
results can be when compared to the actual physical problem. 

8.4 Asymptotically anti-de Sitter spacetimes 

In discussing the cosmological constant we have considered only the case A > 0. 
This is the case of current interest in cosmology. On the other hand, from the 
completely different viewpoint of high energy physics, there has been intense 
interest in the case A < 0. See |84| . 

The expression (|97p for A < defines a solution known as Schwarzschild- 
anti-de Sitter. A Penrose diagramme of this solution is given below. 




The timelike character of infinity means that this solution is not globally hy- 
perbolic. As with Schwarzschild-de Sitter, Schwarzschild-anti-de Sitter can be 
viewed as a subfamily of a larger Kerr-anti de Sitter family, with similar prop- 
erties. 

Again, as with Schwarzschild-de Sitter, the role of the wave equation is in 
some contexts replaced by the conformally covariant wave equation. Note that 
this corresponds to (|106p with a negative fi = 2A/3 < 0. 

Even in the case of anti-de Sitter space itself (set M = in (jOT])). the question 
of the existence and uniqueness of dynamics is subtle in view of the timelike 
character of the ideal boundary I. It turns out that dynamics are unique for 
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(|106p only if the /i > 5A/12, whereas for the total energy to be nonnegative one 
must have /i > 3A/4. Under our conventions, the conformally covariant wave 
equation lies between these values. See [51 [^. 

Open problem 8. For suitable ranges of fi, understand the boundedness and 
blow-up properties for solutions of i\106'\i on Schwarzschild-anti de Sitter and 
Kerr-anti de Sitter. 

See [inni W\ for background. 

8.5 Higher dimensions 

All the black hole solutions described above have higher dimensional analogues. 
See [ZZl 1120] . These are currently of great interest from the point of view of 
high energy physics. 

Open problem 9. Study all the problems of Sections \8.1\\8.4\ in dimension 
greater than 4. 

Higher dimensions also brings a wealth of explicit black hole solutions such 
that the topology of spatial sections of T-L^ is no longer spherical. In particular, 
in 5 spacetime dimensions there exist "black string" solutions, and much more 
interestingly, asymptotically flat "black ring" solutions with horizon topology 
S^^S"^. See [77]. 

Open problem 10. Investigate the dynamics of the wave equation Ogip = 
and related equations on black ring backgrounds. 

8.6 Nonlinear problems 

The eventual goal of this subject is to study the global dynamics of the Einstein 
equations (HD themselves and in particular, to resolve Conjecture 15. II 

It may be interesting, however, to first look at simpler non-linear equations 
on fixed black hole backgrounds and ask whether decay results of the type proven 
here are suflacient to show non-linear stability. 

The simplest non-linear perturbation of the wave equation is 

ng^p = V'{^) (107) 

where V = V{x) is a potential function. Aspects of this problem on a Schwarz- 
schfld background have been studied by [T?I1 [Ml fTT5] . 

Open problem 11. Investigate the problem (\10T^ on Kerr backgrounds. 

In particular, in view of the discussion of Section 18.31 one may be able to 
construct solutions of (|107p with V = + IV'T: for M > and for arbitrarily 
large p, arising from arbitrarily small, decaying initial data, which blow up in 
finite time. This would be quite interesting. 
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A nonlinear problem with a stronger relation to (U) is the wave map problem. 
Wave maps are maps ^ : M. ^ N where M. is Lorentzian and N is Riemannian, 
which are critical points of the Lagrangian 

£($) = Jmi^ 

In local coordinates, the equations take the form 

where T^^ denote the Christoffel symbols of ^at. See the lecture notes of Struwe [145j 
for a nice introduction. 

Open problem 12. Show global existence in the domain of outer communi- 
cations for small data solutions of the wave map problem, for arbitrary target 
manifold M , on Schwarzschild and Kerr backgrounds. 

All the above problems concern fixed black hole backgrounds. One of the 
essential difficulties in proving Coniecturc l5.1l is dealing with a black hole back- 
ground which is not known a priori, and whose geometry must thus be recovered 
in a bootstrap setting. It would be nice to have more tractable model problems 
which address this difficulty. One can arrive at such problems by passing to 
symmetry classes. The closest analogue to Conjecture 15.11 in such a context is 
perhaps provided by the results of Holzegel [M], which concern the dynamic 
stability of the 5-dimensional Schwarzschild as a solution of (|4]), restricted un- 
der Triaxial Bianchi IX symmetry. See Section [2.6.41 In the symmetric setting, 
one can perhaps attain more insight on the geometric difficulties by attempting 
a large-data problem. For instance 

Open problem 13. Show that the maximal development of asymptotically flat 
triaxial Bianchi IX vacuum initial data for the 5-dimensional vacuum equations 
containing a trapped surface settles down to Schwarzschild. 

The analogue of the above statement has in fact been proven for the Einstein- 
scalar field system under spherical symmetry |401 163j . In the direction of the 
above, another interesting set of problems is provided by the Einstein-Maxwell- 
charged scalar field system under spherical symmetry. For both the charged- 
scalar field system and the Bianchi IX vacuum system, even more ambitious 
than Open problem 1 1 31 would be to study the strong and weak cosmic censorship 
conjectures, possibly unifying the analysis of [45l [58l [59] . Discussion of these 
open problems, however, is beyond the scope of the present notes. 
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A Lorentzian geometry 

The reader who wishes a formal introduction to Lorentzian geometry can con- 
sult |91j . For the reader familiar with the concepts and notations of Riemannian 
geometry, the following remarks should suffice for a quick introduction. 

A.l The Lorentzian signature 

Lorentzian geometry is defined as in Riemannian geometry, except that the 
metric g is not assumed positive definite, but of signature (-,+,...,+). That is 
to say, we assume that at each point p e A^"+^[^ we may find a basis of the 
tangent space TpAi, i = 0, . . . ,n, such that 

g = -Go <Si Go + ei <Si ei + •■■en ® +e„. 

In Riemannian geometry, the - in the first term on the right hand side would 
by +. 

A non-zero vector v e TpA4 is called timelike, spacelike, or null, according 
to whether g{v,v) < 0, g{v,v) > 0, or g{v,v) = 0. Null and timelike vectors 
collectively are known as causal. There are various conventions for the 0- vector. 
Let us not concern ourselves with such issues here. 

The appellations timelikc, spacelike, null are inherited by vector fields and 
immersed curves by their tangent vectors, i.e. a vector field V is timelike if 
V{p) is timelike, etc., and a curve 7 is timelike if 7 is timelike, etc. On the 
other hand, a submanifold T, c M is said to be spacelike if its induced geometry 
is Riemannian, timelike if its induced geometry is Lorentzian, and null if its 
induced geometry is degenerate. (Check that these two definitions coincide for 
embedded curves.) For a codimension-1 submanifold E c Ai, at every p 6 M, 
there exists a non-zero normal n'^, i.e. a vector in TpA4 such that g{n, v) = for 
all V € TpTi. It is easily seen that E is spacelike iff n is timelikc, E is timclike iff 
n is spacclike, and E is null iff n is null. Note that in the latter case n e TpT,. 
The normal of E is thus tangent to E. 

A. 2 Time-orientation and causality 

A time-orientation on {A4,g) is defined by an equivalence class [K] where 
K is a continuous timclike vector field, where Ki - K2 if g{Ki,K2) < 0. A 
Lorentzian manifold admitting a time-orientation is called time-orientable, and a 

'''^It is conventional to denote the dimension of the manifold by n+ 1. 
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triple {A4,g, [K]) is said to be a time-oriented Lorentzian manifold. Sometimes 
one reserves the use of the word "spacetime" for such triples. In any case, we 
shall always consider time-oriented Lorentzian manifolds and often drop explicit 
mention of the time orientation. 

Given this, we may further partition causal vectors as follows. A causal 
vector V is said to be future-pointing if g{v,K) < 0, otherwise past-pointing, 
where X is a representative for the time orientation. As before, these names 
are inherited by causal curves, i.e. we may now talk of a future- directed timelike 
curve, etc. Given p, we define the causal future J^{p) by 

J^{p) = p u {(7 € 7W : 37 : [0, 1] ^- : 7 future-pointing, causal} 

Similarly, we define J (p) where future is replaced by past in the above. If 
S c M is a set, then we define 

A. 3 Covariant derivatives, geodesies, curvature 

The standard local notions of Riemannian geometry carry over. In particular, 
one defines the Christoffcl symbols 

K\ = ^5^" (Si' 5a A + dxg^a - dagux), 
and geodesies 7(t) = (a;"(t)) are defined as solutions to 

+ T'^^x^'^x^ = 0. 

Here 17^1, denote the components of g with respect to a local coordinate system 
^/i^ gfif (^gQotes the components of the inverse metric, and we arc applying 
the Einstein summation convention where repeated upper and lower indices are 
summed. The Christoffcl symbols allow us to define the covariant derivative on 
{k,l) tensor fields by 

k I 
Ayi...Uk _ a Ayi...Vk , Vr^* Ayi...p...Vk _ V T'' 4'^i---!^fc 

1=1 1=1 

where it is understood that p replaces Ui, /i^, respectively in the two terms on 
the right. This defines (fc,Z + 1) tensor. As usual, if we contract this with a 
vector V at p, then we will denote this operator as Vt, and we note that this can 
be defined in the case that the tensor field is defined only on a curve tangent to 
V at p. We may thus express the geodesic equation as 

V^7 = 0. 

The Riemann curvature tensor is given by 

DP -/9,r'' - j_ pa pP _ pa pM 
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and the Ricci and scalar curvatures by 



Using the same letter R to denote aU these tensors is conventional in relativity, 
the number of indices indicating which tensor is being referred to. For this 
reason we will avoid writing "the tensor i?" . The expression R without indices 
will always denote the scalar curvature. As usual, we shall also use the letter R 
with indices to denote the various manifestations of these tensors with indices 
raised and lowered by the inverse metric and metric, e.g. 

RpuXp — 9 paRi/Xp 

Note the important formula 

Wc say that an immersed curve j ■ I ^ A4 is inextendible if there does not 
exist an immersed curve 7 : J ^ where J d I and 7I/ = 7. 

We say that {M,g) is geodesically complete if for all inextendible geodesies 
7 : / ^ A^, then / = R. Otherwise, we say that it is geodesically incomplete. 
We can similarly define the notion of spacelike geodesic (in)completeness, time- 
like geodesic completeness, causal geodesic completeness, etc, by restricting the 
definition to such geodesies. In the latter two cases, we may further specialise, 
e.g. to the notion of future causal geodesic completeness, by replacing the con- 
dition / = M with Id (a, 00) for some a. 

We say that a spacelike liypersurface c M is Cauchy if every inextendible 
causal curve in A4 intersects it precisely once. A spacetimc {M.,g) admitting 
such a hypersurface is called globally hyperbolic. This notion was first introduced 
by Leray [TT^ . 

B The Cauchy problem for the Einstein equa- 
tions 

We outline here for reference the basic framework of the Cauchy problem for 
the Einstein equations 

- Ig^^'^R + Ag^- = SttT^.. (io8) 

Here A is a constant known as the cosmological constant and T^i, is the so-called 
energy momentum tensor of matter. We will consider mainly the vacuum case 

Rp. = ^9pu, (109) 

where the system closes in itself. If the reader wants to set A = 0, he should feel 
free to do so. To illustrate the case of matter, we will consider the example of 
a scalar field. 
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B.l The constraint equations 

Let E be a spacelike hypersurface in {A4,g), with future directed unit timelike 
normal N. By definition, S inherits a Riemannian metric from g. On the 
other hand, we can define the so-called second fundamental form of E to be the 
symmetric covariant 2-tensor in TE defined by 

K{u,v) = -giVuV,N) 

where V denotes an arbitrary extension of v to a vector field along E, and V 
here denotes the connection of g. As in Riemannian geometry, one easily shows 
that the above indeed defines a tensor on TE, and that it is symmetric. 

Suppose now {A4,g) satisfies (jlOSp with some tensor T^^. With E as above, 
let gab, V, Kab denote the induced metric, connection, and second fundamental 
form, respectively, of E. Let barred quantities and Latin indices refer to tensors, 
curvature, etc., on E, and let n^(p) denote the components of the puUback map 
T*M ^ T*E. It follows that 

R + {K^y - K^Kl = 16^^Tf,,n^'n'' + 2A, (110) 

VbK^ - VaKl = le^n^T^.n'^. (Ill) 

To see this, one derives as in Riemannian geometry the Gauss and Codazzi 
equations, take traces, and apply (|108p . 

B.2 Initial data 

It is clear that (jllOp , (jllip are necessary conditions on the induced geometry of 
a spacelikc hypersurface E so as to arise as a hypersurface in a spacetimc satis- 
fying (jlOSp . As wc shall see, immediately, they will also be sufficient conditions 
for solving the initial value problem. 

B.2.1 The vacuum case 

Let E be a 3-manifold, g a Riemannian metric on E, and K a symmetric covari- 
ant 2-tensor. We shall call (E, g,i4r) a vacuum initial data set with cosmological 
constant A if ()110p - ()llip are satisfied with = 0. Note that in this case, 
equations (|110|) - (|111|) refer only to E, g, K. 

B.2. 2 The case of matter 

Let us here provide only the case for the Einstein-scalar field case. Here, the 
system is ()108p coupled with 

□gV = 0, (112) 
T^. = d^^d.tl; - ifiM-^V^VVa^. (113) 



97 



First note that were S a spacelike hypersurface in a spacetime {M.,g) satis- 
fying the Einstein-scalar field system with massless scalar field ip, and n'' were 
the future-directed normal, then setting I])' = n^d^cj), = we have 



-((ii;')' + V>V„4), 



where latin indices and barred quantities refer to S and its induced metric and 
connection. 

This motivates the following: Let S be a 3-manifold, g a Riemannian met- 
ric on S, K a symmetric covariant 2-tensor, and i|) : S ^ M, : S ^ M 
functions. We shall call {T,^g^K) an Einstein- scalar field initial data set with 
cosmological constant A if (|110|) - (|111|) are satisfied replacing Tf^^n^n'^ with 
+ V>Va^P), and replacing H^T^^n'^ with ^l^'Vai^- 

Note again that with the above replacements the equations (|110|) - (|111|) do 
not refer to an ambient spacetime . See |36j for the construction of solutions 
to this system. 

B.2.3 Asymptotic flatness and the positive mass theorem 

The study of the Einstein constraint equations is non-trivial! 

Let us refer in this section to a triple {Y,,g,K) where E is a 3-manifold, g 
a Riemannian metric, and K a symmetric two-tensor on S as an initial data 
set, even though we have not specified a particular closed system of equations. 
An initial data set (E,g,_fr) is strongly asymptotically flat with one end if there 
exists a compact set K. c Y, and a coordinate chart on E \ /C which is a diffeo- 
morphism to the complement of a ball in M^, and for which 



where 6ab denotes the Euclidean metric and r denotes the Euclidean polar co- 
ordinate. 

In appopriate units, M is the "mass" measured by asymptotic observers, 
when comparing to Newtonian motion in the frame 6ab- On the other hand, 
under the assumption of a global coordinate system well-behaved at infinity, M 
can be computed by integration of the tg component of a certain pseudotcnsoiQ 
added to Tq. In this manifestation, the quantity _E = M is known as the total 
energy^^ This relation was first studied by Einstein and is discussed in Wcyl's 

^*This is subtle: The Einstein vacuum equations arise from the Hilbert Lagrangian C(g) = 
/ R which is 2nd order in the metric. In local coordinates, the highest order term is a 
divergence, and the Lagrangian can thus be replaced by a new Lagrangian which is 1st order 
in the metric. The resultant Lagrangian density, however, is no longer coordinate invariant. 
The quantity 4g now arises from "Noether's theorem" |124) . See |49) for a nice discussion. 

'^^With the above asymptotics, the so-called linear momentum vanishes. Thus, in this case 
"mass" and energy are equivalent. 



gab = 




kab = Oi{r 
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book Raum-Zeit-Materie |154| . In one looks at E for a family of liypcrsurfaces 
with the above asymptotics, then E is conserved. 

A celebrated theorem of Schoen-Yau [1371 1138] (see also fl56j ) states 

Theorem B.l. Let {T,^g,K) be strongly asymptotically flat with one end and 
satisfy (f770|) . (f777|) with A = 0, and where i^n^n'^ , W^T^yU^ are replaced by 
the scalar /.j and the tensor Ja, respectively, defined on E, such that moreover 
H > \/J°\Ja- Suppose moreover the asymptotics are strengthened by replacing 
02(r-i) by Oiir-^) and oi{r-^) by Oair-^). Then M > and M = iff 12 
embeds isometrically into M'^^^ with induced metric g and second fundamental 
form K . 

The assumption ^ > ^J°Jq holds if the matter satisfies the dominant energy 
condition [91| . In particular, it holds for the Einstein scalar field system of 
Section IB.2.21 and (of course) for the vacuum case. The statement we have 
given above is weaker than the full strength of the Schoen-Yau result. For the 
most general assumptions under which mass can be defined, see [H]. 

One can define the notion of strongly asymptotically flat with k ends by 
assuming that there exists a compact K, such that E \ /C is a disjoint union of 
k regions possessing a chart as in the above definition. The Cauchy surface 
E of Schwarzschild of Kerr with < |a| < M, can be chosen to be strongly 
asymptotically flat with 2-ends. The mass of both ends coincides with the 
parameter M of the solution. 

The above theorem applies to this case as well for the parameter M associ- 
ated to any end. If M = for one end, then it follows by the rigidity statement 
that there is only one end. Note why Schwarzschild with A/ < does not provide 
a counterexample. 

The association of "naked singularities" with negative mass Schwarzschild 
gave the impression that the positive energy theorem protects against naked 
singularities. This is not true! See the examples discussed in Section [^.6.21 

In the presence of black holes, one expects a strengthening of the lower 
bound on mass in Theorem IB. II to include a term related to the square root of 
the area of a cross section of the horizon. Such inequalities were first discussed 
by Penrose |127| with the Bondi mass in place of the mass defined above. All 
inequalities of this type are often called Penrose inequalities. It is not clear what 
this term should be, as the horizon is only identifyable after global properties of 
the maximal development have been understood. Thus, one often replaces this 
area in the conjectured inequality with the area of a suitably defined apparent 
horizon. Such a statement has indeed been obtained in the so-called Riemannian 
case (corresponding to if = 0) where the relevant notion of apparent horizon 
coincides with that of minimal surface. See the important papers of Huisken- 
Ilmanen |95] and Bray [25] . 

B.3 The maximal development 

Let (E,5,iir) denote a smooth vacuum initial data set with cosmological con- 
stant A. We say that a smooth spacetime {M-^g) is a smooth development of 



99 



initial data if 

1. {Ai,g) satisfies the Einstein vacuum equations dJ]) with cosmological con- 
stant A. 

2. There exists a smooth embedding i : E such that {M.,g) is globally 
hyperbolic with Cauchy surface and g, K are the induced metric 
and second fundamental form, respectively. 

The original local existence and uniqueness theorems were proven in 1952 by 
Choquet-Bruhat [33] In modern language, they can be formulated as follows 

Theorem B.2. Let (T,,g,K) be as in the statement of the above theorem. Then 
there exists a smooth development (A4, g) of initial data. 

Theorem B.3. Let Ad, A4 be two smooth developments of initial data. Then 
there exists a third development A4' and isometric embeddings j : A4' A4, 
j : M' ^ M commuting with i, i. 

Application of Zorn's lemma, the above two theorems and simple facts about 
Lorentzian causality yields: 

Theorem B.4. ( Choquet-Bruhat- Geroch f35f) Let {J^,g,K) denote a smooth 
vacuum initial data set with cosmological constant A. Then there exists a unique 
development of initial data (M^g) satisfying the following maximality state- 
ment: If (M.^'g) satisfies (0) with embedding i, then there exists an isomet- 
ric embedding j : A4 ->■ M. such that j commutes with i. 

The spacetime {M..,g) is known as the maximal development of {T,,g,K). 
The spacetime A4 n J'*'(T,) is known as the maximal future development and 
A4 n (T,) the maximal past development. 

We have formulated the above theorems in the class of smooth initial data. 
They are of course proven in classes of finite regularity. There has been much 
recent work in proving a version of Theorem IB. 21 under minimal regularity as- 
sumptions. The current state of the art requires only g e K e H^^'^. 
See P^ . 

We leave as an exercise formulating the analogue of Theorem IB. 41 for the 
Einstein-scalar field system (|108p , (|112p , (|113p , where the notion of initial data 
set is that given in Section IB. 2. 21 



B.4 Harmonic coordinates and the proof of local existence 

The statements of Theorems IB. 21 and IB. 31 are coordinate independent. Their 
proofs, however, require fixing a gauge which determines the form of the metric 
functions in coordinates from initial data. The classic gauge is the so-called 
harmonic gaugcPI. Here the coordinates are required to satisfy 

□gx^ = 0. (114) 

^^Then called Foures-Bruhat. 
^^also known as wave coordinates 
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Equivalently, this gauge is characterized by the condition 

5'"'r^. = o. (115) 

A hnearised version of these coordinates was used by Einstein |74] to predict 
gravitational waves. It appears that de Bonder [7D] was the first to consider 
harmonic coordinates in general. These coordinates are discussed extensively in 
the book of Fock |H1]. 

The result of Theorem |BT3] actually predates Theorem IB. 21 and in some form 
was first proven by Stellmacher [143] . Given two developments {Ai,g), {M,g) 
one constructs for each harmonic coordinates x^, adapted to S, such that 
= ff/ji'j clxQuf = dxg^v along E. In these coordinates, the Einstein vaccum 
equations can be expressed as 

^,r = QZxtr9^9^9^''dp9''^ (116) 

for which uniqueness follows from general results of Schauder |136| . This theo- 
rem gives in addition a domain of dependence propertvF^ 

Existence for solutions of the system (|116|) with smooth initial data would 
also follow from the results of Schauder |136| . This does not immediately yield 
a proof of Theorem IB.2( because one does not have a priori the spacetime 
metric g so as to impose (I114p or (I115|) ! The crucial observation is that if 
(IllSp is true "to first order" on S, and g is defined to be the unique solution 
to dnH), then dnSI) will hold, and thus, g will solve (fTU5)) . Thus, to prove 
Theorem IB.21 it suffices to show that one can arrange for (jllSp to be true 
"to first order" initially. Choquet-Bruhat (33] showed that this can be done 
precisely when the constraint equations (|110p - (|llip are satisfied with vanishing 
right hand side. Interestingly, to obtain existence for (|116p . Choquet-Bruhat 's 
proof [33] does not in fact appeal to the techniques of Schauder |136| . but, 
following Sobolev, rests on a Kirchhoff formula representation of the solution. 
Recently, new representations of this type have found applications to refined 
extension criteria [103j . 

An interesting feature of the classical existence and uniqueness proofs is 
that Theorem IB. 31 requires more regularity than Theorem IB. 21 This is because 
solutions of (|114p are a priori only as regular as the metric. This difficulty has 
recently been overcome in jl29| . 

B.5 Stability of Minkowski space 

The most celebrated global result on the Einstein equations is the stability 
of Minkowski space, first proven in monumental work of Christodoulou and 
Klainerman |51| : 

There is even earlier work on uniqueness in the analytic category going back to Hilbert, 
appealing to Cauchy-Kovalevskaya. Unfortunately, nature is not analytic; in particular, one 
cannot infer the domain of dependence property from those considerations. 
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Theorem B.5. Let (T,,g,K) be a strongly asymptotically flat vacuum initial 
data set, assumed sufficiently close to Minkowski space in a weighted sense. 
Then the maximal development is geodesically complete, and the spacetime ap- 
proaches Minkowski space (with quantitative decay rates) in all directions. More- 
over, a complete future null infinity can he attached to the spacetime such 
that J-{I^)=M. 

The above theorem also allows one to rigorously define the laws of gravita- 
tional radiation. These laws are nonlinear even at infinity. Theorem IB. 51 led to 
the discovery of Christodoulou's memory effect |31] . 

A new proof of a version of stability of Minkowski space using harmonic 
coordinates has been given in |113| . This has now been extended in various 
directions in [M]. The original result [HI] was extended to the Maxwell case in 
the Ph.D. thesis of Zipser |157| . Bieri [13] has very recently given a proof of a 
version of stability of Minkowski space under weak asymptotics and regularity 
assumptions, following the basic setup of [5T] . 

There was an earlier semi-global result of Friedrich [53] where initial data 
were prescribed on a hyperboloidal initial hypersurface meeting . 

A common misconception is that it is the positivity of mass which is somehow 
responsible for the stability of Minkowski space. The results of [113j for this are 
very telling, for they apply not only to the Einstein-vacuum equations, but also 
to the Einstein-scalar field system of Section IB. 2. 21 including the case where the 
definition of the energy-momentum tensor (|113p is replaced with its negative. 
Minkowski space is then not even a local minimizcr for the mass functional in the 
class of perturbations allowed! Nonetheless, by the results of |113| . Minkowski 
space is still stable in this context. 

Another point which cannot be overemphasized: It is essential that the 
smallness in (IB.5|) concern a weighted norm. Compare with the results of Sec- 
tion [Ml 

Stability of Minkowski space is the only truly global result on the maximal 
development which has been obtained for asymptotically flat initial data without 
symmetry. There are a number of important results applicable in cosmologi- 
cal settings, due to Friedrich [55], Andersson-Moncrief [3], and most recently 
Ringstrom [134j . 

Other than this, our current global understanding of solutions to the Ein- 
stein equations (in particular all work on the cosmic censorship conjectures) has 
been confined to solutions under symmetry. We have given many such references 
in the asymptotically flat setting in the course of Section [2] The cosmological 
setting is beyond the scope of these notes, but we refer the reader to the re- 
cent review article and book of Rendall |132[ 1133] for an overview and many 
references. 

C The divergence theorem 

Let {A4,g) be a spacetime, and let Sq, Si be homologous spacelike hypersur- 
faces with common boundary, boimding a spacetime region B, with Si c J^(So). 
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Let rig , denote the future unit normals of Eq, Si respectively, and let de- 
note a one-form. Under our convention on the signature, the divergence theorem 
takes the form 

L p,<+ h'p,= L p,<^ (117) 

J z^i Jo J ^0 

where all integrals are with respect to the induced volume form. 
This is defined as follows. The volume form of spacetime is 

V^^detgda;" ... da;" 

where det g denotes the determinant of the matrix gap in the above coordinates. 
The induced volume form of a spacclikc hypcrsurfacc is defined as in Ricmannian 
geometry. 

Wc will also consider the case where (part of) Si is null. Then, we choose 
arbitrarily a future directcdnuU generator nf for Si arbitrarily and define the 
volume element so that the divergence theorem applies. For instance the diver- 
gence theorem in the region TZ{t' , t") (described in the lectures) for an arbitrary 
current then takes the form 

where the volume elements are as described. 

Note how the form of this theorem can change depending on sign conventions 
regarding the directions of the normal, the definition of the divergence and the 
signature of the metric. 

D Vector field multipliers and their currents 

Let be a solution of 

□ gV=0 (118) 

on a Lorentzian manifold (A^, g). Define 

Tt,.{^) = d^,i,dA - ig^,a"V9aV' (119) 

We call r^i, the energy-momentum tensor of i/'IlB Note the symmetry property 

T = T 

The wave equation (jllSp implies 

V^T^. = 0. (120) 
Given a vector field F^, we may define the associated currents 

Jli^) = V^T.M) (121) 

^^Note that this is the same expression that appears on the right hand side of (I108II in the 
Einstein-scalar field system. See Section IB. 2. 21 
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K"" = ^TT^.T'^^CV') (122) 
where tt^ is the deformation tensor defined by 

The identity ([HP]) gives 

Note that Jj^i'tp) and K^{tp) both depend only on the 1-jet of tp, yet the 
latter is the divergence of the former. Applying the divergence theorem (|117p . 
this allows one to relate quantities of the same order. 

The existence of a tensor T^,y('0) satisfying (|120p follows from the fact that 
equation (IllSp derives from a Lagrangian of a specific type. These issues were 
first systematically studied by Nocther [124] . For more general such Lagrangian 
theories, two currents J^, K with V^J^i = K, both depending only on the 1-jet, 
but not necessarily arising from T^^ as above, are known as compatible currents. 
These have been introduced and classified by Christodoulou [JB]- 

E Vector field commutators 

Proposition E.0.1. Let tjj be a solution of the equation of the scalar equation 

and X be a vectorfield. Then 

HgiX^) = X{f) - 2 ^7r"^VaV/3V - 2 (2(V" ^7r„ J - (V,, ^O) VV- 
Proof. To show this we write 

X(D<,V) = /:x(5"^VaV/3V) = 2^^"''VaV/3^ + ff^'^/Ix ( Va V/3 V) • 

Furthermore, 
and 

□ 

F Some useful Schwarzschild computations 

In this section, {M.,g) refers to maximal Schwarzschild with M > 0, Q = 
M/S0{3), I^, J^{T^) are as defined in Section [SH 
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F.l Schwarzschild coordinates (r, t) 

The coordinates are (r, t) and the metric takes the form 

-(1 - 2M/r)dt^ + (1 - 2M/ry^dr^ + r^dcrgs 

These coordinates can be used to cover any of the four connected components 
of Q \ Ti"^. In particular, the region J^{I\) n J^(I^) (where 1% correspond to 
a pair of connected components of X"^ sharing a hmit point in the embedding) is 
covered by a Schwarzschild coordinate system where 2M <r<oo, -oo<t<oo. 
Note that r has an invariant characterization namely r(a;) = y'Area(S)/47r where 
S is the unique group orbit of the S0(3) action containing xW\ 

The hypersurface {t = c} in the Schwarzschild coordinate region J^{I^) n 
J'*'(I^) extends regularly to a hypersurface with boundary in A4 where the 
boundary is precisely n H^. 

The coordinate vector field dt is Killing (and extends to the globally defined 
Killing field T). 

In a slight abuse of notation, we will often extend Schwarzschild coordinate 
notation to V, the closure of J"(I^) n J^(I^). For instance, we may talk of 
the vector field dt "on" H"^, or of {t = c} having boundary "H^ n H^, etc. 

F.2 Regge- Wheeler coordinates {r*,t) 

Here t is as before and 

r* = r + 2Af log(r - 2M) - 3M - 2Af logM (123) 

and the metric takes the form 

-{1 - 2M /r){-dt'^ + {dr*f) + r'^das2 

where r is defined implictly by ()123h . A coordinate chart defined in -oo < r* < 
oo, -oo <t < oo covers J~(IJ) n J+(I^). 

The constant renormalisation of the coordinate is taken so that r* = at the 
photon sphere, where r = 3Af. 

Note the explicit form of the wave operator 

□g^ = -(1 - 2M/r)-\dfij - r-^dr*ir^dr*^)) + f^fA^ 

where ^ denotes the induced covariant derivative on the group orbit spheres. 

Similar warnings of abuse of notation apply, for instance, we may write 
dt = dr* on v.'*'. 

F.3 Double null coordinates {u,v) 
Our convention is to define 

u=-it-r*), 
2 

Compare with the Minkowski case Af = where the SO(3) action is of course not unique. 
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v=-(t + r*). 

The metric takes the form 

-4(1 - 2M/r)dudv + r^da§,-2 

and J^{I\) n J^(I^) is covered by a chart -oo < u < oo, -oo <v < oo. 

The usual comments about abuse of notation hold, in particular, we may now 
parametrize Ti^ n V with {00} x [-00, 00) and similarly T-L' n V with (-00, 00] x 
{-00}, and write dv{-oo,v) = dt{-oo,v), du{-°°,v) = 0. 

Note that the vector field (1 - 2Af/r)~^9„ extends to a regular vector null 
field across T-L^ \ T-T . Thus, with the basis 9u, (1 - 2M jrY^du, one can choose 
regular vector fields near T-L^ \ without changing to regular coordinates. In 
practice, this can be convenient. 
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